For Referetice 


NOT TO BE TAKEN FROM THIS ROOM 


Gx apnis 
UNITASTTATIS 
AUBERTAEASIS 


AAT 


8 Wn om 
ave at e ¢ 
MAY hota a 
a 1 x i 
Wey ‘ / 
Tale Mee vi i 
ad va a | ep 
Nana? 
, ive 
4 h 
A Hee FOS ts YH 
my ay “i 
f 4 m 
1 7 
; 
‘ j f 
Y a 
, oo [ hs 
fa ae: 
i J Aa Tt, 
4 1 a 
i ¢ sy i 
. ti re Sip? 
Ni eae he WL é a Mi 
Wy i e AVAL? 
‘ i y ¥ 
! Li i eee 
r = Pa 
v Le 
i AG 
a ‘ Las 4 
Wy eee! y " 
ea 
: i 
d avin 
{ al } f } 
i ‘ Ay ‘ 
ian ; ; 
i 
en) A 
i , LF aa 8 
- 1 i ® 1% 
‘ h aa ; 
r Whe, a 
DN aan 
H 
: ‘ / v 1 : 
A t | 
; oh 
7 N 
R 1! We) 
U ‘ . 
nr i at Oi 
f a , 
os ANA 
i a 
“ ;  ) 
; r 


a | 
Sh 
r 
oh, 
ee 
»> ¥ 
¥ 
\ 
f i a 
Vlg 5 4 
aA fy 
a] , 
A a 


Ly 
i x 
ee 


vis 


2 


THE UNIVERSITY OF ALBERTA 


APPLICATION OF FUNCTIONAL ANALYSIS OPTIMIZATION TECHNIQUES 
TO THE ECONOMIC SCHEDULING OF HYDRO-THERMAL ELECTRIC 
POWER SYSTEMS 


by 
Pi 
G, MOHAMED EL-AREF EL-HAWARY 


A THESIS 
SUBMITTED TO THE FACULTY OF GRADUATE STUDIES 
AND RESEARCH IN PARTIAL FULFILMENT OF THE 
REQUIREMENTS FOR THE DEGREE OF 
DOCTOR OF PHILOSOPHY. 


DEPARTMENT OF ELECTRICAL ENGINEERING 


EDMONTON, ALBERTA 
PAUL 97 ¢ 


ABSTRACT 


In this thesis problems of optimum scheduling of hydro-thermal 
power systems are discussed. The scheduling problems are solved by 
the use of functional analysis, and in this case, the minimum norm 
formulation is employed. 

Optimal schedules are derived for the classical all-thermal 
power system problem. Here the scheduling equations are shown to be 
equivalent to earlier results obtained using other optimization 
techniques. The problem of hydro-thermal systems scheduling when the 
hydro plants are located on separate streams is treated. One, a system 
with fixed head hydro-plants and negligible transmission loss is 
considered. Two, a fixed head hydro plants system where transmission 
losses are included is presented. Three, a power system is considered 
where head variations at the hydro-plants are not negligible. In each 
case the optimal schedule is obtained. The problem of actually 
implementing the optimal schedules is discussed. The scheduling 
equations obtained are shown to be equivalent to those obtained using 
other techniques. 

New scheduling equations are developed for the common-flow hydro- 
plants of the system. Here the time taken by water to flow from the 
upstream plant to the downstream plant is taken into account. Also, 
the tail-race elevation at the hydro-plants is considered. 

In the above problems the general loss formula is employed to 


represent the electric network of the system. However, in the final 
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chapter of this thesis, the problem of the optimal hydro-thermal load 
flow is investigated. Here the exact model of the electric network 
is employed. Reliability considerations are also incorporated in the 
formulation. Finally, it is shown in the last chapter how the 
efficiency variations and trapezoidal reservoirs effects can be 
included in the formulation. 

The computational aspects of the obtained scheduling equations 
are discussed. Practical examples are also given to illustrate the 


results obtained. 
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CHAPTER I 


INTRODUCTION 


1.1 Background 


A prime objective in the operation of a power system is to achieve 
optimum economic dispatch. This is the problem of scheduling the 
generation at various generating plants. Here the system's power demand 
is to be supplied at the lowest possible power production cost. Normally 
this is planned for an appropriate time interval. The optimum operation 
of a power system will also depend upon restrictions imposed by factors 
other than operating economics. Possible decreases in power production 
costs of only fractions of a per cent are still of vital concern to the 
electric utility industry. In addition, the capability to solve the 
economic dispatch problem is extremely useful for the planning and design 
of future equipment additions to power systems. For these reasons, the 
economic dispatch problem has been the subject of extensive research [1,2]. 

In economic dispatch it is customary to consider the operating costs 
only. This ignores expenses of capital, labour, start-up and shut-down 
related to the length of the outage period for a certain unit. It is 
essential to have an accurate knowledge of the manner in which the total 
cost of operation of each available energy source varies with the instant- 
aneous output. In most cases, the fuel supply available to the conventional 
thermal plant is not a limiting factor in the operation of the plant. In 


such cases the appropriate price to use for economic dispatch is the 
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current cost of incoming fuel adjusted for handling costs, maintainance 
cost of fuel handling facilities, etc. ‘In this thesis it is assumed 
that the fuel cost is a second order polynomial of the power output [3]. 

The hydro-thermal optimization problem is different from the all- 
thermal one. The former involves the planning of the usage of a limited 
resource over a period of time. The resource is the water available for 
hydro-generation. Most of the hydro-electric plants are multipurpose in 
nature. In such cases it is necessary to meet certain obligations other 
than power generation. These may include a maximum forebay elevation 
not to be exceeded due to flood prospects and a minimum plant discharge 
and spillage to meet irrigational and navigational commitments. Thus 
the optimum operation of the hydro-thermal system depends upon the 
conditions which exist over the entire optimization interval [4]. Many 
systems with large water storage capacity will require a year for the 
optimization interval. Another system may have run-of-the-river plants 
with only a small or moderate storage capacity. An optimization interval 
of a day or a week may be useful in this case [5]. 

Other distinctions among power systems are the number of hydro 
Stations, their location and special operating characteristics. The 
problem is quite different if the hydro stations are located on the same 
stream or on different ones. In the former case, the water transport 
delay may be of great importance [6, 7, 8]. An upstream station will 
highly influence the operation of the next downstream station. The 
latter, however, also influences the upstream plant by its effect on the 
tail water elevation and effective head. Close coupling of stations by 


such a phenomenon is a complicating factor. 
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1.2 Scope of the Thesis 


In this thesis, optimum generation schedules will be developed 
for a hydro-thermal electric power system. The scheduling problem is 
solved by use of functional analysis where the minimum norm formulation 
is employed. The solution found here is guaranteed to be the unique 
optimal solution [9]. The power system considered contains an arbitrary 
number of plants. 

Aspects of functional analysis optimization techniques that will 
be applied to the power system problems are presented in Chapter II. In 
this presentation it is chosen to start with the simplest possible problem 
and then gradually lead up to the more complex problems. 

In Chapter III, we consider an all-thermal electric power system 
with negligible transmission losses. This problem is a classical 
one [10]. It is included here to show the application of the minimum 
norm formulation in the case of additive linear transformations. 

Chapter IV is concerned with hydro-thermal system with hydro plants 
on separate streams. The electric network representation relies on the 
active power balance equation. Section 4.1 covers the background and 
previous work done in this area using other optimization techniques. 
Attention is focused in this chapter on the hydraulic head variations in 
the case of vertical sided reservoirs and transmission losses. In 
Section 4.2 a system with fixed head hydro plants and negligible 
transmission losses is considered. The assumption of negligible trans- 
mission losses is relaxed in Section 4.3. Finally, the more general 
problem of variable head hydro plants with transmission losses in the 


active power balance equation is presented in Section 4.4. The results 
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obtained are shown to agree with Kron's equations [11]. The practical 
implementation of actually finding the optimum generation schedules is 
illustrated by way of two examples in sub-sections 4.3.5 and 4.4.6. 

Here functional analytic computational search techniques are employed. 

In Chapter V the problems of common-flow between plants are 
considered. Here the time delay of flow between hydro plants on the 
same stream is taken into account. The effect of tail water elevation 
on the effective head at the hydro plants is considered. The problem 
of the power system with variable head hydro plants on the same stream 
is the subject of Section 5.2. Here a practical example showing the 
computational aspect of the problem is given. The computational scheme 
employed here is again a functional analytic one. A general power 
system in as far as the relative locations of the hydro plants are 
concerned is treated in Section 5.3. 

A more complex problem in the area of economic scheduling of power 
systems is that of optimal load flows. This is the subject of Chapter 
VI. Here a more realistic representation of the electric network is 
adopted. Starting from the basic model, a suitable form of the load flow 
equation is derived to facilitate the mathematical formulation. All of 
the electric variables of the system are considered in this chapter. 
Here we also include the more ambitious objective of system reliability. 
Also inequality constraints imposed on the network variables are con- 
sidered. Furthermore, consideration is given to variable efficiency and 
trapezoidal reservoirs hydro plants at the end of this chapter. 

Only the titles and the broad outline of the problems considered 


in this thesis are mentioned here. A more detailed description of each 
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of the problems and its relationship to the previous work in this area, 


will be found at the beginning of each chapter. 


Lani, 
ry i id r wa 
he BNA at t B; Ln 


7 


r % a ‘ ZY ie a 


CHAPTER II 


THE FUNCTIONAL ANALYTIC OPTIMIZATION TECHNIQUE 


2.1 Background 


During the years optimal control theory was developed, powerful 
general solution methods were introduced. These are based on the now 
widely known "Maximum Principle" and "Optimality Principle". Parallel 
to the development of these, starting in 1956, attempts have been made 
to introduce methods of functional analysis into the study of optimal 
control problems. 

At first it seemed that the methods of functional analysis 
applied only to a very restricted class of problems. But in spite of 
this, the number of studies using the ideas of functional analysis has 
increased. In solving optimal control problems, by using the Maximum 
Principle, or by reduction to the Euler equations, these methods do not 
show how to select the initial conditions required for solving the 
adjoint system. The methods of Dynamic programming and the approach 
that leads to the Hamilton-Jacobi equations do not have this deficiency. 
However, the solution of functional equations is not an easy problem [13]. 

One of the typical features of the functional analysis approach 
is that it yields necessary and sufficient conditions for the existence 
of solutions. This fact makes it possible to study the qualitative 
aspects of optimal processes. Moreover, this approach is free of the 
concrete nature of the system. Thus many formulations hold for systems 
that are distributive, digital, composite, nonlinear or biological. Of 


course, results obtained on the basis of an abstract formulation must 
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then be given concrete identification in its various physical forms. 

Below, we give a survey of certain works where the methods of 
functional analysis are used in solving problems in the theory of 
optimal processes. It is not the author's intention to give a complete 
exposition of the application of functional analysis to the theory of 
optimal processes. In fact, works dealing with the abstract minimum 
norm formulation will be our main concern. 

Investigation of the problem of approximate solutions to first 
order ordinary differential equations, led D.S. Carter [14], in 1957, 
to the earliest minimum norm formulation. Carter's problem was 
concerned with obtaining an element of a specific Banach space. The 
image of the element sought under a first order linear differential 
operator was to be of minimum norm. The norm adopted was the maximum 
norm. The element was to satisfy a two-point-boundary condition. 

In 1962, W.T. Reid [15] extended Carter's results to the case of 
an nth order differential operator. This was achieved by reducing the 
problem to a problem in the theory of moments. The general results of 
the Hahn-Banach theorem were then applied to the reduced problem. 

A minimum norm problem in Hilbert spaces was considered by A.V. 
Balakrishnan [16] in 1962: Given a compact, bounded, linear operator 
mapping Hy into Ho and with Hy and Ho being Hilbert spaces an element 
of a given ball in Hy is sought. The norm of the difference between the 
image of this element and a given element in Ho is to be a minimum. A 
sequence of elements in Ho was obtained and was shown to converge to the 
desired element. 

A wide class of minimum norm problems was considered in 1962 by 


L.W. Neustadt. In his paper [17], Neustadt employed a variational 
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approach to find the minimum norm element. This was set - the Banach 
spaces of the L, type (1 < p<). The system satisfied a linear 
integral operator and the problem was reduced to minimizing a functional 
of a new variable. This variable is in many ways analogous to the co- 
states of the Pontryagin's Maximum Principle. The method of steepest- 
descent was suggested for implementing the final optimal. 

G.M. Krane and P.E. Sarachick [18] considered a minimum norm 
problem in 1963. This was essentially the same as Neustadt's problem. 
The only exception was that the element (control) sought was to belong 
to a specified ball in the Ly Space under consideration. Holder's 
inequality was used extensively to specify the optimal solution. 

While most authors chose to consider relatively well specified 
systems, W.A. Porter chose a more abstract approach to the optimization 
problem. Knowing that diverse systems of equations can be associated 
with linear transformations, Porter (1964) considered a problem 
involving a linear transformation on a "Hilbert space". The cost 
associated with an element of the Hilbert space was given by the Hilbert 
space norm [19]. Later together with J.P. Williams [20, 21] he extended 
the results of this abstract problem to cases involving Banach spaces. 

It was noted that the results of these approaches are applicable 
to systems of discrete, continuous and composite types. These results 
can be utilized for various optimization problems [22, 23, 24]. The 


work of this dissertation involves a recognition of some of these areas. 


2.2 Some Functional Analysis Concepts 


The discussion of this section is aimed at displaying some of the 
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concepts and symbols which are utilized in the next See ar Basic 
functional analysis concepts are left to the references of which 
Poe425..co 2/5 25] arewel lstnted. 

The next section will be concerned with linear spaces which have 
norms defined upon them. A norm (commonly denoted by ||.||) is a real 
valued, positive definite (||x|| > 0 for x # 0), absolutely homogeneous 
({[ax|| = [al-||x]]), and subadditive (||x + y|| < ||x]]+]]y]|) function. 

A transformation is a mapping from one vector space to another. 

If T maps the space X into Y, we write T: X+Y. If T maps the vector 

XeX into the vector yeY, we write y = T(x) and y is referred to as the 
image of x under T. Alternatively a transformation is referred to as an 
operator. 

The transformation T: X+Y is said to be linear if for every x,, 

Xo eX and all scalars 04 and a, one has T(a4X, + WnXq) = a4 T(x, ) + aT (X5). 
The linear operator T from a normed space X to a normed space Y is said 

to be bounded if there is a constant M such that ||Tx|| < M||x|| for all 
XeX. The normed space of all bounded linear operators from the normed 
space X into the normed space Y is denoted by B(X,Y). 

A functional is a transformation from a linear space into the space 
of real (or complex) scalars. A functional f on a linear space X is 
linear if for any two vectors x, yeX and any two scalars a and g there 
holds f(ax + By) = af(x) + Bf(y). A linear functional f on a normed 
space is bounded if there is a constant M such that |f(x)| < M||x|| for 
every xeX. The smallest such constant M) is called the norm of f. The 


norm of the functional f can be expressed as 
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[|] | = sup| f(x) |/| [x] | | (2.1) 
x #0 


Given a normed linear space X, one can define bounded linear 
functionals on X. The space of these linear functionals is a normed 
linear space x". The Space X" is the normed dual of X (alternatively 
x" is called the conjugate space of X), and is a Banach space. If 
X is a Hilbert space, then X = x" Thus Hilbert spaces are self-dual. 
The normed dual as of x" is called the second dual space of X. A 
normed linear space X is said to be reflexive if X = qiie Any Hilbert 
space is reflexive. 

Let X and Y be normed spaces and let TeB(X,Y). The adjoint 


e * * * e ° 
(conjugate) operator T : Y +X is defined by 
* * 
<Kyiuy> = <1X,y > 


An important special case is that of a linear operator T: H+G where H 
and G are Hilbert spaces. If G and H are real then they are their own 
duals and the operator Ta can be regarded as mapping G into H. In this 


case the adjoint relation becomes 
* 
<5 Vomn eX, ley 


Let G and H be Hilbert spaces and TeB(G,H) with the range of T 


being closed. Define the set M as M = {x eG: ||Tx, - y||=min||Tx - y|[}. 
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Let x eM be the unique vector of minimum norm. Then the pseudo-inverse 
Operator To of T is the operator mapping y into its corresponding Xo 
as y varies over H. 

The Minkowski functional p of a convex set K in a normed linear 


Space X is defined on X by: 
p(x) = inf{r: Sek r>0} C2n2) 


In a Banach space B, the set U = {xeB: ||x|| < 1} is said to be 
the Closed Unit Ball in B. The boundary aU of U is called the unit 
sphere aU = {xeB: ||x|| = 1}. The Banach space B is called rotund 
(strictly convex) if 3U contains no line segments. Any convex set K in 
a rotund space has at most one minimum element. The Banach space B is 
said to be smooth if at each point of oU there is exactly one supporting 
hyper-plane of U. 

A vector xeU is called an extremal of feB. if x satisfies fOxe= |) fl]. 
At most one extremal of f exists if and only if B is rotund. If xeB then 
feU is an extremal of x if f satisfies f(x) = ||x||. At most one 
extremal of x exists if and only if B is smooth. If xeB, the Hahn-Banach 
theorem [9,28] guarantees the existence of at least one extremal f in Be 
Accordingly, if B is reflexive, every f # 0 and feB. has a unique 
extremal in B if and only if B is rotund. Every x # 0 and xeB has a 
unique extremal in B. if and only if B™ 1s,.smoo thie. inadl Ly 4 <i Bess 
reflexive, rotund and smooth and we denote by x the unique extremal of 


x in Be then x = TIT for all x # 0 and xeB. 
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2:3) The Minimum Norm Problem 


The main results of this dissertation are based on previous 
analysis by Porter and Williams [19, 20, 21] of an abstract 
minimum norm problem. For background purposes this section considers 


this problem which may be formulated as follows: 


Let B and D be Banach spaces. Let T bea 
bounded linear transformatton defined on 
B wtth values in D. For each & in the 
range of T, find an element weB that 
satisfies 

Seog! 
whtle minimizing the performance tndex 


Tay = rel | 


The solution of this problem as obtained by Porter and Williams 
will be given in the form of theorems. 

Theorem 1 :(Existence 

The minimum norm problem formulated above has a solution for 
every bounded linear transformation defined on B if and only if B is 
reflexive. 

Theorem 2 :(Uniqueness 

For every gD, a unique optimal solution exists if B is reflexive, 
rotund and smooth. 


Theorem 3 : (Characterization) 


With the conditions of theorems 1 and 2 satisfied there exists 
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° ° * ; 
a unique unit ¢. eD for every gD which defines the optimal u denoted 


4 
by ur, as: 
aot ee 
up = TE = ple)T 6, eyes) 
o, is determined by the conditions 
2. Either: 


a. <nsb,> < [p()]'<,9,> for all neC 
or 
b. [IT oel [p(e) 1" '<e,6,> 
Theorem 3 defines the optimal solution in an implicit form. To 
eliminate the ambiguity the following transformations are defined: 
1. Define the transformation K: Bak by the relation 
K(f) = ||f||f for every feB- (2.4) 
Thus K is an extremizing norm restoring operator. 


2. Define the transformation J: DD by the relation 


J(g) = ||T'o||T(T's) for every oeD™ (2.5) 
Thus 
J(¢) = T(K(T'6)) 6) 
or 
J = TKT” (2.7) 


Thus theorem 3 can be restated as 


Theorem 3 :(Modified Optimal Characterization 


The unique optimal u eB is given by: 


where the pseudo-inverse operator TS is given by 
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(229) 


In Hilbert spaces the operator K reduces to the identity operator 


so that 
Te =T oT iz 


provided that the inverse of Tp exists. 


(2.0) 


An extension of the results of the minimum norm problem discussed 


in this section is as follows: 


Let B, D, T and & be as in the minimum 
norm problem. Let a be a gtven vector 


tn B. Then the unique u,cB satisfying 


E 
E.= fu 

whteh minimtzes the performance tndex: 
J(u) = ||u - a | 


ts given by 


ue = r'[e - Tal+ a 


(217) 
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CHAPTER III 


A THERMAL POWER SYSTEM 


3.1 Background 


This chapter is devoted to the problem of optimizing the operation 
of an exclusively thermal electric power system. Much of the early 
investigations in the economy loading field were centred upon steam- 
electric generating units operating in parallel. This resulted in 
the development of the theory of incremental loading [10, 29, 30]. The 
method leading to this theory successfully employed simple calculus 
techniques. 

Achieving minimum total cost of supplying the power requirements 
of a system requires an accurate knowledge of the operating cost 
functions. A cost function represents the manner in which the total 
cost of operation of a generating unit varies with its output. The total 
cost of operation includes the fuel cost, cost of labour, supplies and 
maintenance. However, no methods are presently available for expressing 
the latter as a function of the output [1, 10, 31]. Arbitrary methods 
of determining these costs are used. The most common one is to assume 
the cost of labour, supplies and maintenance to be a fixed percentage of 
the incoming fuel costs. It is common practice to obtain the operating 
cost function by establishing the input-output curve for the plant 
considered, then adjusting for the cost of the fuel per unit input to 


the plant. Throughout this dissertation it is assumed that the cost 


my 


Ce poe 


TT a ah ae eae we uD ae ae Be A i 
i eG pt ae vee eo * oe TG oe : sin ' \) Deny \ hike a OF 
sine ae i 
. Teva > ae SNE a ce Ae | 
oe 5 a hg i + Y ri , " vy sy i” % up 
5 f \ % a § s on " ta an 7 r \ 
fi = zt ) of reas AN 
5 2 > ? 4 ' 7 i an 
} + { ay eae Md i“ 


notte iego sit ontsimt3qo to natonge ai ig iy ne a 


a ae an 
f Tne) | iin 


; yii69 oft to dau © bmateye a 1awoq Stats) Fenris ‘Uovtey 
~(ipede noqu bevinea ew bist? on tbeol wronoa8 ay nt 


at bed }fueen atm Riyciliis BE oni dnyege. atin vents ‘ 


ip 


pee 2 sfgat2 tovotane pteasoaue oom ahh of at 
Ah ene Pits \ 
Diy 


Tetod aia often eninge “ont, ‘einsebiqay 


‘is3ot ont -thiggio att. a Ai 2aiysy. Fin pate 
“ba ‘esifqgue cwodet to $200. $209 lout arth 


“entezayaxs 708 aldetteve mazar, ste 2bo,tjan or “on 


abortson Soaysieaa te OL aa) tugtue att +0 poh 
at — | 


Tame he 9 


ee | emiees of 2h 80 omg! $208 ‘oT sbeau aie 


*. 


weer vont és a sonanedat on bis 


git 


‘ i 


i a ie y 
y n “ r in af pniasiieg 
j a Ua : 


9 Here a ne om 
tomule2 abet wok 59% zeth 


Se 
a. Tana ie 


= 
an 


functions may be approximated by second order polynomials of the 
instantaneous output power of the thermal plants [3, 32]. 

The problem considered in this chapter deals with the case when 
all sources of generation are either located at the same bus or close 
enough so that transmission losses may be neglected. In this case the 
instantaneous system power demand is equal to the sum of all available 
instantaneous generations. Accurate forecasting of the hourly power 
demand is thus an important factor in any dispatch strategy. This 
depends on both weather information and historical patterns. Procedures 
for active power demand forecasting are in existence [33]. In this 
thesis it is assumed that the daily power demand of the system is known 


beforehand. 


3.2 Statement of the Problem 

Consider a power system with m thermal plants on the same bus 

supplying the load centre. Assume that the system's power demand is a 
known function of time over the optimization interval. The problem is 
to determine the generation of each plant as a function of time over the 
optimization interval under the following conditions: 
| 1. The total operating cost over the optimization interval is a 
minimum. 

2. All operating costs attributed to the fuel cost at the thermal 
plants (adjusted for labour, supplies and maintenance) are 
given by: 

at 
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3. The total instantaneous power generation in the system matches 


the power demand. 


P(t) = d, ira 08 2) 


o 
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3.3 Mathematical Formulation 


The object of the optimizing computation is: 


wes 3 


_ Fils, (t))et (3.3) 
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subject to the constraint given by (3.2). Define the m dimensional 
control vector function u(t) by: 


Ve Reag ena (3.4) 
le m 
Define next the m dimensional vectors L and M as: 


u(t) = col.[P, (t), P. 


iGo 


= col.[8,, Bos sees Bad (2.5) 
COL Pils Stixcs al tn) (3'6) 


= 
T 


Finally define the mxm diagonal matrix B by: 
Bre diagLy,> Yo? sees ea (307) 
Substituting (3.1) i, @3-3) fandarsing (3.4), (3.5) and (3.7) 
transform the cost functional of (3.3) to: 


i 
J, (u(t)) Hh L'u(t)dt off u(t) B u(t)dt (3.8) 


Note that the a. 'S of (3.1) are dropped from (3.8), since these are 


independent of the control vector u(t). Furthermore, substitution of 


(3.4) and (306)/1n’(3.2)-yields: 
t) = Mu(t) (3.9) 
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In order to cast the problem in a minimum norm formulation, the 
control vector u(t) is taken as an element of the Hilbert space 


Lh" t0,T,]. This is the space of the m-dimensional vector valued square 


integrable functions defined on [O,T,]endowed with the inner product 


defined by: 


) 
f 

e¥(t), u(t)> = v(t) B u(t)dt (3.10) 
O 


ool tte ta te mt 8 


for every V(t) and u(t) in L 
positive definite for (3.10) to be a valid inner product definition. In 
the problem at hand, the y's are positive which ensures the validity of 
10). 

Furthermore, the given power demand function Pp(t) is considered 
an element of the space L,[0,T,] of square integrable functions defined 
on [0.T,], with the inner product definition: 


E 
we) say) t) > SGM (2a) 


for every x(t) and y(t) in L,[0,T,]. 
With these definitions it is seen that (3.9) defines a bounded 
linear transformation T: LS") ro, Te bL, [0.7 eI. Also the cost functional 


of (3.8) reduces to: 


Jo(ult)) = <Wou(t)> +] [u(t) || (3.12) 
where V is given by: 
vi = Lip (3.13) 


Furthermore, (3.12) can be written as: 
2 2 
Jo(u(t)) = |fu(t) + (W/2.)] [2-1 1v/2- 1] (3.14) 
Since V is independent of the control vector u(t), it is only needed to 


consider minimizing: 
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: (3.15) 


Jy(u(t)) = [Ju(t) + (w/2.)I| 
Also, since a norm is a positive scalar, then minimizing the norm is 
equivalent to minimizing its square. Thus the final cost functional 


that needs to be considered is given by 


J(u(t)) = |Ju(t) + (v/2.]| (3.16) 
subject to 

Py(t) = T(u(t)) (3.17) 
or 

T(u(t)) = M'u(t) (3.18) 


The problem is now formulated as a minimum norm problem whose form 
was given in section (2.3). The optimal solution vector u,(t) is 


obtained in the next section. 


3.4 The Optimal Solution 


In view of Section 2.3, the optimal solution to the problem 
formulated in the previous section is given by 
ue(t) = T'EPp(t) + T(W/2.)] - (W/2.) (3.19) 
The first step in characterizing the pseudo-inverse transformation tT is 


* 
to obtain the adjoint transformation T . This is obtained by utilizing 


the adjoint equation: 


<Pp(t),T(u(t))> = <T (Pp(t)) ult)» (3.20) 
The left-hand side inner product is evaluated in LoL0. 1] as 
if, 
Py (t) ,T(u(t))> =f Ppt) (a u(t) Jet (3.21) 
(9) 


The right-hand side inner product in Us" 300,T is given by: 


* 


T * 
<T™ (P(t) su(t)> ale (Py(t))]'B u(t)dt (3.22) 
8) 
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Thus the equality (3.20) reduces to 
T 


1 
* = 
J cogent - wTar"e, (ene u(tiat = 0 (3.23) 
0 
The equality (3.23) is satisfied for 


* = 
[T"(P)(t))] = M'B™'P, (t) (3.24) 
This yields an expression for the adjoint transformation i as 
T'[P)(t)] = B'M Pp(t) (3.25) 


The second step is to evaluate the operator J(Pp(t)) using the 


relation (2-7), and the equations (3.9), (3.17) and (3.24) to obtain 
3(Pp(t)) = MB 'M PL (t) (3.26) 
Note that (m'B7 M) is a quadratic form which has the scalar value 
m 
MiB M=) 3 (3.27) 
seen ee 
Thus the inverse operator J is obtained as: 
m 
a MP (t)) = Pity (3.28) 
1 DL YG 
Finally the pseudo-inverse operator Te is obtained as 
m 
ia bass ie 
T[P)(t)] = B M Pott /d, 7 (3.29) 


Note that substituting equations (3.6) and (3.7) in (3.29) yields: 


T'[P)(t)] = ee aca ans Te 

fh % tm Yi 
(3.30) 

The vector V in (3.19) is given by (3513)5 “Substituting. (3.5) 


and (3.7) the following is obtained: 
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Thus substituting (3.30), (3.31) and (3.32) in (3.19) the optimal 
vector given in (3.19) is obtained. Here recalling the definition of 
u, (t) in (3.4), the optimal power generations are obtained as: 


m m 
Pe ey EP 05D Bir Wry D (1/y;))} 


G4 1=] 


- HM eeieyel i) odsei lees asm (3.33) 


It is worth mentioning here that the optimal solution obtained 
here can be shown to satisfy the optimality conditions of the classical 


theory of variational calculus [34]. 
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CHAPTER IV 


POWER SYSTEMS WITH HYDRO PLANTS ON SEPARATE STREAMS 


4.1 Background 


In this chapter the problem of obtaining optimal generation 
schedules for hydro-thermal electric power systems is considered. The 
system considered here is characterized by the presence of hydro plants 
that are not on the same stream. According to their order of complexity, 
three optimal operation problems are formulated and solved in the next 
three sections. 

The optimization procedure is based upon predicted future system 
demand and a forecast of available water resources. Due to accuracy 
considerations present day forecasting techniques are only reliable for 
Short time intervals. Thus the deterministic optimization problems 
considered here are categorized as short-range economy dispatch problems. 

The inclusion of the system's transmission losses in the power 
balance equation is based on the well established loss formula [10]. It 
is noted that such a function as the transmission losses cannot be 
expressed in terms of only generator powers in an exact manner. There 
are a number of approximations involved in the loss formula derivation. 
However, it produces close answers with errors of up to a few percent [34]. 
Very sophisticated methods of calculating the loss formula coefficients 
exist and are being used [35, 36]. 


The operating costs of the hydro-electric plants do not change 
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with the power output. This would lead to an optimization schedule 
assigning the power demand to all existing hydro-plants. In practice, 
however, the amount of water available over the optimization interval 
is limited. Thus a constraint on the volume of water discharge over 
the optimization interval is imposed. 

A brief presentation of previous investigations of hydro-thermal 
power systems economic operation will be given. Various optimization 
techniques have been used in the past. Among these are the variational 
calculus principles, the methods of dynamic programming and the 
Pontryagin's Maximum Principle. 

In 1940, J. Ricard [37] obtained a set of operating schedules for 
a hydro-thermal system with no losses. His work was continued by W.G. 
Chandler, P.L. Dandeno, A.F. Glimm and L.K. Kirchmayer [38] in 1953, 
who included transmission losses but with constant hydraulic head. The 
latter method was improved in 1958 by A.F. Glimm and L.K. Kirchmayer [11] 
by including variable head plants. The work of Kron, who developed 
equivalent equations, was reported in [11]. 

A set of scheduling equations was developed in 1953 by R.J. Cypser 
[6]. These were developed under the assumption that variations in 
elevation and plant efficiencies can be neglected. J.J. Carey [39] 
suggested an approach which will linearize Cypser's equations. C.W. 
Watchorn gives in [4] a set of equations to be satisfied in order to 
achieve maximum economy. 

It should be noted that the above mentioned investigations employed 
the Euler equationof the calculus of variations to obtain the scheduling 


equations. 
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The work of R.A. Arismunander [40] in 1960 and ae 1962 [41] 
jointly with F. Noakes dealt with short-range optimization of hydro-thermal 
systems. He employed all the necessary and sufficient conditions for 
optimality of the variational calculus. In addition to this he proved 
the equivalence of all previously developed equations. 

In 1962, J.H. Drake, L.K. Kirchmayer, R.B. Mayall and H. Wood [7] 
presented a dispatch formula based on the calculus of variations. This 
formula is restricted to the case where all the hydro-plants operate 
with constant head. The system considered has series plants, multiple 
chains of plants and intermediate reservoirs. 

L.K. Kirchmayer and R.J. Ringlee presented a dispatch formula for 
a hydro-thermal power system in[42]. Head variations were considered. 

The formula applies for power systems having one hydro-plant. 

In a discussion of the paper in reference [7], C.W. Watchorn 
points out the importance of considering variable head for the optimization 
of such systems. In separate discussions of the same paper, C.W. Watchorn 
and R.A. Arismunandar, both point out that a river time delay of a couple 
of hours is highly important for accurate optimization of many power 
systems. 

B. Bernholtz and L.V. Graham [43] presented a dynamic programming 
solution to the hydro-thermal optimization problem. L.K. Kirchmayer and 
R.J. Ringlee [42] have also a presentation of some results and conclusions 
with regard to dynamic programming solution. 

The application of both the Pontryagin's Maximum Principle and 
Dynamic Programming to the hydro-thermal dispatch problem was considered 


by E.B. Dahlin in 1964 [5] and later, jointly with D.W.C. Shen [44]. 
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The general dispatch formulas obtained were applicable to H wide class 
of systems. These were the systems with hydro-plants having fixed 
head, varying head, and hydraulic coupling, both with and without river 
transport delay. 

The economic operation of a simplified model system and the long 
range operation of a multi-reservoir system is considered by I. Hano, 
Y. Tamura and S. Narita in 1966 [3]. They employed the Pontryagin's 
Maximum Principle to obtain the scheduling equations. 

The various factors involved in hydro-thermal coordination and 
their interrelation required for optimum generation are discussed by 
C.W. Watchorn in 1967 [3]. In a discussion of C.W. Watchorn's paper, 
G.S. Christensen points out that the scheduling equations obtained 
constitute only a necessary condition for optimality. He suggested the 
Steepest descent method to search for the global optimum mode of 


Operation. 


4.2 Power System with Fixed Head Hydro-Plants and Negligible Transmission 
Losses 

In this section, a hydro-thermal electric power system is considered. 
It is assumed that the transmission losses and head variations at the 
hydro plants are negligible. Imposing a constraint on the volume of 
water discharged over the optimization interval is equivalent to 
constraining the hydro-energy over the same interval. The formulation 
of the problem under this condition will fail to define the optimal 


hydraulic generations. This difficulty arises if one uses either of 
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the variational principles or the functional analytic technique adopted 

in this thesis. The main reason for this is the absence of second 

order terms of the hydro-generations in both the cost functional and 

the constraints. To avoid this difficulty in this particular problem, 

it is assumed here that the integral of the square of the hydro 
generation over the optimization interval is a prespecified constant [45]. 


This implies the presence of an upperbound on the hydro-energy available. 


4.2.1 Statement of the Problem 

A hydro-thermal electric power system is considered. The system 
has n generating plants, of which m are thermal. A prediction of the 
system's future power demand is assumed available. The problem is to 
determine the generation of each plant in order to minimize the operating 
cost of the system under the following conditions: 


1. The operating costs at the ith thermal plant is approximated 


by: 
2 
FyUPg (t)] a, + B;Ps (tf) + VPs, (t) $/hr (4.2.1) 


2. The total generation in the system matches the load. Transmission 
losses may be neglected. 
3. The integral of the square of the hydro generation over the 


optimization interval is a prespecified constant. 


4.2.2A Minimum Norm Formulation 


The object of the optimizing computation is 


NG 
i 2(t) }dt (40222) 
find n | ae + p,P. (t) ~ WP ee 2 
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Subject to the constraints: 


m 
Po(t) = 2 Po (t) +2 Ph (t) (4.2.3) 


i; 

f pest it Kk: i= m1,...4n (4.2.4) 
1 

6) 


The constraints given by (4.2.4) may be included in the cost 
functional defined in (4.2.2). The resulting augmented cost functional 


is then given by: 


I 
fil m 
Matte ahE ls a F, (t)) =f ie. + 2] 
fe) 
n 


510 rokeP, “(td iapae (4.2.5) 
i=m+] i 

Here terms explicitly independent of the power generations are neglected 
in the cost functional. The k.'s are unknowns to be determined such that 
(42.4) AS satisfied. 


Define the nx] column vector control function u(t) as: 


Tinea ines taney 


u(t) = col.[P. Gaaye US 
] 2 m m+] 


vey PL (t)] (4.2.6) 
n 
Define two nx1 column vectors L and M as 
= CO1.[By 5 «++> Bas Os -+es 0] (4.2.7) 
COLE lg u.ceate la (4.2.8) 


j= |r 
i l 


and an nxn diagonal matrix B by: 
Bye diagLy,> Yoo s+3 Ym? K atl? Lien kJ (4.2.9) 
Using these definitions in (4.2.5) and (4.2.3), the problem 
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reduces to minimizing 


t 
a 
Jo (u(t) ) =f ctu(t) + ul (t)B u(t) jdt (4.2010) 
0 
subject to satisfying 
Py(t) = M'u(t) (4.2.11) 


The control vector function u(t) is considered an element of the 
Hilbert space eto). This is the space of all n-dimensional vector 
valued square iitecrabie functions defined on the interval [O,J. The 
Space is endowed with the inner product definition 

<V(t) ,u(t)> ‘ V(t)B u(t)dt (4.2.12) 
0 
for every V(t) and u(t) in the space Lo pl0» Td The power demand function 
P(t) is considered to belong to the Hilbert space L[0.T¢] of the single 
valued square integrable functions defined on the interval [0.1]. The 
inner product definition for this space is: 
<x{t),y¥(t)> -f' x(t)y(t)dt LAS 2 aS) 
0 
for every x(t) and y(t) in L,[0,T-]. 
The cost functional given by (4.2.10) then reduces to 


Jo(u(t)) = <Veu(t)> + | u(t) ||? (4.2.14) 
where 


yl 


=! 


B= (4.2.15) 
The constraint given by (4.2.11) is seen to define a bounded linear 
transformation T: Lp plOsTeloL [0.7 ¢] of the form 

Pa(t) = Tlu(t)] (4.2.16) 


Following the same procedure as in Section (3.2.2), the problem 
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reduces to minimizing: 

J(u(t)) = |lu(t) + v/2]| CESORT 
subject to 

Py(t) = Tlu(t)] (4.2.18) 


This is a minimum norm problem of the form discussed in Chapter 2. 


4.2.3 The Optimal Solution 


According to the results cited in Chapter 2, the optimal solution 
to the problem formulated in Section 4.2.2. is given by: 
ue (ty = T"[P (t) # Tee) > (W/2) (4.2.19) 
The system's transformation T here is the same as that of the 
all-thermal problem solved in Chapter 3. Moreover, the spaces involved 
are essentially the same. Thus the derivation of the optimal solution 


as given in (4.2.19) follows essentially the same steps as given in 


Section (3.2.3). The optimal generations are thus found to be: 
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PMR OnS Sima s Re (4.2.21) 
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It can be seen that the sum of the optimal generations as given 


by (4.2.20) and (4.2.21) yields the power demand function Py(t)- This 
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is in agreement with the constraint given by (4.2.3). The (n-m) 
unknown k,'s can be obtained by substituting (4.2.21) in the corresponding 


constraint (4.2.4). 


4.3 Power System Containing Hydro Plants with Negligible Head Variations 


The problem posed and solved in this section is the same as the 
problem in Section 4.2, except that the transmission losses in the system 
are also included here. The transmission losses are assumed to be 
represented by the General Loss Formula. The General Loss Formula provides 
second order powers of the hydro-generations, which makes it possible to 
take directly the hydro-water draw-down constraint into consideration. 

The results of this section were first reported in [46]. It is 
easily shown that the optimal generations obtained in this section satisfy 
the variational conditions for optimality. 

4.3.1 Statement of the Problem 

Consider a power system with m thermal plants and (n-m) hydro plants. 
A prediction of the system's future load and water supply is assumed 
available for the optimization interval under consideration. The problem 
is to find the active power generation of each plant as a function of 
time over the optimization interval under the following conditions: 

1. The total operating cost of the thermal plants over the 
optimization interval is a minimum. 

2. The operating costs at the ith thermal plant are approximated 
by: 
kat: $/hr (4.3.1) 


3. The total active generation in the system matches the load 
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plus the transmission losses. 

4. The time integral of water discharge from each hydro plant 
is a prespecified constant amount. 

4.3.2 A Minimum Norm Formulation 


The object of the optimizing computation is 


Ts m 
Min ‘ties be 2 FC, (eat (4.3.2) 
1 iF 


The generation schedule sought must satisfy the active power balance 
equation: 
m 
Py(t) = ) P oD 4g Pee Ph be) - P(t) CANS 73) 
The transmission loss is a quadratic function of the active power 
generated by the system plants and is given by [10]: 


n n 
Do Pat IEG RAP) tells Blo ae Ro 4.3.4) 


with Bij 's and BG 's and Kg being the loss formula coefficients which 
are vee to be known, with the property: 


="Be- CPG] s\T, Seay Hi) (483%5) 


Furthermore, the water discharge at each hydro plant must satisfy the 
following constraint on the volume of water used during the optimization 


interval: 


q,(t)dt = b i ='mtT Ss 44m (4.3.6) 
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The ith hydro plant's active power generation Pi, Ge itis given by: 
a 
n 444 (t)h, (t) 
Ph, ft) = T> Than 7 KW (4.327) 
Assuming constant head and efficiency at each hydro plant, the constraint 


given by (4.3.6) is equivalent to the following energy constraint: 


Vs 
i P, (t)dt = K, ie einer ena (4.3.8) 
; a 
(@) 
with 
a ern 
es ‘le ent 
Saas 5 (4.3.9) 


The constraint given by equation (4.3.3) can be added to the 
integrand of the cost functional given by (4.3.2) using theunknown function 


6(t), so that a modified cost functional is obtained: 


: 
feogM 
dls. (ts Pp, (t)) = [By (P (t)) + oft) (Pp Ct) 
} yp (t) 
- PEs At) c= Deu yrs 
a Si fet nis 
P (t)) jdt (453.10) 


Substituting for P, (t) from (4.3.4) and Fekete) from (4.3.1 )£4(4.3710) 
j 


can be written as: 


Ts m 
Jy(P, Ct)» Py (ED) = { [2 (ey * CE)EB;g - WD)Pe Ct) 
(@) 
n 
+r) BBN 17 Ph (t) 
TyePe(e(t)e oft) Ef 
, tat 6 
: ibytts, : i=) j=1 
P. (t)B:P.(t) Jae (4.3511) 
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Here terms explicitly independent of the power generations are dropped. 
Define the nx1 column vector : 


u(t) = POUR Mr ty sees Py (t),...5P) (t)] 


m+] n 
(qoe 12) 
L(t) = col.[(8, + C,0(t))s..0s m+po(t)e+++2C,6(t)] 
(2 selon 
where 
C. = Be G - | 1) “waleshow all (4.3.14) 
And the symmetric (nxn) matrix 
B(t) = (b;;(t)) (4.3.15) 
with 
b; ;(t) a B; 9(t) Ivy So acharet «tefl (i 7 j) 
b. ;(t) 7 Yj + B.,6(t) eid > Mm 
b..(t) = B,,e(t) i = m+l,...,n 
then (4.3.11) becomes: 
fr fl 
J,(u(t)) = f {kL (t)u(t) + uo (t)B(t)u(t) dt (4.3.16) 
0 
If we let 
v(t) = B'L(t) (4.3.17) 
equation (4.3.16) can be written as: 
T 
f 
J, (u(t)) aif (vT(t)B(t)u(t) + ul (t)B(t)u(t)3dt (4.3.18) 
0 
T 
f 
or 3, (u(t) a ce(u(t) + MED) Tact) (u(t) + YH); 
0 
I 
aie {t) p(t) WY) qat (4.3.19) 
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The last term in the integrand of equation (4.3.19) does not 


depend explicitly on u(t), so that it is necessary only to consider: 
/f = V 
Dy = fie =—yTa(t)(u(t) +2 .ae (4.3.20) 


The problem is now transformed to that of finding a control vector 
u(t), that minimizes the cost functional given by (4.3.20) while satisfy- 
ing the energy constraint given by (4.3.8). Notice that e(t) will be 
determined so that the optimal control vector e hieries the active 
power balance equation (4.3.3). 

The control vector u(t) is considered to be an element of the 
Hilbert space Lp pl0sTelof the n-vector valued square integrable 
functions weaned: on Lost ubose inner product is given by 

<W(t) ,u(t)> five t)B(t)u(t)dt (4.3.21) 


0 
for every V(t) and u(t) in Ls"4L0.T¢. provided that B(t) is positive 


definite. This means that u(t Jel, oe Te] which minimizes: 
J(u(t)) = | [u(t) y+ Xt e) || (4.3.22) 
and satisfies (4.3.8) is sought. 
Define the (n-m)x1 column vector. 
Be cole (Kase ay J (4-323) 


and the nx(n-m) matrix 
M = col.[0,1] (4.3.24) 
with 0 being the mx(n-m) matrix whose elements are all zeros, and I is 


the (n-m)x(n-m) identity matrix, so that the constraints of equation 


(4.3.8) can be expressed as: 


is f wutnee (4.3.25) 
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Equation (4.3.25) can be shown to define a bounded aaa. trans- 
formation T: LD plOsTp PR). 
The real vector space a(n-m) is endowed with the Euclidean inner 
product definition: 
“Kate sake Y (4.3.26) 
for every X and Y in paz) 


Equation (4.3.25) is now written as: 


Eee rLutt)) (4.3.27) 
with T[u(t)] given by: 
ie 
le | wluceyae (4.3.28) 
(9) 


The problem under consideration has been reduced to finding the 
optimal generation vector u(t) which minimizes the cost functional given 


by (4.3.22) such that (4.3.27) is satisfied for the given vector &. 


4.3.3. The Optimal Solution. 


In view of the results cited in Chapter 2, there is exactly one 
optimal solution to the problem formulated in the previous section, namely 


the optimal vector. 


u(t) = Te + Ted - Yel (4.3.29) 


where T" is obtained as follows: 


ie the adjoint of T, is obtained by using the identity: 


< > = <Te,u(t)> 
g-Tu(t)> (nom) T é,u(t) t0,T, (4.3.30) 


(n-m) 


where in R we have: 


<e,Tu(t)> = ¢!Tu(t) (4.3.31) 
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using (4.3.28), the following is obtained: 


or 


where use is made of 
Br '(t) = (B '(t))! 
But in ue p LO; A) we have: 


f 
<B!(t)Me,u(t)> = | (B7! (t)Me)"B(t)u(t)dt 


— oe 


hence, by equation (4.3.30) it is seen that 
Te = Bl (t)Me 
which defines the adjoint of T. 
Next we find the transformation J given by: 
Je = TT 6) 
using equations (4.3.28) and (4.3.36) as 
Te Ta-1 
IS fee (t)MdtJé 
0 


which yields ie 


(4.3.35) 


(Ater36)) 


(4.3.38) 


(ar 3.30) 
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Finally, using the definition of the pseudo-inverse transformation 


Th one obtains: 


Thus the optimal generation vector u,(t) is given by: 


(4.3.40) 
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(4°3-41)} 

The optimal solution obtained here may be written in component 
form. However, this cannot be done unless the number of plants is 
specified and the inverse of B(t) is evaluated. The optimal solution 
contains the unknown function 6(t) which is determined such that the power 
balance equation is satisfied. 

4.3.4 Implementing the Optimal Solution. 

Without loss of generality, consider a power system whose loss 
formula contains zero off-diagonal coefficients, then B(t) can be written 


in partitioned form as: 


B(t) = diag[B..(t), By(t)] (4.3.42) 
with 

Bog (t) = dial (yy + 0(t)B)y)o++-9(yq + 9(t)Bpm)I (4-3-43) 
Let 

L(t) = col.[L,.(t)s Li (t) J (4.3.44) 
with 

Leg (t) = col.[(8, + 0(t)C)>-+-2(8, + 2(t)C,)] (4.3.45) 

Li (t) = col.[6(t)C 45... ,0(t)C, J (4.3.46) 


Substituting in equation (4.3.41), the optimal solution for this 


particular system is obtained as 


where 
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so that the optimal power generations are obtained as: 


[8, + C, o(t)] 
P. (t) = - ——————_ W = Dane git (4.3.49) 
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The last expression contains the unknown function 6(t) which will 
be determined so that the power balance equation (4.3.3) is satisfied. 


Using equation (4.3.4), equation (4.3.3) is written as: 
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then substituting in (4.3.51) for Be (t) and Ph (t) as given by 
E E 
(4.3.49) and (4.3.50) respectively we get 
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where 

x(t) = 07! (t) 

n 
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D =) ise Ts 

fetthl » thee! 

Ce aE 

Li K. + 3 i T° Me Loree ol 


: n ane 
* “Lo i=m+] 4B 
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By) ea ein uae 
pit) = 4B, .° i= 1,...,m 


Thus equations (4.3.49), (4.3.50) and (4.3.52) completely define 
the optimal generation vector required. 

4.3.5 A Computational Example 

The method which is used to find the optimal solution is a scanning 
process that searches for a function x(t) which is a solution to (4.3.52) 
and simultaneously yields a positive definite matrix B(t) (which makes the 


definition of the space 8") 0,7] valid). Furthermore, the optimal 
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generations must be physically realizable in the sense that no negative 
or complex active power generations are permissible. 

The procedure followed in solving (4.3.52) is one of discretization 
which will result in a system of N simultaneous algebraic equations in 
N unknowns (Xy 90+ sXy) which can be solved by any classical iteration 
technique. 


A study was made of a sample system with the following particulars: 


i 
—d 
s 


Number of thermal plants (m) 


Number of hydro-plants (n-m) = 2. 


Loss formula coefficients 


4 a a 


Baa = -1..6000°x 10° =°292000, x10, SBas ae) 56000e. 10" 


515) 
= 0.00 Boo = 0.00 B39 =90.00 


1] 
Big 

Thermal Plant's cost equation: 

S 


F(P.(t)) = a + 4.0 P.(t) + 1.20 x 107° p@(t)  $/Hr. 


The energy constraints on the hydro-plants: 


K 3600.00 MW-Hr. 


2 
2400.00 MW-Hr. 


aw 
iH} 


Example. (1): 
Let 
Py(t) = 400.00 MW Dict culem tend aamcreliinr 
then equation (4.3.52) reduces to a second order algebraic equation in 


x(t), which yields a feasible x(t) given by: 


x(t) = 0.21629 0.< beset, 
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so that we get the following generation allocation by using equations 


(4.3.49) and (4.3.50): 


Optimal thermal power generation P(t) = 160.68 MW(O0 < t < Te) 
Optimal hydro-power generation Pi, (t)” = 150.00 MW(0 < t < Te) 
: cs 
Ph (t) = 100.00 MW(0 < t < Te) 
: ees 


The result of this example agrees with the obvious assumption of 
constant, hydro-generations at the average values satisfying the energy 
constraints and the value of thermal power generation that satisfies the 


power balance equation. 


Example (2): 


Let 
of 
Py (t) = 400.00 MW(O < t < 5) 
T 
= 600.00 MW(5— < t < Te) 
T. = 24 Hrs 


then equation (4.3.52) reduces to two simultaneous algebraic equations in 


two unknowns Xy and Xo where 


x(t) = x) 


2 a 
The values of xy and Xo satisfying these two equations and which 


simultaneously yield a feasible optimal power generation are found to be: 


x > 0.20008 


Xo 0.19349 
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with the corresponding optimal power generations given by: 


P(t) = 249.53 MW, PL (t) = 114.42 MW » Pi (t) = 49.40 MW 
2 3 
T 
Ons tus x 
P(t) = 288.21 MW,P. (t) = 185.58 MW , Py (t) = 150.70 MW 
2 3 
Al 
sex tc Ty 


4.4 Power System with Variable Head Hydro-Plants and Transmission Losses 


The problem considered in this section is more general than the 
ones dealt with in the previous sections. The results obtained in this 
section are shown to agree with Kron's results [11], that is, for the 
case Of one thermal and one hydro plant system. The work of this section 
was reported by the author jointly with G.S. Christensen in [47-48]. 

4.4.1 Statement of the Problem. 

Consider an electric power system with m thermal plants and (n-m) 
hydro-plants. A prediction of the system's future power demand and water 
supply is assumed available for the optimization interval. The problem is 
to find the active power generation of each plant as a function of time 
under the following conditions: 

1. The total operating costs of the thermal plants over the 
optimization interval is a minimum. 

2. The operating costs at the ith thermal plant are approximated by: 

FOP, (t)] = 0 #8; Pet) + ¥q P(t) $/Hr. 
(4.4.1) 


3. The total active generation in the system matches the load plus 
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the losses. 

4. The transmission losses in the system may be represented by the 
general loss formula. 

9. The time integral of water discharge for each hydro-plant is a 
prespecified constant amount. 

6. The ith hydro-plant's reservoir is assumed to have vertical sides 
and small capacity. Head variations are related to the discharge by the 


reservoir's dynamic equation. 


Sha (t) = i,(t) - q.(t) (4.4.2) 
7. The tail-race elevation at any of the hydro-plants does not 
change with the water discharge. 
4.4.2 A Minimum Norm Formulation 
The object of the optimizing computation is: 


FopPebt dt (4.4.3) 
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The generation schedule sought must satisfy the active power balance equation: 


m 
Py (t) = ” P ae Ae one - P(t) (4.4.4) 


The transmission power loss is given by: 


(4.4.5) 


Furthermore, the water discharge at each hydro plant is to satisfy 


the following constraint on the volume of water used over the optimization 


interval: 
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The ith hydro plant's active power generation Ph (t) is given by 
j 


q,(t)h, (t) 
P io ee c 7 = m+] wn (4.4.7) 
j j 
with 
G. = 11.8/n, 


The effective hydraulic head at the ith hydro plant h.(t) may be expressed 
using (4.4.2) a 


t 
hy bhler etesy cay a)do - fotos + h.(0) (4.4.8) 
8) 


5 om til ee are 
Thus substituting (4.4.8) in (4.4.7), h.(t) is eliminated. The expression 


for the hydro-power is then given by: 
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9 (4.4.10) 
- -] 7 = +] 
Eee (S.G.) 1 = mrls.c yn 
(AAT) 
then (4.4.9) reduces to , 
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(4.4.12) 
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It is convenient to introduce the volume of water discharges variable 


Q, (t) by the following definition: 


t 
Q.(t) = [9;(o)de 1 =imtl, .. Qn (4.4.13) 
0 
thus 
q,(t) = Q.(t) i Spine ae oi (4.4.14) 


This reduces the power generation as given by (4.4.12) to 


1 ¢ Q.(t)[N. (t) at E. Q.(t)] (4.4.15) 
1 = tn 
And the constraint on the volume of water discharged over the 


optimization interval given by (4.4.6) reduces to: 
Q, (Te) = b, i = mtl,...,n (4.4.16) 


The problem now is to minimize the cost functional given by (4.4.3) 
subject to satisfying (4.4.4), (4.4.5),(4.4.15) and (4.4.16). Including 
(4.4.4), (4.4.5) and (4.4.15) in the cost functional yields: 

F,[P, (t)] + a(t)[Pp(t) 
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Using the following relations: 


with 


ri(t) = A, (t)N,(t) + n.(t)N.(t) (4.4.20) 


We can write the cost functional (4.4.17) as: 
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Dropping terms in (4.4.21) that are explicitly independent of the 
control functions re Goce Pi (t) and Q. (t) over the interval [0,7] it 
is only necessary to consider minimizing: 
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Celina 
Define the (2n-m)x1 column vectors. 
u(t) = col.[P.(t).p,(t) Q(t)] (4.4.23) 
Mee) = (cot. (Lo (e) Lp (t)) u(t) (4.4.24) 
where 
P.(t) = sees yc (4.4.25) 
P, (t) = SO ere macnn (4.4.26) 
L(t) = col.[{8, -(1-B.  )A(t)},.-. 518, catia ive) acd 
(A427) 
L(t) = col.[t-[n q(t) + (-Brayzyo)A(t) Is... 
f=[n.(t) Fie B att) Is! (4.4.28) 
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L(t) =: col.[-r.,(t)s...s-r(t)] 


where 


B(t) = (bj; (4) 


>) 


is the (mxm) matrix whose elements are: 


b. (t) = ss A(t) B. ie 
ig?) a r(t)B, an T5J 
Beret ba (Ct) ): 
—=sh _» Man 
is the mx(n-m) matrix whose elements are: 
= \ awe 
De r(t)B. | 
nae 
Bei(t)rse(bes -(t)) 
—ns ‘Ins 
is the (n-in)xm matrix whose elements are: 
cv = oe 7 = 
ie” r(t)B, | 
J — 
and By (t) is the (n-m)x(n-m) diagonal matrix: 
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Using these definitions (4.4.22) becomes: 


Let 


and we have 


S t) 
0 0 
where vi 
-| 
C(t) = [B.(t) - Bo (t)Bp (t)B, (t)] 
-] 
C(t) = [B,(t) - B (t)B.'(t)B, (t)] 
-] 
Cont) = [-Bl'(t)Boy (t)C, (t)] 


provided that the inverses in the last equalities exist. 


(4.4.39) can then be expressed as: 


vl (e) = Ev (e) fu) Ct) [yg (t)] 
where 

Motte Li (oyce(t) + L(t yc, (t) 

it 1 ih 

vit) = i t)eta(t) +L, (t)c, (¢) 

ip ener 

Vale) = Lolt)ey(t) 
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Furthermore, (4.4.38) can be written as: 


- E42 a(t) Wt) at (4.4.50) 


The last term in the integrand of (4.4.50) does not depend explicitly 


on u(t), so that one needs to consider only minimizing: 


Tr, 
V(t 
Jo(u(t)) Ee! {fu(t) + =t)57 B(t (t) + = ) (4.4.51) 
0 
subject to the constraints. 


Q; (Te) = b. 7= mt] ja 45n (4.4.52) 
Define the (n-m)x1 column vector: 
b = col.[b .15++-»b,] (42453 } 
and the (2n-m)x(n-m) matrix M by: 
¥ d 
= Bk yas. (4.4.54) 


0 being the (n-m)xn matrix whose elements are all zero, and I being the 
(n-m)x(n-m) unity matrix. 
Thus the constraints of equation (4.4.53) reduce to: 
eae 
b “f Mu(t)dt (4.4.55) 
0 
The control vector u(t) is considered an element of the Hilbert 
Space a “ro, nf cl of the (2n-m) vector valued square integrable functions 


defined on [0,7] whose inner product is given by 
MG 
t 
<u(t) u(t)> = fv" (e)B(t)u(t)at (4.4.56) 
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[0, Te], provided that B(t) is positive 


The given vector b is considered an element of the real space Rin-m 


with the Euclidean inner product definition: 


“Nel = xly (4.4.57) 


for every X and Y in ay 


Equation (4.4.55) then defines a bounded linear transformation 


ie Boe ote, This can be written as: 


= T[u(t)] (4.4.58) 
and the cost functional given in (4.4.51 reduces to 
Jo(u(t)) = ||u(t) + +i \\2 (4.4.59) 
Finally, it is only necessary to minimize 
) = |Ju(t) + aH) | (4.4.60) 
subject to 
= TLu(t)] 
(n-m) | 


for a given b in R 


4.4.3 The Optimal Solution 


Tne optimal solution to the problem formulated in the previous sub- 


section using the results of Chapter 2 is given by: 


u,(t) = T'[b + T(v/2)] - (W/2) (4.4.61) 


where rT’ is obtained as follows: 


me the adjoint of T, is obtained using the identity: 


<—,1U> = a a tes 
pein g8es rg, tl (4.4.62) 
Let 
. 4.63 
3 COl.[E spose + sbqJ (4.4.63) 
* 
2 4.4.64 
Te = col.{T,. 1,» Ty) ( ) 
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with 


Ty iS col.[t,, meee td (4.4.65) 

Th Fecot Et > --.5 th] (4.4.66) 

Ty r COl.Lu yiseees WA (4.4.67) 
(n-m) 


then in R , the inner product of the left hand side of (4.4.62) is: 


<E,Tu> = a f Mu(t)dt (4.4.69) 


where Q; (0) = 0 according to (4.4.13). 


And in LS 


2n-m) 
2.0 


[0.1], the inner product of the right hand side of (4.4.62) 
is 


. 
on oe “foe u(t)dt (4.4.70) 
O 


Using (4.4.64), (4.4.30) and (4.4.23) this reduces to: 


+ Ty| ByQidt (4.4.71) 


Thus the identity of (4.4.62) using (4.4.69) and (4.4.71) reduces to: 


+ T,'B.Q(t)}dt (4.4.72) 
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Equation (4.4.72) is satisfied for the choice: 


1) el t e[0,T¢] (4.4.73) 
T(t) = 0 t e[0,T¢] (4.4.74) 
Hes | t e[0,T,) (4.4.75) 
iat) : ce ey (4.4.76) 


Thus Te is given by: 


* aan 
T Le] = = & (4.4.77) 


where 0 is the nx(n-m) matrix whose elements are all zeros, and u(t) is the 


(n-m)x(n-m) diagonal matrix given by: 


u(t) = diaglo,,(t),....6,(t)] (4.4.78) 
e.(t) = 0 t e[0.T,) 
2 
> Ey f(T.) Boe (4.4.79) 


The operator J is next evaluated as: 


J(z) = TET] (4.4.80) 


Using (4.4.54), (4.4.55) and (4.4.77) this reduces to: 


Te : 0 
p i 
or Te 
- d 
J(é) = qe u(t) —& dt (4.4.82) 
O 


so that J(é) = [ulTe) - ulo)]e (4.4.83) 
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but by (4.4.79) u(o) = 0, then 


J(E) = ulTe)é (4.4.84) 
This yields: 


Ie) = w Tee (4.4.85) 


Finally, the pseudo-inverse operator vt is obtained from the 


definition 
Tis) = Tbe (4.4.86) 
using (4.4.77) 4.4.85) this yields: 
eo es x (4.4.87) 
or 
re |as (4.4.88) 


From (4.4.54) and (4.4. ie one obtains 


T(V(t)/2) = ZLVe(T) - Vq(o)] (4.4.89) 
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The optimal solution is now found by substituting (4.4.89), (4.4.88) and 


(4.4.46) in (4.4.62); 


Ree) = - FV, (t) (4.4.90) 
é 
1 , 
af a> aC) (4.4.91) 
Wt) 
Q, (t) = u( Tare] + 5 (ol 1) : Vo Dynes i 
(4.4.92) 


which can be reduced by using (4.4.79): 
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V 
since Q, (0) = - ae GC. 


Thus according to (4.4.57), (4.4.48) and (4.4.49) the optimal solution 


aS: 
edt} gon ACL (t)L, (t) + Gece) (4.4.95) 
AC pC (tLe (t) + o)(t)Ly (t)] (4.4.96) 
1 ae lea 
Q-(t) = - pig (t)Ly(t)] (4.4.97) 


It is noted that the optimal solution involves the unknown functions A(t) 
and n,(t) which are to be determined such that the constraints (4.4.4) and 
(4.4.15) are satisfied. 

4.4.4. Comparison with Kron's [11] Equations 

Consider an electric power system with one thermal and one hydro plant. 
The optimal solution for this system is given by (4.4.95), (4.4.96) and 
(4.4.97) where all vector quantities reduce to one element vectors, since 
m = 1 and n-m = 1 in this case. 


Here we have using (4.4.27), (4.4.28) and (4.4.29) 


L(t) Sate G A(t) - Biol} (4.4.98) 
L(t) = - {n,(t) + a(t)[] - Bog 1} (4.4.99) 
Lg(t) = - ro(t) (4.4.100) 


And using (4.4.41), (4.4.42), (4.4.43), (4.4.44) and (4.4.45) 
B 
22 
: 4.4.101 
C(t) Ty F XUEYB) | Bpp ~ 2CE)B] 289] ( ) 
(yy ef: d(t)B, 4) 


CA) = 
n' ) Mt) La YO t Bay B55 - r»A(t B, 9B. (4.4.102) 
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Co, (t) i ge 1] + MEE, C(t) (4.4.103) 
C s(t) =.- Ma. (t) (4.4.104) 
hs Boo S ie 

= Z 
C(t) = Eon (t) (4.4.105) 


Thus the optimal solution is: 
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[{8, - r(t)[1 - Byy]} + Baie) 


+e b= Bool} ] (4.4.106) 
q Ly, + A(t)B, 4] 
Fie FE 2h (tiny ty op eee, 
By oA(t) 

- 7 + At)By, {By - r(t)[1 - Bi od} 

- {n,(t) + r(t)[1 - 1} (4.4.107) 
SEN A ay, (4.4.108) 
E EBf, (t) 


Furthermore, eliminating n,(t) between (4.4.106) and (4.4.107) yields 
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where 

A= (yy + A(t)B, 4) Boo - (t)By Bo, (4.4.110) 

Equation (4.4.109) can be rewritten as: 
aa Re bt) + [ey ~ alt) SALEESeae EIR Tali 
2atig 22) 1] 
2A A(t)Bo, 
een Ee i 
22\*) W E 

or 


2, (Oly + A(t)B),] + Ca, - a(t)(1 - By)] 


+ a ee = 0 CASA TT) 
This equation should be satisfied if P, (t) and P, (t) are to be optimal. 
E 
The optimal volume of water discharge is reduced to: 
1 : no(t) : 
Q, (t) oe E,tNa(t) ay No(t) pacer. (4.4. 12) 


This results from applying (4.4.20). 


Furthermore (4.4.106) and (4.4.107) can be written as: 
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[no(t) + a(t)(1 - Bog)] (4.4.114) 


Multiplying both sides of (4.4.113) by Bio» (4.4.114) by B45 and adding, 
the following is obtained. 
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Chee TS) 
If one lets 
D, (t)S_(t) = (1 - Boo) ~ Piss - Booth (t) 
(4.4.116) 
then no(t) = - d(t)D, (t) (4.4.17) 
no(t) = -[i(t)D,(t) + a(t)D, (t)] (4.4.118) 
so that the optimal volume of water discharged given by (4.4.112) becom 
] : d(t)D, (t) 
Q, (t) = EptN2{t) * Nott) x7Ey DCE} + AEG, (E) (4.4.19) 


Kron's scheduling equations [11] are: 
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SOPL ee eevee m at =p (4.4.121) 
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This is exactly (4.4.111), thus the optimal solution obtained in sub- 
section (4.4.3) is Kron's first equation. 
Furthermore, one finds that 


oP 

Lb 
eS PB 

oP 


And by (4.4.116) this reduces to: 


aPy 
D(t) =] - SP (4.4.125) 


h 
Thus Kron's second equation (4.4.121) transforms to: 


== 2Bier (ti) = 2B 


22° h P(t) 


20 L24S 


aot + Sha(t)D(t) —"] = 0 (4.4.126) 


By (4.4.7), (4.4.14) one obtains 
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he wie 
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And by (4.4.15) 
i 4.12 


Using (4.4.127), (4.4.128) and (4.4.11), equation (4.4.125) becomes: 


EQA(t)O(t)Qo(t) + FeLAlt)D(t)INp(t) - E,Qo(t)I]= 0 


or 


No (t)a(t)D{t) . 
Qo(t) = EANo(t) + greypEy aa ye) (4.4.130) 


But this is precisely (4.4.119), so that the optimal solution obtained 


is the second equation given by Kron. 
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4.4.5 Implementing the Optimal Solution 


The optimal solution obtained in subsection (4.4.3) contains the 
unknown functions n,(t) and A(t). These can be determined by substituting 
the optimal solution in the corresponding constraint equations. The 
resulting equations in n (t) and A(t) are generally nonlinear. In order 
to get deeper insight into how the actual solutionis obtained, the 
following simplifying assumptions are made: 


1. The system is characterized by a loss formula where: 


a. Bi, = 0 i7#j 902 lhc egn 
b. EG = 0 lr eeilecrccct 
Cc Ko = 0 


2. The rate of natural water inflow to the reservoirs is constant. 
These assumptions will represent no loss of generality. 


The first assumption implies the following in (4.4.30). 


Bo (t) aal®) (4.4.131) 
By (t) =.0 (4.4.132) 
B.(t) = diagl(y, + A(t)Byy)s-+ 05m + A(t)Bog)] (4.4.33) 
By(t) = diagla(t)B nay) (meq)ee- + 2t)Ban] (4.4. 134) 


And in (4.4.27) and (4.4.28): 
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Thus the optimal solution given by (4.4.95), (4.4.96) and (4.4.97) 


reduces to: 
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The hydro-power constraint given by (4.4.15) is 
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Differentiating both sides of (4.4.155) yields: 


Q; (t) = FAN, (t) +N, (t)x, (t) +N, (t)X,(t)] 
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E i 
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Recalling (4.4.10) for N.(t) and using (4.4.11) one obtains: 
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Let 


a.(t) = 2B, .E. Tea) 1 = mtl,...,n (4.4.163) 


Then (4.4.162) reduces to: 


nj(o) -1 
1+ see exe Ox (¢) dt] 
0 


i. (t) 
be) Se ee eee t 
‘ a, (t)x;(t) Tyrer “40 
Tivaunele. can (4.4.164) 


Under the assumption of constant water inflow i, (t) = 0, thus (4.4.164) 


reduces to: 
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The optimal Q;(t) given by (4.4.146) reduces to 
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At t = Te, equation (4.4.165). is: 
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Thus solving (4.4.149) and (4.4.165) subject to the boundary conditions 
(4.4.172) and (4.4.173) completely defines the optimal schedule. It is noted 
that (4.4.149) is an algebraic equation in X(t) and the (n-m) unknown 
functions n.(t). On the other hand (4.4.165) is a set of (n-m) first order 
nonlinear differential equations in x, (t), where x. (t) and n.(t) are related 
by (4.4.153). 


In the computerized search for the above mentioned unknown functions, 
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by utilizing the physical significance of each variable involved. In 


addition to this, restrictions on the variables can be obtained so that 


the functional analytic formulation adopted is a valid one. 


Consider the optimal thermal power generation expression of (4.4.144). 


If the thermal power generated is to satisfy the following practical 


limitation: 


cuba (tC) exe i Gt osm 


Max. 
Thus using (4.4.144) one obtains 


By Pedy P 


Min. i i Max. 
eee SCF= *“A(tB) ecco Ge 
Les eBiihMin, weg ii Max, 
Let 
B : “viPMin, 
y = Max [ ] 
SS ae ee inal A AT 
By + CyyPuay ; 
‘ = Min [ 
Max =] sm V- 2B iP Max 
Then 
Amin < A(t) Max 


46 OU 


(4.4.174) 


4 4 at75) 


~-— 


(4.4.176) 


Use 


4.4.177) 


Furthermore an expression for the effective head h(t) can be obtained 


utilizing (4.4.8), (4.4.10) and (4.4.146) as: 
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x. (t) = - j = m+tl,...,n (4.4.181) 
Tat.) 


Here V(t) is the volume of water stored in the reservoir. In the case 


when this volume is restricted between upper and lower bounds given by 


< V.(t) < V m+1y...9n (4.4.182) 
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Furthermore, if the hydro-power is restricted such that: 
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This combined with (4.4.185) yields 


n, (t) > 0 T= mtl,..csh (424,189) 


The restricted search area in the (n-) phase plane is shown in 
Fig. (4-1). It is noted here that (4.4.177) and (4.4.189) guarantee that 
the matrix B(t) in the inner product definition is positive definite. 

It is worth mentioning here that it is assumed in the problem formulation 
that the inequality constraints are not violated. Thus these constraints are 
not included in the cost functional to obtain the optimal solution. The 
optimal solution is then implemented in a way that confines the search area 
to the regions of the space of the unknown functions where the inequality 
constraints are not violated. This agrees in principle with the nonlinear 
programming approach to this type of problems. In the Kuhn-Tucker method, 
unknown multipliers (called the Kuhn-Tucker Multipliers) are associated with 
the inequality constraints for inclusion in the cost functional. These 
multipliers are set to zero as long as the inequality constraints are not 
violated [49]. Thus an optimum solution is obtained by scanning the whole 
space of the unknown functions. If the solution obtained violates any in- 
equality constraint, the corresponding variable is set to the nearest value 
that does not violate the constraint. The main difference between that 
method and the method adopted here is that the search region in the latter 
is smaller than that of the former. This reduces the computing time 
considerably and leads to good estimated values for the unknown functions. 


Consider Equation (4.4.165) which can be rewritten as: 


x. (+) = f.[x.(t), Vaan n.(o)] 1 =m lee. s.n (4542.190) 


where 


8 aah ses "aunt banstety Son. oy rer 
| ee anghtw foe’ eerree ans nteddo, oF rao 3800 ve at 
8916 flovsa2 and, bent tod Ses VBW 6) Of bednsnpgi on ate 


eu appa ans aver anotd ant wordy atts TO, 93k 
_ sbi non a3 a Fu algisning ny 2957106 eka ae oi vd 
\sbongem plouT- itu sat aT ama l on to. ant aie of tt ty 
“fd beteisoze5 S48 (a9 atiqra tum Hout neu aif botfso). 2h 


- gzadT Thott anu $200 Sit nt “nokevfont. 1, angtercanes 


“ton 318 aintsrierios wt feupant ailt 28) enol eh MBS ( rs 32 a 


_ stot ai ‘eoineo Nd bentetdo. at noisuton ne. wT .f 


- ygme: ssid zsobbay. tT emt ait % nly we 
| ‘onan oft nor 2aut BV. | beseni tes boop aba bas 
oe nsssiey © od: ASS ota oe = ay Hor oe 


he teu | ’ hy be _" ; me — e i, 
a » recta FU yn eho 


Ff; [x;(t),a(t) on, 


iE, 
(o)] = -{[] + {Ln (oexetf x," (tat) 
0 


/r(t)}1/a, (t)x, (t)1-1 


(e-em (4.4.191) 
Let 
af;(s) = t50x;'1)(s),a(s)ong(o)] - #,Lx, @)(s),a(s) sn, (0) ] 
7 = mtl,.:.,n (4.4.192) 
Then using Equation (4.4.191) this reduces to 
1 (1) (2) 
ee = ye pra rae ee oe )) 
+e aa(s)xy' (six, oes) gee ee 
(x, §2) (5) Ndo 
+ Risa, Ngo 
ACs) Rt 
Seeley aes: 
(x (o)) ‘do 
- (le el ja ot reer 
(424° 193) 
The following inequality is true for all s: 
1 (1) (2) 
re < Exc pn (S) = Ke cts I 
| A .(s)| fa MP ge isixy Uo 1 1 


69 


Eee es 


saint i 7 raat! ha ivy 


70 


S 
Here the fact that expt{ Ex; (2) (6) "ldo < 1 was used. 
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Furthermore (4.4.198) provides exactly a Lipschitz condition given 


by 
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This means that Picard's iteration process given by: 


t 
eter ine, pts x, (0) + [lx (™(s).a(s).n 


j 
fe) 


is guaranteed to converge to a solution of (4.4.190) for given n. (0) and 
US). 

4.4.6 Practical Application 

A computer program was written to solve (4.4.165) and (4.4.148) for 
x.(t) (i = m+1,...,n) and a(t) to obtain the optimum generation schedule 
Figure (4-2 ) shows the flow chart for this program. An initial estimate 


of the function A(t) is made. This was taken as 


oe P AMin 


where Aan is given by (4.4.177). For each hydro plant, the initial value 


t eL0,T,] 


of n.(t) was estimated as: 
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‘Min 
A solution to (4.4.165) was then obtained for each hydro plant. This 


was done by utilizing the Picard's algorithm given by (4.4.201) and (4.4.202). 
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The value of x. (Te) obtained here was compared with the boundary condition 
given in (4.4.172). If the error in this step is large, the estimated 
value n,6)(0) was modified to the correct direction that minimizes the 
error. When all the x, (t)'s were obtained the corresponding n.(t)'s were 
evaluated using (4.4.154). 

Thus equation (4.4.148) becomes a (2m+2)nd order algebraic equation 
in A(t) for every t. This is solved for A(t) in the region given by 
(4.4.171). The A(t) obtained was then taken as the initial estimate 


instead of du and the process is repeated. If the difference between 


in 
two successive evaluations of X(t) was less than a prespecified amount the 
iteration process was stopped. The last step is to evaluate the optimum 
generation schedules as given by (4.4.144), (4.4.145) and (4.4.146). 

This program was applied to a hydro-thermal system with four thermal 
and three hydro-plants whose particulars are summarized in Tables 4.1 and 4.2. 


The optimal generation schedules and the assumed power demand curves are 


shown in Fig. (4-3). 
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Table 4.1 


Thermal Plants' Particulars 


Table 4.2 


Hydro Plants' Particulars 
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Figure 4.1 The Search Area in the f-n plane. 


L(t ) 


Ue oe hia) 2a sn ae ‘d 


ee Oe | De ee 
N i ye 2 tee cf 


a it n/t aoa 
a . ae ; - " | . pm ¥ i! WW a 5 1 sit ie. , i? ia a : vou ry 


rey 


a fe: 
yaa ee ve 
| 


iy 


Initiallize 


(9) (4) 


Select 
nj (o) 


Solve (4.4.165) f 
x.(t) using C.M.P. 
4.4, 


get n,(t) by (4.4.154) 


Solve (4.4.148) for 
] 
r 


Figure 4.2 Computer Program Flow Chart. 
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Figure 4.3a The System's Power Demand and Optimal Generations. 
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Figure 4.3c Optimal Thermal Output and Transmission Losses. 
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CHAPTER V 
POWER SYSTEMS WITH COMMON-FLOW HYDROPLANTS 


9.1 Background 


The problems considered in this chapter are characterized by the 
presence of common-flow hydro-plants in the system. These problems are 
of a complex nature. The hydraulic coupling between plants on the same 
stream is a factor that is not present in the problems dealt with in 
Chapter 4. The importance of including the time delay of flow between 
coupled plants in the optimization scheme was pointed out by C.W. Watchorn 
and R.A. Arismunander in separate discussions of [7]. 

Among the early contributions to this problem is H.A. Burr's work 
[50]. He developed loading schedules for a two-plant common-flow hydro 
system but the assumptions made were too simplifying. Later, P.R. Menon 
[51] used the Euler equations for constructing sets of minimizing sequences 
for a three-plant hydro-thermal system. It is noted here that the system 
considered by Menon wasof low-dimension, that is, a small number of hydro- 
plants were considered. Also this was a long-range scheduling problem. 

In [44], E.B. Dahlin and D.W.C. Shen treat the problem of a power 
system with hydro plants on the same stream using the Pontryagin's Maximum 
Principle. They used a river flow model which introduced a large number of 
differential equations and boundary conditions. This was a definite 


contribution to the theory of economy scheduling. Unfortunately this model 
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made the problem more difficult to analyze numerically. 

A more recent related work is that by R.H. Miller and R.P. Thompson 
[52]. Their work is concerned with the Pacific Gas and Electric Company 
hydro-thermal system. A linear programming approach is used for solving 
the long range scheduling problem. A set of inequality constraints on the 
reservoir's storage and head variations are imposed. However, the time 
delays of flows were not taken into consideration. 

It is evident that the need still exists for an optimum scheduling 
scheme for common flow systems. In the following two sections two distinct 
problems are discussed. In Section 5.2, the problem of a system with a 
general number of hydro plants on the same stream is treated. A more 
general situation is discussed in Section 5.3. Here a system characterized 
by multiple chains of hydro-plants is considered. In these problems, the 
time delay of flow between plants on the same stream is taken into 
consideration. Also, the effect of the tail-race elevation on the effective 
hydraulic head is considered. 

The time-delayed control systems field is still an active one. One 
of the main references on this subject is [53]. In this reference a dynamic 
programming approach is used to obtain the optimal solution for some certain 
problems. For an extension of the classical calculus of variations results 
to control problems with delayed arguments, reference [12] is an important 


contribution. 


5.2 Power System with Variable Head Hydro-plants on the same stream. 


An extension of the results obtained in Chapter 4 is made here. The 
problem of a hydro-thermal electric power system with variable head hydro- 


plants on the same stream is discussed in this section. The time delay of 
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flow between the hydro stations is taken into consideration. As the 
effect of tail-race elevation on the effective hydraulic head is 
considered. The results of this section were reported in [54,55]. 

0.2.1 Statement of the Problem 

An electric power system with m-thermal plants and (n-m) hydro-plants 
on the same stream is considered. The hydraulic part of the system is 
shown in Fig. (5-1). A prediction of the system's power demand and water 
supply is assumed available over the optimization interval. The problem 
is to find the active power generation of each plant as a function of time 
under the following conditions: 

1. The total operating costs of the thermal plants over the optimiz- 
ation interval is a minimum. 

2. The operating costs at the ith thermal plant are approximated by: 
zi 


pyre 8) = aq. + p,P. (t) +; Yeh es ty) $/Hr (522 a) 


: i i 
3. The total active generation in the system matches the load plus 


the losses. 

4. The transmission losses in the system may be represented by the 
general loss formula. 

5. The time integral of water discharge for each hydro-plant is a 
prespecified constant amount. 

6. The effective hydraulic head at the ith hydro plant is equal 
to the difference between the forebay elevation y; (t) and the tail-race 


elevation yy (t), thus 
j 


h.(<) = y,(t) - yy (t) 1S Mees san (5.252) 


7. The fovebay elevation y;(t) is given by: 


Lg i 


: nant wmataye itd 40 so ottuerbut att, 


eta bis bosinad “aWon aimee oad 0 ni 


ok 


7 nT Tavaaiet nofdesinisao, | 


“simttgo sid Yev0 atneta ‘femiads edt to ‘et2e3. 


is 


! ae 


ce ah fe * . ogee | 
vie, . a a. , Rae 


Rr 


i 


e p 
Bee 
my 


edt xd bainazeygen od ‘et n nataye att at 2922 


foups 2t ict 3 tee ast ait ts 
 8964- “Tis ort bas 1 Gh ae te 10 
fi ‘ Pia tre 


(S62) 


Figure 5.1 
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; suye ot By s.(t) i = m+l,...,n (5.2.3) 


This relation is true for vertical sided reservoirs. Yas and 8 are 


1 ys 


constants corresponding to the forebay geometry. 
8. The tail water elevation varies with the rate of water discharge 


according to the relation. 


yz _(t) =yr + By q(t) 1g= ait le. 3. ,n spartan) 
j jo j 


Yt and Br are Known constants corresponding to the tail-race geometry. 
j 


10 
Thus substituting (5.2.3) and (5.2.4) in (5.2.2), the effective head h, (t) 


is given by 


h.(t) yarn By 34 (t) - By a; (t) eve. ee. a Wore oF 


where 


Sere eae i = mly...an (5.2.6) 
9. The reservoir's dynamics are described by 
Smeg (E) = Fineg CO) + naga Cty aq) ~ Spey ft) 
i = 2,...,(n-m) (Sez) 
stp it) = ire) - Ine 6) (59230) 


Here, the time delay of water discharge between two consecutive hydro-plants 
is assumed to be a constant tr. 
5.2.2. A Minimum Norm Formulation 


The object of the optimizing computation is 


Te om 
Hing) ie F.[P. (t)]Jdt (5.2.9) 
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The generation schedule sought must satisfy the active power balance 


equation. 
(t) =} 
Pett) = se eae) Beant) =" Pet) (O22 ono 
ele ae : 
The transmission power loss is given by: 
ion n 
P, (t) a 7 P,(t)B; sP.(t) wee BP; (t) - Kio (5.2.11) 


Furthermore, the water discharge at each hydro plant is to satisfy the 
following constraints on the volume of water used over the optimization 


interval. 
If 
q.(t)dt =b. diate ecues ol (Se2et2)) 
O 


The ith hydro-plant's active power generation P(t) is given by: 
j 


P, (t) = ~—+— i = m+l,...4n eanic) 


where the G.'s are the efficiency constants of the hydro plants. Substituting 


(5.2.5) in (5.2.13) the hydro powers are given by: 


Gag q (t) 
: oy: + 8B ; ie (t)- 8, s(t) ] 


t 
Dye (t) = Spaj(0) +f iggy (Ide (5.2.15) 
A O 
Yes (tery) =f Gq (0-74 Id (5.2.16) 
t O 
Qnty (t) -| Anti (odo (ore) 
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Then integrating the reservoir's dynamic equations (5.2.7) and (5.2.8) 
yield: 


sath mbt q(t) - 04 (4) (5.2.18) 


then (5.2.19) yields 


meee) = Vag yh tot yy) = 9 
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(5.2 2a) 
and for the (m+1)st plant: 
Ge (t) re Olli te uel = ke eo 
Ph CP) ; Co] yore m+] m+] 
m+] er 
Rar amt) | (5.2.23) 
Cal m+] 
Let 
ee ae * Toon eee (Ney) 
Ansa (t) = = gay Cones * Ay. Dua (D1 


85 


| @ a fate pa (vs 
| a 
« ! eer: a ae - 
ves tg) faa 
. 2 ‘ ag | 


a 


ss ; ‘ E ( ST ; an. OS: aah By 7 jh f | 
| , iby a “ aH " x Le) in 7 7 =-s a 
e | es Ove ee iy eM RD = an 
’ y 


, Be ee ae Aig 
mo 


bree. Set 


8 i 2) 


‘i | co 


: A us . h. P 


86 


y 
+3 , . 
A ie fe tes eno) tS. 2225) 
m+7 
B 
ine 
+ 
Oe alam i= 1,...,(n-m) (5.2.26) 
Cn 


then the hydro power given by (5.2.22) and (5.2.23) becomes: 
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The problem now is to minimize the cost functional given by (5.2.9) 
suppect tO Samistying. (512.10), (5..2.lO), f5.c.l7 ta (See. cos (5.2.21) and 


(5.2.12). The constraints (5.2.21) are rewritten in the equivalent form 


cee ec 


i-l 7 27 2.55,.40n-) (once cu) 


Let 


1) 1 Sys (-m) (5.2.31) 


hee ee? | 
RA eee th : 7 aS 
ite . 7 : 

tie | a" i B te 
; i ; 
a \ Te 
Wy ne us 
4 ‘ t* t 


kn 


Pl’ eee 


i v e 
i Wy ict 
4 t ; es i 
U ne 4 me 
% SG Toy ee : iy ES ook Alea hu Aon a ey wu 
i i dir 1 te eR ey 4 j us _ ie 
t ie i a =. hy A Bliss La 
i a La, Aik Nelo) ee 3 
1 1a ¥ i a i uae 
&, an i a 


cond (65-2) be GS-8) Ww navi Seah: 


te 


8 re aka MVE yin iiike eat oh 


ee Peel " cre, (iad - (a). wit’ 
‘ ah nS wry ean RO ne CS ee 


Eis atk a A ie 


| | | ee ee ef 
: (S$ .S.2 0 (}) pm? em 7 
; er tpn” 


I aN Ye qh — at) AON hain re ee wD: = (3) Me epee Oe, ee 
(8545.2), (m Meereed as) Pia ee saad 7 ant 


it 7 Ty - 


| ‘ a aE amas 7 Udo aay dia 
(O08 8) yd navtp jenoftonut feo ot osiitintin os 2h on otdona 
eae ee es ohn ae aD aie MH eas 


| au) ei png ai be : 
big (18.8.2 i +(8S.8.8) (SES, ) + (3F. ay on &. or 


PN tel a 
nati tape ‘ate USS 2D etntendine 


Meee 
Ane 
a ie 


mot trate ups "ot mt! 


me 


1 MRSS) ATA sinlese UB 


; us y : _ L : , i , (ay is mS 
oa : ; : “x : ~ ea: if 
| er ae : om a t a7 eT ee Ace y 
7 7 : - : q wy . At J yy f 
an y : : E . ie Ee ie 
i 7 a i ; . A mn The : 
4 ry en i * i » t, : 7 i 
n 1 5 : or, _~ i foe 
ye ‘i om ; wh a 
_ i VY ; oy aa 1 , 
ei eo, oe 6° 
\ ; i y _Y ty % ° 
r ‘ ¥ 
: 


a | | } { ; 
7 (ob.8.2) asi 


a 
yr) 


87 


Note that Vnti-p6tsty_y) is a known initial condition function. Then 
(5.2.30) becomes: 
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Moreover the constraint (5.2.17) is rewritten in the equivalent forms. 


Gait) = Ot) (5.2.34) 
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An augmented cost functional is formed as: 


6 
J(U(t)) 8) J. (U(t)) (Ove.00) 
i=o 
where Jo is given by (5.2.9), and U(t) is the control vector defined by: 


U(t) = col.[P(t) Wo, (t)Wayo(t)s--- W(t) I (5.2.37) 
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a em m+ 
(5.2.38) 
Wagqdte) socole(quys (t)s0ee,(t)] (5.2.39) 
Weg (t) = col.[q.,.(t),Q...(t).x..(t)] 
i = 2,...,(n-m) (5.2.40) 
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where A(t) on 4 (t)s m44(t) and Key (t) are unknown functions of time. 


‘ 
Jy corresponds to the power balance equation (5.2.10) and (5.2.11) 
after dropping Py(t) which is independent of the control vector U(t). 
Jo and Jy correspond to the constraints (5.2.27) and (5.2.28). Ja and 
J, correspond to (5.2.34) and (5.2.35) respectively. Finally Je corresponds 
to (5.2.29). Note that in (5.2.45) and (5.2.46) equivalent forms of the 
constraints (5.2.17) and (5.2.21) were used so that a valid norm definition 
can be obtained. 

It is more convenient to reduce Jos J3 and Je to modified values by 


performing integration by parts and neglecting terms explicitly independent 


of the control U(t). Then the following is obtained 
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i | 
f n-m D 
Te ues Y Tants (Engg tot Dat (5.2.52) 


Interchanging the integration and summation operations yields 
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(5.2.59) 
Integrating the last term of the integrand in (5.2.59) and neglecting terms 


explicitly independent of the control, one obtains 
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Supstitucing, (5.219), (5.2.41), (5.2047), (5.2.48), (5.2044) (52-49) 

and (5.2.60) ‘in: (5.2.36) for Jas ij = 0,...,6 respectively yields 
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and the square symmetric matrix B(t) as: 


B(t) = diaglB,(t) Baay(t) Bpyolt)>-0B,(t)] (5. 
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Hy (t)  Bagandtt) [1 - Big | (5 
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Bo aght) = Mei lt) = ACEDLT ~ Bega) (5, 
2 met) q 6) A (tiny (t) + may it) eS) 
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Agee = Ay tein te) (5 
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Oe a ee Ta 
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tt 
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Boag ft) = (by; (t)) i= 
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Bich) ts nnd V1 on 
M43 (t) 


m+j m+ 
th (t) 
be“(t) = n— +r (t) 
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This reduces J[U(t)] to: 


Let V 
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= vey B(t) Wei at 


The last term in the integrand of the last expression does 


explicitly on the control U(t). 
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f V V(t | 
Jfu(t)] | tu(t) + Sey "B(e)ru(t) + = dt (5.2.88) 
0 
subject to 
T 
i nt 4 (odo = bei nj elt eee aes) (5.2.89) 
0 
Define the (n-m)x1 column vector: 
eo Aaa (5.2.90) 
and (n-m)x(4n-3m-1) matrix Kas: 
T 
Tt 
K ah 
Ft 0 0 Ses 
KI 
—~n 
(5.2.91a) 
oe 5.2.91b 
with Kee Ort Tila 0] ( ) 
lal © 5.2.91 
‘sate aan bali ( c) 
This transforms the constraint (5.2.89) to: 
I ¢ 
b | klu(s)ds (5.2.92) 
0 


The control vector U(t) is considered an element of the Hilbert space 
A eeu te of the [4n-3m-1] vector valued square integrable functions 


defined on [0,1] endowed with the inner product definition: 
Te 
t 
<V(t),u(t)> -{ Vv (t)B(t)u(t)dt (5.2.93) 
0 


ary a 
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hen ae 
| s{@.S.8) 
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motion stasigsaat ‘sno 


97 


for every V(t) and u(t) in Lenn. T I, provided that B(t) is positive 
definite. 
The given vector b is considered an element of the Real Space pin-m) 


with the Euclidean inner product definition: 


SK aca Oe (5.2.94) 


for every X and Y in Rin-m) 
Equation (5.2.92) defines a bounded linear transformation 
[re ae rk This can be writen as: 
enna (5.2.95) 


and the cost functional given in (5.2.88) reduces to: 


Itu(t)J= |Ju(t) + 4)? (5.2.96) 


Finally, it is necessary only to minimize 
V 
Itu(t)] = |fuce) + =e), (5.2.97) 
subject to b = T[u(t)] 


for a given b in p(n-m) 


5.2.3 Ihe-Optimal Solution 


The optimal solution to the problem formulated in the previous 


subsection, using the results of Chapter 2 is: 


u(t) = T'[b + rtd] : seh (5.2.98) 


where T" is obtained as follows: 
ee the adjoint of T, is obtained using the identity: 


* 
Seow nem Lew (4n-3m-1) 


Loop [0,1] (5 52(99)) 


Let 
& = COl.[E  yo+++ 28nd (5.2.100) 


* 
i lei aig col.[T, sTingy tme2?* + nd (522.101) 


, ¥ 40 tndrato 0 ns “bovobtane oy at ‘i 0: 2 
| | nor statteb soon 9 


anita (se. 
ne: 
a | ; hd : 25° nod ry od) 1169. ai i ~ nit ‘ 


Pcabede 


ete ais 


ei ; . x y “fo sala ete at sesaitt baby pork & 3 


ae fn 


. sot 2 
(heeS.8)) 


Faso 


by 
os a 


auofvaig oid mf batsfumiot motdorg 


whe = col[t,,...,t,] (5.2 al02) 
T = col.[t ot ] (5.2.103) 
aie nt] nt] 
Tee colt Ale ae ] (5.28 104) 
ccs nti mei Xmed 
In pin-m) the inner product of the left-handside of (5.2.99) is: 
ai 
<€,Tu> = a Eshie: (5°2.105)) 
0 
This reduces to T 
n f 
<E,Tu> = ) E. | q. (o)do (5.2, 106) 
i=m+] 
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In Cet ot nogt 1, the inner product of the right-handside in (5.2.99) 


ie 


. 
‘ fe 
<T’Eu> att (tc) "B(t)u(t)at (5.2.107) 
O 


Using (5.2.37), (5.2.67) and (5.2.102) this reduces to 


<TE 5u> aha 


p(t) Idt (5.2.108) 


Thus the identity of (5.2.99) using (5.2.106) and (5.2.108) yields: 


1 T 

fn-m f t a : 

- P(t) + a pe lh 

JE, Pein (280 fo tte cae) eT Baas 

oa 2 
Wag (t) Jat (5.2.109) 
Let 

ae 5.2.110 
Omt [oe] 0] ( ) 


aah (ors 
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ne Re Doe ae 


ae al 
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ae io Daria a i= 2,...,n-m (2d 
then (5.2.109) becomes 


Tene ‘ ly T n-m T 
| - 

j oy oa gtines ode [ taper +2 Titi Bned () 

0) 0 

Wo4 (t) dt tbc lz) 
Equation (5.2.112)-is satisfied for: 

T(t) = 0 (5.2.113) 
=p ease} 

T ig] . 

Te (Gb) sd Bi ee oe (5.2.114) 

m+4 inthe 


From (5.2.77) the following is obtained: 


ai j -] 
BGS) ssadiag](C., pmag (t))° sf-(de yo (tsny) 
h(t) 
m+] =] 
£8 eo hi (5.2.115) 


thus one obtains from (5.2.114) for i = 1 and (5.2.110): 


b 
I © +] 
ib Aantet hes [en tey 0] (5.2.116) 
Let | 
-| rule 


then by (5.2.114) and (5.2.111) one obtains: 


T 


= ° 3 " : b ] t 
Tae it) Ubi 11g (O Pein (t) ming. | )] 


1 = Zac cea (5.2:418) 


This means now that TE as given by (5.2.101) is completely defined 


at “103 vorietvee ab {800.5.2) 08 
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Oy (5.26113) si. 225116) sand a(5.2 11a). 
The operator J is next evaluated as: 
J(e) = TET] (5.2.119) 
Using (5.2.91), (5.2.92), (5.2.95) and (5:2.101) this is given by: 


Jf Vs 
1 
J(é) = col.[b | dt,b i a (Cate ...29) 
= m+] CT me] Ot) m+2 Ve 
’ : (52120) 
or 
onan (522, 12)) 
with 


8 ie 
‘ ] 
A = ciagl { ay tt | a (tidty..<)) GeZe.12Zp 
2 Cot ne] t 3 Vinee 


Thus the inverse operation I is given by: 


PT Pw! & (5.2.123) 
where i" 5B | 
it = diaal( | LC ea t)]dt)” ay a dt) leeeane 
0 0 


(Be Ze lees 
Finally, the pseudo-inverse operator T’ is obtained from the 
definition: 
IE =) Leg] 


Ucing Aoe2e101).- (52.113) and (5.2. 116) this: becomes 


y Sa i (5.2.125) 
fee seco. | tye Pies ed 
ire pW] Wh 

with 
t! = col.[0,...,0] (5.2.126) 
4 b 
+] 
on = col [ a 9 0] 
+] 
i Goer (8) ye at ey tia 


25" ‘ie at ive tx ‘ a iy 


Nene a. AN mY ioe : 
pia Nd nowte ah efi F (F9t.s § & | one (8 


«enol at it 


Let 
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ra : ale ; (t) 
“Wo. = col.[t! t PUN Uae eee 


(mi) ,."(mi), a) Ce) 
(m+i) 
ere 


e m+7 ) : 
t . po | T= -2505 
(m+i) 1 ayy t 


(m+7 ) 


b ty) 


-a 
3 . m+i VV nei) 


m+ Te 
i a] (t)dt 
vs Mies hele 


[5 
il 
|o 
+ 
—-f 
— 
MII 
we 


a cOl.[niyysss- on] 


then applying (5.2.92) one obtains: 


T 


f 
: 1 
Mn > mts “| 2 (mei), (dt 
6) 


Pouce 


pore) 


(oes 


where it is assumed that V(t) given by (5.2.87) is obtained in the 


partioned form: 


This gives the 


form as: 
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(52.137) 


p(t) = = W(e/2 Papers: 
Y (mer), ©) ml 4 a = 
¥ 0 
Feat) = TE Aig ET Peg Ain Ty 
mae [1/n,,, (t) Jat 
‘ (§.2.138) 
Y(me1)6¢) 
Q(t) = - pk (5.2.139) 
Re i: U 
ty Oe weer ) i huis nh ) 
m+7 iP 
i (tat 
Aa) 
me (5.2.140) 
ia Mm+i “12 me (t) 
mas fos 2 weg 
eee) : ie 
fan (t)dit 
a (m+ ) 
ieee orn (5.2.141) 


ty abe) ep ae (t) 
x Gea) = - Wp, Set) (5.20142) 


5.2.4. The Modified Optimal Solution 

In formulating the problem at hand, pseudo control variables q(t) 
and x(t) were introduced. It is possible to eliminate these variables 
together with the multiplier functions m(t) and ratt) associated with them. 


First we establish some relations between some of the variables encountered 


in the previous subsection. 
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Using (922.66) ,4€5.2.72) Athi 27eys (5.2.74 iy (592887) and CBee) 
in (5.2.136) we obtain: 


monet), S Aes neg e)* Paiin ergdang te) 
(iS 2S wee (5.2.143) 
2.40 
V(mei),‘*) ; Momei), aay m 
ie e2 5 eration (5.2.144) 
a3 (t) 
Vomti)6®) = V(meiy, &*) aya 
(m+i) 
i= 2,.0.,n-m (5.2.145) 


substituting (5.2.144) and (5.2.145) in (5.2.141) and using (5.2.140) 


respectively one finds: 


oe ee 
ae ee Liza F) 9, ,(t) i = 25...4n-m (5.2.146) 
m+1 VV mei) m+1 
seule? 
xe (t) = A ae ft) 1 = 25eeesmem (5.2.147) 
mt+j TT ned) m+1 


Invoking the constraint equation (5.2.21) using (5.2.146) and (5.2.147) 


One obtains: 
q (t) 
t emti 


1 = 2on cane (oscel4e) 


From (5.2.78) through (5.2.85) and (5.2.117) it follows that 
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0. $2) “te sa ‘ 


me (8, §. a 


(2pf. 8.8 TN eae 


(ODT.$.2) onteu brs: (BE. S8ua 


: . } WL \ Le , aN, ie 


Pee jh el eel aan ie 


({P1.9.2) ‘bos (abL.Si¢ 


104 


md. | 

ie Ne ‘ rE bmi ey (boty) I/dy g(t) 
12202)).0 nem (5.22. 149) 

mh (yj ee 
V2 (mei) mm é mimes E) 1/2445 (t) 
(2 nem es ose 50) 
ty eee 
Ahaha = [reg (tim g(t) + Ban. (t)( ne + 
re Ge) * dmeie] (tots ) 1/24, (t) 


1 2523 ine) S22ul5)) 
where bees ht) is the determinant of the matrix Bag (t). Thus (5.2.148) 


is rewritten as: 


Bienes t) 
Y <q(tst,_7) = a(t) (S721 52) 
m+7-1 i-] E 

AR) 


and (5.2.146) as: 
Ene (t) + Bei Mined (tJ 


(Cor west oe ea Te q; At) (5 2aubs ) 
mi i (c) coma es Canoe Ts 


Substituting (5.2.143) in (5.2.140) one obtains: 


2q, _ (t) 
m+1 a 
ayy (t) e Aq (t) May; (t) as Preity tots) Cati (5.2.154) 
mt+j 
with 
mi iu (5.2.155) 
m+i 
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The expression for ayy (t) from the inverse operation given by (5.2.149) 
m+7 


after substituting for A_,.(t) yields: 


m+7 


m+7 m+ 


226 Sone) = 
a7 (t) m+7 mt+i Dh ok) = ae 


(mg(t) + Boyan, .(t)]° 


meity (boty) 


iZ 2 
Basti mit) 
. 2 e e 
Ar 44 ct) (5.22156) 
Substituting (5.2.156) in (5.2.154) and using (5.2.152) and (5.2.153) 


we have 


2CneiMmey (td, Ct) : And (EI Ot) ‘i Pintity (t2t 5) 


“Beams (t) Yinei-1Et4-1) * Um g(t) + Bae ghey (t) 


eat Nay (Discarloze) 


Let 


(4) = 2 estes (EG, Ct) + Aggg Ct) Mg (t) 


= 8B be) ued r ey) (5.2.158) 


then differentiating (5.2.157) and substituting (5.2.158) and (5.2.57) one 


Nn : 
m+1 m+ 


obtains: 


(+) die tne o(t)Q,  (t) #182 mo tt)a, CH) 


m+i mt 


d 
ek a 20mg (TG 9T4 inti (ttTy) = 0 


m+7 
tex Te - 1; (5.2.159) 


and 


ny 


106 


(Co) FEB ene eG 


Gea uy Bat itmey (t)Q eaten 


d : 
+ dele (IQ, (tI paOh a Pe eae 7 12.160) 


Substituting (5.2.56) in (5.2.153) and using the result in eliminating 
Me (t) from (5.2.159) and (5.2.160) these yield: 


SIM en li blir pelea (5.2.161) 
and 


Loewy er Bed neg (tI On gy (t) il Ope woah Te “Ts (3c Or) 


Rewriting (5.2.152) advanced one plant and with time lead T, as: 


Seay | eae (ttt. )q (tt+r.) = ZV 2 (t+r. ) 
mt+i+] m+it] 1 emtit] 1 m+i+] i 
[Une g (Te ots) + Ceo (5.2.163) 
One can eliminate Yr aataayy (ttt) from (5.2.161) to obtain: 
Zeit) + Bang hmeg (EQ, Ct) + Bra ieaMneiat (Pty) 
q- (ttr.) = 0 ele eee (5.2.164) 


m+i+] 


Using (5.2.158) in (5.2.162) and (5.2.164) one obtains: 


Fel 2C pases (40 CE) eee ilnette — met nea 


emt 


(toty y)]t Bayh g (OQ, Ct) + Basia nape (ttt) 


‘Een gaa! ats Ri aed (5.2.165) 
n+it+l] 
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and 


Geen (80, te) t Ang CEng CE) - Bane 
mr] 


Pera ave ene Uae ty) Wren atic. (t)0 


mints ieee 


m+ 


1c e(Te ~ ts Te] C5n2e L66) 
The equations (5.2.165) and (5.2.166) thus obtained depend on 


n_,-(t) and Q neq (t)- The multipliers m_,.(t) and r_,.(t) were eliminated 


m+ m+} 
as a result of invoking the corresponding constraints. This holds for 
i = 2,...,n-m or hydro-plants affected by plants which are upstream of 


them. Note that for i= n-m, the terms inr do not exist as in the 


mt+i+] 
cases i = 2,...,;n-m-1. However, it can be shown that (5.2.165) and (5.2.166) 
hold for i = n-m by repeating the same analysis. 


FOr ="1.- asing™ (5<c.70), (o.c./71) 5 1(0«2.77) and (ove.0/) one 


obtains: 
Vane yy 6) = [Ang (tng (t) + mia (t)I/C.. 39,45 (t) 
(Soc. ten) 
nel 9)" = = [ih y(t) Bago (tsty)I/[o,49(t>ty) 
bh alt 
+B Petty (5.2.168) 


Thus (5.2.138) and (5.2.139) become: 


Crete T(E, | (t) : Any ED Pippy CE) * Mey (t) = Oney 


(5.2.169) 
hoCh epee ce) 

Q t) Se Oe ee aC) (5.2.170) 
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Differentiating (5.2.169) and substituting (5.2.170) for hey (t) one 


obtains 


rg cretee Ele CW) armen ce 


E 


m+] 
: Petey (EQ, Ct) bi PeaMmeg(ttr)Q, (tery) i. 
a e[0,T, -ty] (5.2) 
BELA (EG Ct) i Any tng (t)] 
Peas g()0C eg (ENO) tele =n ai mea ono IT) 


gE 
Note that in (5.2.171) and (5.2.172) elimination of rae (ttt) was made 
possible by utilizing (5.2.153) for i = 2. 
We are now in a position to rewrite (5.2.165), (5.2.166), (5.2.171) 


and (5.2.172) in a unified fashion as: 


d . 
Fetes Mmei (tO, Ct) Heuer coliensy G2) 


+ B ( thong tq (boty etG_7) =" 


med ea 


= hanna (Sez 70) 
with the boundary conditions: 
Q. f= 20 7 245 oe 0c (502.474) 
m+i 
et aig ge i-= 1,033 nm (See kes) 


where 
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oy at 


ie th 
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Gmeq(toty) = BapoMneg(ttty a, (tty)  t eL0,Te-ty] 


m+2 


=y | 1b elTe-tysTeJ 


(552.126) 
Ge § boty oth) = Beier neiey (tty )9 (tii) 


Smtitl 


d 
= dae eines Mg eYey-phtet yd 
ne el0.t, 1] 


= 8B (ttr.) 


(tt. )q 
mit] m+i+l emtitl 


d 
F de One Mined SE) CY a (Gd eTG 1 


£ e(t._1.Te-t;_7] 
Seg meet neg CE) Vinay 784-7241 


tiQneveadtaty titles shee estate 


fie Pee nated (Cerra e) 


dil taze cae ton fal M(t) op tot mds 


secfOse ) 


n-m-1] 
Bn (t) vp 4 (Tp me] ® nome) 


py t lta mle 
(5.2.178) 


ate 
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It is interesting to note that (5.2.173) describes the optimal 


interaction of the volume of water discharged Q (t). The solution of 


E : 
m+ j 
this equation will depend on n_,.(t) which is to be determined such that 


m+4 
(5.2.28) is satisfied. The value of Ney it) depends solely on the (mt+i)th 
plant characteristics. Moreover, the solution will also depend on 
Gtq (tet; oTs_7)- This depends on the behavior of the plants immediately 
upstream and downstream from the (m+i)th plant. 

The equations (5.2.137) and (5.2.173) will be referred to as the 
modified optimal solution. It can be shown that this solution satisfies 
the necessary and sufficient condition for optimality for control problems 
with time delayed arguments as given by D.K. Hughes in [12]. 

9.2.9 special Cases of the Problem 

It is the purpose of this subsection to show the optimal solution for 
some specific cases of the problem considered in this section. First we 
give the solution for the case when the water flow time delays between hydro 
stations on the same stream are negligible. The results obtained here are 
identical with those obtained in [7]. It is worth mentioning that this 
special case was formulated and solved without introducing the pseudo- 


control X44 (t). This is possible since x_,.(t) was introduced to facilitate 


m+j 
dealing with the time delays. The results obtained both ways are identical. 
Next, the solution when tail-race variations are neglected is given. 

Finally, the optimal solution for plants not on the same stream is obtained. 
It is shown that the solution here is identical with the results of Chapter 


4. In all these cases we are dealing with the hydraulic discharge equation 
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A. Negligible Time Delay Case 
In this case wie ‘Orsorthat (a2, ST) yie lds 


t 
Vint -1 (t 20) = [ape j(od0 t $ 0 
10) 


and 
Vnti-1 (020) = 0 (5.2.179) 


Also (5.2.176), (5.2.177) and (5.2.178) transform into: 


She (90) = Brgatnga Ct), (t) (5.2. 180) 
Trt E2020) = Breia ames (tag Ct) 
- TB, Mei (#0, (1) 
195h25..5-5(NeM=1) (5. 2210)) 
g,(ts0) = - $e{B,n,(t)Q,  (t)] (5.2.182) 


Thus the optimal discharges are given by (5.2.173) using (5.2.180) through 
PS cae) 


B. Negligible Effect due to Tail-race Elevation 


tw this’case C_.. = 0, since 6 = 0.. Thus (5.2.173) is rewritten 
m+ es 
for this case as: 


d_ 
Q, (t) = - aaa oe ea oT. at dt 


m+i m+i mt+i 
ee One Uae oe a 3 
(Angg(t)magg(t))J i = 1,...,n-m  (5.2.183) 
where the g_..'s are given by (5.2.176) through (5.2.177) in the case when 


time delays are present or (5.2.180) through (5.2.182) when time delays 
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are negligible. 
C. _Hydro-plants not on the same stream 
In this case two observations are made. The first one is that in 


(5.2.7) the coupling term q (t-t._)) will be absent. The second 


mt+i-1| 
observation is that one can assume that the T.'S approach infinity, so 
that the discharge from one hydro-plant would not effect the others. 
These two observations can be stated as: 


lim q 


T .>00 
1 


ec) = 0 12-1, nem) (5.2.184) 


m+ 


This leads to: 


DT a? = 0 1.2. 190; .gdnem (5. 20185) 
and 
Vien Uapeg bot 4) = 0 i= 1,...,(n-m) (5.2.186) 


F 
The last two equations result from applying (5.2.184) to (5.2.16) and 
(5.2.32) respectively. Applying the last results to (5.2.176) through 


(5.2.178) one obtains 


Gaim ohn thea a? =e 0 it) lsieeshen (5.2.187) 
F 


Thus (5.2.173) reduces to 


+B (0 6 (ey Sis a a een (5.2.188) 


mi nti 


Furthermore, if tail-race elevation is neglected, then (5.2.188) reduces 


to: 
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1 = 1,...,n-m (5.2.189) 
It is our intention now to show that (5.2.189) is identical with 
(4.4.97). To this end we rewrite (4.4.97) as: 


q(4) (t) = meg OE) 


ae ee) (5.2.190) 


and using (4.4.20) this is given by: 
(t) ‘ 
gy tt (t) = IN, ys (t) ey belie) (5.2.191) 


m+7 Emi 


Also by using (4.4.10) and (4.4.11), (5.2.190) reduces to: 


t 
a4) (t)=s, oh, (0) + Jina; (o)de 
6) 


>m+i ‘nti mi 
44 (t) 
Hee ins 
Note that S; here is that of Chapter 4, which means that oF h. (0) 


mt+i =mt+i 
is in fact the initial storage Si 44 (0) of this chapter. Thus we have 


nee) 

4 mt+7 : 

gee) tt ) = [Spe CO) * Tye (OD) - peep tes Ct) 
(5.2.192) 
Using (5.2.15) this is given by 
het) 

4 _ m+7 : 

qe) Ct) = Day it) = heat D4 (t) (5.22193) 


rv) Sees) ees (5.2.194) 
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Observe that in Chapter 4, On; Was taken as zero. Thus (5.2.194) and 


F 
(5.2.193) are identical. This proves that the results of section 4.4 
are a special case of the results obtained in this section. 

9.2.6 Implementing the Optimal Solution 

The modified optimal solution obtained in (5.2.4) contains the 
unknown multiplier functions n. (t) and A(t). These are to be determined 
such that the constraints (5;2.10), (5.2.27) and (5.2328) are satisfied. 
In this subsection, the reduction of the equations specifying the optimal 


schedules for the general system is given. 


Let the (nxn) symmetric matrix: 

Brit) = (b.. (t)) (5.2.195) 
be the inverse of B,(t) as given by (5.2.76). Note that elements of 
B,(t) are dependent on A(t) so that one may write 
HL GN as) (5.2.196) 


From (5.2.64) and (5.2.196), the expression for V(t) given by (5.2.133) 


and (5.2.87) in the component form (5.2.134) is: 


(5.2.197) 
and the optimal power generations given by (5.2.137) component-wise 


are: 


Pp. (t)=-V_ (t)/2 Ta tren cree (5.2.198) 


Prequbtidig: = WCE) 2 i = mtl,...5n (522-1039 
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Here it is noted that the power generations depend on the n, (t) and A(t). 
The optimal power generations should satisfy (5.2.10), thus the 

following relation must hold true: 

n 


n n 
Pyit) +5 } Vp. (t) + a) Vp, (t)B; sp. (t) 


i t=0 gel P; 


= 0 (5.2.200) 


This is an algebraic equation in A(t) and the (n-m) functions n.(t). 
The optimal hydro-generations should satisfy (5.2.27) and (5.2.28), 
thus 


V (t) = 2[A 


Wt )Gaa it) +)Bos Qe (PIG. y(t) 


ae 
tae y(t) (5.2.201) 


Peds t2EA_,.(t)0...(t) +850 (tO. (t) 


- B 


(t,t, 1) 


ma teeny Cree 


$C QM. (t)] f= 2ynlinem (6.2.202) 


m+7 
Here Vatie](teta_7) are given by (5.2.32). 

The problem of implementing the optimal solution is thus reduced to 
that of finding [2(n-m)+1] unknown functions Q,(t)[ism+1,....n], 
n.(t)Li=mtl,....n] and x(t). These are obtained by solving simultaneously 
(ee2.173), (5.2.200), (5.2.201) and (5.2.202). ‘These are exactly 
[2(n-m)+1] equations in the unknown functions. Due to the nonlinearity of 


these equations, one inevitably resorts to iterative techniques. 


Consider equation (5.2.173), rewritten as 


my 
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+ > b) 
mt+j ios | 2C a 2C a4 ni enti), 
1 = 1,...,n<m (522.2033) 
Let 
Deere) 
_  mti ss 
Png it) = nie 7 = 1,...,n-m (5.2.204) 
Ly fils Fe Repaid (5.2.205) 
. = = hb, n- » bye GUD 
m+} eC 4 
Gyq (tots Ts tod OE (t))= Gey (tots stay) 
A ..(t) at) 
+ oF + aie Pq it) i= 1 ,c.e6n-m (5.25206) 
m+7 m+i 


Then (5.2.203) reduces to: 
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fait) = Ri s(t)Z a(t) + & (t) 


1q-ald.... (Hea (5.2.209) 
with 


0 ] 
Rei ft) = 
Emed mei (ED “Omg (t) 
1 to »(n-m) 
LG, ‘ col. (0. “Gi i = Lee ih) (S220) 
Here the boundary conditions Q (o) = 0 and Q (Tz) “S9pe et canebe 
omit j Smti f per 
written as: 
M Zmeq 60) +N Zag Tp) = Ge 
ae ee eri (5.2211) 
with 
10 ot | 
Mé=8/8) pa) (5,215 7) 40 ob (5.2.212) 
C+ = col.[0,b,,<] ilcieall bee (5225293 5 


Equation (5.2.209) with the boundary condition (5.2.211)are next transformed 
into an integral form. The results obtained here rely heavily on those 
given by P.L. Falb and J.L. De Jong in [56]. The following is an integral 


representation of (5.2.209) and (5.2.211): 


‘ 
t 
bint s[T.-t] ft[T.-s ] 
] m+1 f as apa 
261) (+) = =, sf coats Foal Cele +| Ty Fney (S)ds 
t 0 3 t 


| Dy bennotentns #xon snstits gu ecblaiti visbad 


(80S.8.2) | (MEAP heat 
‘ rR | a ay 
ooh et cue CF) ss 
phe : . vT: ‘ 
fortdesanel St 


(O7E9:8) | Reo. et ee F Caaf 0 Feo = =e 


Pe i i 
a pa 


gt) 30 bars 03 (0 ) : Pe nots tho sind a 


oa 


oo 
= 
— 
EA 
- 
if? 
. 
Ta 

tee 
SS 
=a 

a 
> 
— 

ii 


ieterorae y's! 
“wis. , ieee A ehh 
. Nay Gan Ok Se iplan 


th a ee ae, ee 


"Seon no > vi tse vier ont bantaste | esteem 


“AlSS, 2) 


my f 18,5 8) 4 oe : 


ii 
4 


118 


Here 
fine 5) = SypegPmnj (S)Zpe4 (8) + omnes (S)Z (a (5) 
+ Ga ts) 1g= "1a... (N-m) (5.2.216) 


Note that in (5.2.215), satisfying the boundary condition (5.2.211) is 
guaranteed during the search for the required solution. For the sake of 
Simplicity, the practical application of the suggested computational 
scheme is shown by way of an example. This example is concerned with a 
one-thermal, two hydro-plants on the same stream system. This is the 
subject of the next subsection. 

o.2./_ Practical Application 

A computer program was written to solve (5.2.200), (5.2.201), 
(5.2.202), (5.2.214) and (5.2.215) to obtain optimum generation schedules 
for a sample system. The sample system is characterized by a diagonal 
loss matrix (i.e. By; = 0; i #3). The number of hydro-plants is two and 
they are on the same stream. There is only one thermal plant in the system. 


For this system, the equations describing the optimum mode are: 
ds 
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t , 
3 
ZA P(e) = rfb3t + J stty-t1¢g(s)ds + | tLT--s]f,(s)ds] 
) Us 
(5.2.219) 
2) t T 
2, (ye T,13 + J. sf.(s)ds + | (Te-s)f,(s)ds] 
: (5.2.220) 
Define 
Coo = col. [y;(t) yo (t) s¥3(t) y(t) (5 2e22 0) 
where 
yj, (t)-4 241) (+) (5.2.222) 
y(t) = z57(¢) (5320223) 
y3(t) = 281) (+) (5.2.224) 
y(t) = z\2) (4) (5.2.225) 


Then equations (5.2.217) through (5.2.220) define the operator 

equation 

¥(t) = T(¥(t)) (5.2.226) 
This equation can be solved iteratively using a modified contraction 
mapping algorithm of the form: 

PD te wre (5.2.227) 
Here I is the identity operator, U is any operator which is linear and 
[I-U] is invertible. The convergence conditions for this algorithm are 
given in Appendix A. 


The initial estimate of the unknown variables is taken as: 


et) = b5t/T, (5.2.228) 
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SAMPLE SYSTEM CHARACTERISTICS 


Thermal Plant: 


ase 4.0 


Loss Formula Coefficients: 


oo 
The Hydro Plants: 
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Figure: 5.2 
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Optimum Schedules for the Sample System. 
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Figure 5.3 Optimum Head Variations and Rate of Water Discharge 


for the Upstream Plant 
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Figure 5.4 Optimum Head Variations and Rate of Water Discharge 


for the Downstream Plant. 
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Variation of Relative Error with Number of Iterations. 
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yy? (t) = b3t/T, (552-230) 


yao (t) gebasT (5.2.231) 


A sample system whose particulars are summarized in Table (5.1) 
was used to test the method. Optimum loading schedules obtained are 
shown in Figures (5.2), (5.3) and (5.4). Figure (5.5) shows the 
variation of the relative error between successive approximations with 
the number of iterations. The operator U for this test problem was 


Raker as U = -0.95q4 


2.3 Power System with Multiple Chains of Variable Head Hydro-plants 


This section is concerned with the case of a power system with series 
plants (on the same stream), multiple chains of plants and imtermediate 
reservoirs. The variety of models that can be considered from a theoretical 
standpoint is infinite. However, a practical model is chosen in this 
section. Here the formulation adopted is applicable to any practical 
system with a larger number of hydro-plants. The results obtained here 
are reported in [57]. 

5.3.1 Statement of the Problem 

A hydro-thermal electric power system is considered. The system has 
one thermal plant and eight hydro-plants. The hydraulic portion of the 
system is shown in Fig. (5.6). The operating conditions require minimizing. 


T 


J - [te Pe Gace? 2(t)]dt Gn 
: g g 


while satisfying the following constraints: 


1. The power balance equation. 
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Figure 5.6 General Layout of the Hydro-Plants for the System in 
Saction 5.3. 
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8 9 9 | 
Pa(t) = ie + i sire we iby Ps (t)B, P(t) 
9 
ee Bae Rais) (Sone) 


2. The volume of water discharged at the hydro-plants during the 
optimization interval is a prespecified constant: 


T 


f 
I q,(c)do = b, ere ck (5.3.3) 
O 


3. The upstream plants' active power generations satisfy: 
Ppa (t) + Ay(t)aj(t) + Bay (t)0;(t) + Cya,°(t) = 0 
leSuehee 535058 (5.3.4) 

These are upstream plants so that all of them can be represented by 
(o52-c3)% “Note that A.(t), B. and C. are given by (5.2.24) through 
e220): 

4. For the intermediate plants; the situation is illustrated in 
Fig. (5.7). Let there be (i-k) plants upstream from the (i+])st. The 
flow from each plant has a transport delay of ta(Jaks...57) to the 


(i+1)st plant. Then the (i+1)st reservoir's dynamics are expressed by: 
5a 25 
Integrating (5.3.5) one obtains 
S547 (t) =D, 47 (t) + ey x.(t) - Q. 47, (t) (52326) 


where 
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Figure 5.7 The ith Hydro-Plant's Reservoir. 
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t 
a(t) =[ 4; (odo (5.3.8) 
O 
iG 
x; (t) =| q, (0-1, )do ate ae (5.3.9) 
0) 
Let tera 
vi(tyt,) a h q,(s)ds t< a 
e j= kysceei (5.3.10) 


which is a known function of time from the previous history of the 


system. Then (5.3.9) reduces to: 


x s(t) = os (t.t5) pet, 


I aek se 1 (5.8.12) 


The hydro-power generated at the (i+1)st plant is given by: 
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+ Cay Qiay (EY + By ag [EIA CE) = 0 Sama) 


In the system at hand, plant number 4 is downstream from plants 
number 1, 2 and 2. Also plant number 6 is downstream from 4 and 5, and 7 


is downstream from the 6th plant. Thus applying (5.3.14) we have: 
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Ph it) + Ag(t)q,(t) _ Byqq(t)[x, (t) t X(t) + X(t) ] 


4 
+ Bagg (t)Qg(t) + Cyag(t) = 0 (5.3.15) 
Pelt) + Aglt)ag(t) = Byag(t)xg(t) + x(t) 
+ Begg(t)Qe(t) + Ceae(t) = 0 Grey 


Pho(t) + Ay(t)an(t) ~ Byay(t)Exg(t)] + Cpap(t) 


+ Bogo(t)Q,(t) = 0 (5. 3.174 


5.3.2 A Minimum Norm Formulation 


The cost functional to be minimized is: 
ge 


f 
a a [a + eP, (t) + yP, 2(t)]dt (5.3.18) 
i J g 


subject to satisfying (5.3.2), (5.3.3), (5.3.4), (5.3.8), (5.3.11), (5.3.12) 
(5.3.15), (5.3.16) and (5.3.17). Including al] these constraints except 
for (5.3.3), (5.3.8), (5.3.11) and (5.3.12), we have to minimize 
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a ie era é iat 2 ay yt ta). (6. e: a 


“tepaxe ieee samme fi patbufont 
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[x,(t) + x.(t)] - Bon7(t)q,(t)x_(t) jdt (ons. 9} 


boyecteo (5.3.3), (5.3.8), (5.3211) and Gon Seay) 
Consider (5.3.8), the equivalent of which is: 
Sate) 0, (t) Te etl eon (5.3.20) 
and the equivalents of (5.3.11) and (5.3.12) are 


x(t) = vj (tst,) test (5.3.21) 
ey (tate te Ov (nave Utes) 


+ Q;°(t-t5) are (6533228 


The extra functional to be added to (5.3.19) is (so that (5.3.20) and 


(5:53:21) are taken care of): 


veo 
T T 
aye 4,2 Boy f 
~ 1 ids r 
[voto fo 


[0,2(t-r,) + 20;(1; 5140; (t-r;) dt 
(hx3o23) 


Here terms explicitly independent of q.(t)s Q. (t) and x, (t) were neglected. 


The functional given by (5.3.23) can further be reduced to 
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6 Tent. 
2 7 
rs (t)x, Gains = f r, (str. ) 
1=] . 
[0,“(s) + 2v.(t. st. )Q,(s) Ids (ope a5) 
Let 
p.(tst.) a 2y.(tz st, )r.(ttr.) 0 < t < eae 
aS 0 tenes < t s Te 
i=l, 56 (5.3.25) 
= 0 Te-t, <ab. < Te 
lie 8 36 (S320) 


then the augmented cost functional given by Jy of (5.3.19) plus agro 
of (5.3.24) is given by: 


Vs 
Jol.) “J [(e - a{t)(1 = Byo))P, (t) 
0 g 
8 


+ Dtng(t) - a(t)(1 - Byy)]P, (t) 


j 


6 
+ Litn,(t)a,(t) + mj(t) + p,(tsry)J9,(t) 


+ 
It 100 


[n.(t)A,(t) + m,(t) Ja, (t) 


ee 
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Dh (t)0, (t) +P, 2 Ct) 
1=] g 
on 
Ot) dy iy P5 (t)B, P(t) 
8 8 Ban. (t) 
2 
#2, Cert) af ttt — 0, 7(0) 


Define the control vector by: 


u(t) = col.[P(t) Wy (t), Wo(t)s Wy(t)s...sWo(t)] (5.3.28) 


with 
Oe) COUP. Cb) wet tyes On 
P( ) co cP, | ) he 54 
W.(t) = col.[Q.(t).q,(t)] aml Zecea a: (5.3629) 
Wa(t) = col. [Qq(t) saq(t) x, (t) x5 (t) »x3(t)] (5.3.30) 
We(t) = col.[Q¢(t) ag(t) x, (t) x, (t) J (5.3, 31.) 
Wo(t) = col.[Q,(t) s47(t) x6 (t)J (Ge 


and define vector L(t) by 


bit) = col. (L(t) ky, (t)e--vky (0) (5.3.33) 
where 


L(t) = col.[{n, (t) - A(t)(1 - Byp)}seees 


{no(t) - r»(t)(1 - Bag) t.{8 - r»(t)(1 - Bg) 3] 
(53. 34.) 
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ner | 
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p ner 


Ly. (t) = col. Lh, (t).mj(t) + n,(t)A,(t) + py(ts7,)] 
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35} 


36) 


Tn=. 1e25o55 ‘5. 
by {t) = col. [thy (t) sm, (t) + ng (t)A,(t) 

‘s py(tst,) 305040] ‘ie 
Ly (t) = col. [ihe (t) sme (t) + ng(t)Ag(t) 

z, Pe(t.t¢) 0,0] (5. ; 
byt) = col. [th (t) sm. (t) * n5(t)A,(t) 0] tO. oe 
by $t) = col. [tig(t) smo (t) + no(t)Ag(t)] AC 

Let the square symmetric matrix B(t) be given by: 
B(t) = diag[B,, (t eBngp) sso 42) (Se 
with BA(t )= (b; 5(t ) 9x9 (508 
by (t) = r(t)B,, ds fend §9one liwtege eee 
py Pate. ae spe ee (5 
B h(t) 
By, () = diag[-(—5 + 6.(t.t;)),Cin;(t)] 
es aor (5 
B, (t) = diagl-( Fo + eg(tstg)) By  (t)] (5.3.44) 
4 G4 
By nad ty Bn, (t) B 
AC a4 eee 
Ue “ihe 0 
By (t) =| Bae 
CsXy ¥ 7 0 rat) 
a a 0 0 ¢ 
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Bu tey 4 axagp eee 
We a fag et *g(tt6)) By , (t)] 
- (5.3046) 
Bele (tc) Bate ie) 
ery ta, Pees Bells 
Badia) 
B i wa bee r,(t) 0 
VX, 2 4 
Beng (t) 
ee 0 ro (t) 
(5/3. 400 
Bona (t) 
By At) = diag[- mercer By (ec (5.3.48) 
GsX7 
a Bons (t) 
Con, (t) ue : 7 | 
B (ie BLA 510249 
W»X5 eis : : re(t) 
__ Befg(t) 
By,{t) = diag[- a> Cong (t)] (5.3.50) 


Using these definitions, the cost functional given by (5.3.27) 


reduces to: 


ig 
J(u) sl (Llu(t) + ul (t)B(t)u(t) dt (5.3.51) 
O 


Let 
j Tt) = Ptr ee (t) (5.3.52) 


then J in (5.3.51) reduces to 
iu 


f 
yu) =f Ceut) + SEs Ta(tyue) + Gh 
e) 
- OUD Bt) Wt) at (5.3.53) 


The last term in the integrand of (5.3.53) does not depend explicitly on 


u(t), so that it is only necessary to consider 
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From (5.3.40) we have 


with 


then 


Let 


and we have 


j 


j 


B"!(t) = diag[ By (t) By 1 (t) 4+ By “1(4)] 
.52) is rewritten component-wise as: 
W(t) = COT. [Vy (t) aly (toe ++ ody, (19] 
Bye cole LV) “(bya eet); Ven a): 
—- hy ne - 
Wy €t) = col.[V Pies Ve T=) 1k243,0.50 
en CO et eat everest) s\snmc te van mts) a 
ana oy eas Waa 
V,, (t) = SN Bray. ESV Cv c 
Wy CE) = col. Ly (Day (tay, Ct) ty, (8) 
V = ea ea en) 1G 
Wy (t) col.[V  (t) Woo! )oMy )] 

eh Tove 
Wit) = L, 8, (t) 

i = 
a) - Ly, (t) By. (t) 

= 
8, (ty) et (C;5(t) gq 
Tee = diag[- GT ay one 


e.(t.t.) 


ea: Mrs 
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(5.3.54) 


(553455) 


(5y3e56) 


(S.3.07) 


(532580 


(5.3.64) 


(5.3.65) 
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] = ] -] 
eM, (t) diag[- hi aaa 8 - (t)] 
Oe ees Oy(tsty) 4 
(553 


Det) (be en (5.3. 


GX, 4. 


is a symmetric matrix whose relevant elements are given by: 


ae RE he 
Wy (t) oe [Cyn,(t) = ml 1 me (5 3 
1] Ve 1 
Ban, (t) 
12 1] 
B,n,(t) 
My (e) = ero (th, My (t) (5 3 
Ts 1] 
B,n,(t) 
=, Gnd 
My (t) = ete Wy (t (Gus 
14 ] 
Next from (5.3.46) 
Bert) = diagl- cee ee 
We Taal glietb) » Wax 
a Pepe nF O¢(tst¢) (5 3 
with 
B, (EN(t) 2 toy h(t) (5.3 
6 ij 


being a symmetric matrix whose relevant elements are given by: 
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6 
i, (te) e= rE W (t) (5500/5) 
12 1] 
Bene (iG) 
Pt) = 8 by (t) (5.3.76) 
63 are(t) “6 
And from (5.3.48) 
Bi, Tht) sodiagl ———4 Bl (t)] (58 a) 
vy (0) Saale Brey» 8 4 
aX, 
with 
B, e(t}.= (aaa pl) (5.3.78) 
aX7 ij 
where rg(t) 
by (t) = 7 (5.3.79) 
‘1 Bron, (t) 
Con, (t)re(t) - Z 
Bon+/,2 
by (t) Dee _ (5.3.80) 
2 By nz (t) 
Cyn, (t)rg(t) - gL — 
Also from (5.3.50) 
mg Z ] 
B t) = diag[- : J (553,.8))) 
By ( Baig(t) * Cang(t) 


Now the components of the vector V(t) are obtained as follows: 


Using (5.3.34) and (5.3.64) in (5.3.62) one obtains: 


8 
V (t) == y [n. (t) a A(t) [1 > Bs oll, j(t) 


oer (aX .6.2) 
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At) az A(t) @ _ Big) IC; (t) 


Bk - 4(t)(1 - B «Os 
[e - a(t)(1 - Bota (t) (5.3.83) 
Substituting (5.3.35) and (5.3.65) in (5.3.62) we have: 


th. (t) 


Wy €t) = col.[- aR, ay 
6,(t.t.) 


m,(t) + n,(t)A,(t) + p;(tsx,) 


F 
C.n,(t) 
ihvetoa WANs as) (5.3.84) 
Using (5.3.36), (5.3.66) and (5.3.67) in (5.3.62) one obtain 
4(t) 
(t) = - ————_——_ (5.3.85) 


+ 6y(tst,) 


Vy (t) = Emg(t) + ng(t)Ag(t) + Pqltscq) Ibu (t) 


42 rs 
OM, (t) fseqice) 
12 
Ho Ge te tt) (5.3.87) 
N43. Wao ee 
ae) 
Pee ae 43 : (5.3.88) 
Waa Wao Ma 
by (t) 
(t)e= View 0b) rei (5.3.89) 
Wes Wao Wy t 
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‘ sa ane 
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And from (5.3.37) and (5.3.78) in (5.3.62): 


hs 0 tare er a (5.3.90) 
a: O¢(tst¢) 


Vy (t) = Emg(t) + ng(t)Ag(t) + pg(tseg)]b, — (t) 


62 
64; 
by (t) (5.3.91) 
(t) ) le 
“i Rt. (5.3.92) 
We3 N62 We t 
7 
by, _() 
2 13 
Vii ‘=e Vy (t) Seco (5.3.93) 
64 Bee 
1 


Also substituting (5.3.38) and (5.3.77) into (5.3.62) one obtains: 


vot) = drab (5.3.94) 
Wo ~ «BAL (t) at 
2 
tA, (t t 
Vi (t) = uc aaaiae ie Ls aL (5.3.95) 
h2 a] 
Con, (t)r¢(t) aan 
Vwi GE) teeV we (53%) 
Wo3 Woo 2r_(t) Pron 


Finally from (5.3.39) and (5.3.81) in (5.3.62) the following is obtained: 


agit) g(t) + ng(t)AgCt 
Sit ) = col.[- Baiig(t) ° Sa (5.3.97) 


Thus the components of the vector V(t) are determined by the above equations. 


The problem formulated thus far is that of minimizing (5.3.54) subject 


to (5.3.3). Define the 8x1 column vector: 


‘ae 
a Ni 


14] 


B Foon} sucesbad (5.3.98) 
and the (8x21) matrix k! as: 


ie a0 
T T 
K! =19 
ee) ee (5.3.99) 
T 
Ke 
with 
kK, = [0,1] | i = 1,2,3,5,8 (5.3.100) 
d 
K, = [0,1,0,0,0] | (5.3.101) 
: 
ke = [0,1,0,0] (5.3.102) 
Kj = [0,1,0] (5.3.103) 
Then (5.3.3) reduces to 
ae 
af Ku(s)ds (5.3.104) 
O 


The control vector u(t) is considered an element of the Hilbert space 
Pe One| of the 21-vector valued square integrable functions defined on 


[0,7] endowed with the Us product definition: 


a 
<v(t),u(t)> = [Vi (e)e(t)u(e)dt (5.3.105) 
0 
for every V(t) and u(t) in ia Cortes provided that B(t) is positive 


definite. 


The given vector b is considered an element of the Real space R 


with the Euclidean inner product definition: 


>= 


2X,Yo = xy (5.3. 106) 


| 


for every X and Y in R. 


is 
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Equation (5.3.104) defines a bounded linear transformation 


L: ee organ. This can be expressed as: 


b =stu(t)] (5.3.107) 


and the cost functional given by (5.3.54) reduces to: 
V 
Jfu(t)] = |fu(t) + Shed) |? (5.3.108) 
Finally it is necessary only to minimize: 


J[u(t)] = |fu(t) + at) (5.3.109) 
subject to 
by= al [u(t)] for a given b in R. 


9.3.3 The Optimal Solution: 


The optimal solution to the problem formulated in the previous 


subsection, using the results of Chapter 2 is: 
V V 
ee) ibs Ty Sith (5.3.110) 
—t = 2 2 
where i is obtained as follows: 


ue the adjoint of T, is obtained using the identity: 


<&,Tu> 8 eer *e yu> 
a: pier 
Lo plot el (San) 
Let 
& = col.[éy.--+ sg] (5.3.112) 
* 
ae SUG RESP Uy (5:31 41) 
= col. eae 5.3.114) 
Ty col Lt Pg” ( 
Ty Fs col.[ty ty. ] i = 1,2,3,5,8 (5 Sa115) 
1 Y ”) 
Fm col lt, stata: eu. Cae (5.3.116) 
mA Way” Wao’ Wq3° Mag” “qs 


in by sastmtntm ot Kino 
ard at bade.futmrot maton and 


tee = GOL. Leo Ve, Ebaetcts | 
We Wer “62° Wes” Wea 
Tere cola Pt. ste ore 
Wy Wa W790? Wo3 
. = col.[06,] eee Sk 


o, 

dy = col.[0,£,,0,0,0] 

be = col.[0,€,,0,0] 

9 = col.[0,€,,0,0] 
Then in R® one has 


7 
a er 
<E,tu> 5 if aL o, W, dt 
p8 A i=] 


: 21 
and in Lo plot] 


Meek pi i seals 8 
Te > 44 20, (t) 


Substituting the inverse matrices in (5.3.127) this reduces to: 


ES 
j ; 
ty. = col. LO am cep i = 1,2,3,9,8 
Tj = col.[0,é,b, (t),é,b, (thé by (t), 
“Ma 454 410 443 
Ebi be 
4 Wa 
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Ty = col. [0,é¢b, (t),egb,y = (t) seedy = (t)] =~ (5.3. 130) 
°1 ae ig 


Ty = Pye eh7Pw (t),e5b,  (t)] (5.31881) 


: 7 ie 
This completely specifies TE given DV Gos blo). 
The operator J is evaluated from: 
Je] = T[T"e] (S3al32) 


Using (5.3.104) and (5.3.113) this reduces to: 
t 


‘ f 
] 1 
Ojo = com. [ts i Tins tey tthté { tn. (ty att: 
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Te 
1 
es | Conte) 7 dtdstcg| yg 
O 


O 
if 
] 
onloo 
(67 f by, (t)dt (89 | eae (5.9 alos) 
0 


or JUTE] = AE (5.3.134) 
with Te T, Te 

ie diagl | at, [ carey ttf cay tt: 

us } C jn, (t) oN : 33 


(fo dt (5.3.135) 
(eddts | carey ] 2J- 
O 


Thus the operation ee is given by 


; ee, an Ne 
so seouber' eins (er { 
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-] -] 
WEE, = Weok (5534136) 


This yields the pseudo-inverse operation given by 


oP a = 
We =7 [ole] (5c 301a7) 
ass 
Tel, = 0 (5.3. 138) 
- 
bel = col.[0, = ] 
71 f 1 
ni(t) { ater tt 
O 
j =a 2 355.8 (5.3.139) 
by ea My Eq 
+ A 1] 12 
Tél = CON Os = oe : 
i (t)dt Ju, (t)dt 
) 11 ) 1] 
bY, (t)éE, PW, (t)e, 
13 14 
; 1 (5.3.140 
Le lf 
Jou, (t)dt fou, (t)dt 
0 1] 0 1] 
by (t)ée by (te, 
: an giz 
Tely, = col fo. T ) Te ; 
feu, (t)dt J oy, (t)dt 
0 1] 0 1] 
PW (t)é 
2430.5 (5.3.141) 
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by, (EE, by (ty 
af 1] LZ 
Tély Gol sro? T; F ve ] 
f by, (tat i by, _ (tat 
0 1 0) 1] 
(Seg. 142) 
In the expression for the optimal solution (5.3.110) let: 
V ) 
Big De (= ) (5.3.143) 
This is an (8x1) vector whose components are found to be: 
jf m.(t) + ng(t)A.(t) + ps (t,7,) 
eee i i j i Lie 
nN j Zens (et) 
Wal 
fe) 
1 = 1,270 45 (5.3.144) 
of 
be ] 
ng 7 bg + xf Umg(t) + mg(t)Ag(t) + Pg(tsey)] 
0 
by (t)dt (5.3, 145) 
a 


by (t)dt (5.3.146) 
an 
. 
f [my(t) + no(t)A,(t)] 
ny = bs + : if WS ae (RS.147) 
enters gl 
+ 
f [mg(t) + m(t)A, (4) J 
] 8 8 8 
ne be i a Can ; dt (533.048) 
) 


Thus replacing & components in (5.3.138) through (5.3.143) by components 


of n as given in (5.3.144) through (5.3.148) one obtains 


ay 


a ae et tasnocd ert satsev (1x8) 
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(5.3.149) 
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otal Shan 
(t)ny by (t)n, 
413 414 
uf i 
fev (t)dt [> (t)dt 
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ih Ft Me 
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oF 
if (t)dt 
opera 
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Now the optimal solution given by (5.3.110) is obtained component-wise 


as; 


(t)= - = — (5.3.154) 


= ] _- 
ae Bf, (t) #28, (b a; 1 Sel eee” toon e 5) 


Ty ae CREASES 


met) ten. CEPAML YY Pope Sy] 
q-_{t) ee ZC, (t) 


: 
fm, (t) + n,(t)A,(t) + p,(t.x,) 
ef "eT PCE pate a 
i ) 
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“(0 n,(t) t 
itis wis sO al Des hor) 
by, 
a, ,(t) = > HL E(t) + ng (tyAg(t) + py(tsrg)1 
i 
£ 
by (t) [bg + 2 | by (t)[mg(t) + ng(t)A,(t) + pg(tsr4)Idt 
al oe ad 
+ uF 
i b,  (t)dt 
W 
0 1 


(5632) 584 


(wat .£.8) ae 


6 


ng(t)Ag(t) a peltstel 


149 


(t)[mg(t) + ne(t)A(t) + pe(tsre)]dt 
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Plows) | 
ab = Bray he Aw (5.3.166) 


ee 
el = amaae lia (5.3. 167) 


0.3.4 The Modified Optimal Solution 


Here, the pseudo-control variables x(t) and q(t) will be eliminated 


together with the associated multipliers. Rewrite (5.3.157) as: 


2cjny (40, (t) + m,(t) + nj(t)Aj(t) + py(tory) =e 


LP ERAN ree G53 (533.a0g) 


(5.3.158) as: 
Barer ef) + mg(t) + ng(t)Ag(t) + py(tseg) = eg 
44) (5.3. 169) 
(5.3.159) as: 
Bey 4, {t + mg(t) + ne(t)Ag(t) + peltst~) = eg 
a (5.3.170) 
(50.160) as: 
Bey %,(*) + my(t) + no(t)A,(t) = e7 (5.3.71) 
i 
and (5.3.161) as: 
2Cgna(t)Q, (t) + mg(t) + ng(t)Ag(t) = eg (5.3.172) 


8 
where e; are the constants in the corresponding 9; 4%) expressions given in 


(5.3.157) through (5.3.161). Differentiating (5.3.168) through (5.3.172) 


one obtains: 


< 


is bas! vey a 
Pee ca 


15] 


ge C2Cnj(t)Q, (t) + nj (tA, (t)] + h(t) +B, (tyx,) = 0 


elec ds (George) 

cE [ By, TT 0, (t) o ng(t)A,(t) ] fe th, (t) Te py(tst,) =n() 
11 

(5.3.174) 

a [ Beye, a ng(t)Ag(t)] + the (t) + pe (tt) = 0 

"1 (5.3.175) 


& | Br rere, (¢) + no(t)A,(t)] + to(t) = 0  (5.3.176) 
7 
1] 


ae [2Cgng(t)Q, (t) PME ee ane (ea. 0 ce wa 593.177) 


Substituting for p.(t.t,) from (5.3.25) and (5.3.26) in (5.3.173) through 
(5.3.175) one obtains: 

d e 

aee el yny (t)Q_ (t) + n.(t)A.(t)] + ih. (t) s 2y(t. tT) 


j 
j 


6.(t,t,) =e Pea vWe2ess00 ooo Ivo) 


e,(tsty) = 0 (5.32179) 
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Substituting (5.3.155) for h(t) in (5.3.178) through (5.3.180) the 


following is obtained: 


Sel2Cjn,(t)Q, (t) + ng (t)A,(t)] + B.A, (t)Q, (t) 


1 ‘ 
a | bs 


t 20;(t.t;)[v5 (t,t) + ee = 70 


t= Ls 25050 (5.3.181) 
Sat Ber ee, + ng(t)Ag(t)] + Byny(t)Q, (t) 
4 


11 


+ 20 (ts %)L¥q(tq>t,) + 4, = 0 (5. S3he2) 


ci stay 4 (t) + ng(t)Ag(t)] + Bene (t)Q. (t) 


Me *6 =6 
+ 20¢(tste) Lug (test) + Q(t) = 0 (5.32 163) 
Consider the constraints (5.3.11) and (5.3.12) 
x(t) = v.(tst,) ts a, (S2o.ice 
x; (t) = vs (v4 5t;) + Q.(t-t.), tas t< Ty (5.3.185) 


Rewriting (5.3.185) for a time advance T, one has: 


x; (ttr 5) = vj (tg.74) > Q, (t) Ons tae Te-t. (5.3.186) 


The equations (5.3.162) through (5.3.167) can be rewritten as: 


O28¢ 27 (5.3.187) 


f *j 
where for i = 1,2,3 we have j = 4, for i = 4,5 we have j = 6 and for i = 


|e. 


FS3 


And by definition (5.3.26) 


6. (t,t) = r.(ttr.) OCT Teas (5.3.188) 


Thus (5.3.181) through (5.3.183) reduce to: 


ael2Cjn (#0, (+) + nj (tA, (t)] + B)n, (4) | 


+ Bin, (t+r, )Q (ttre) f=0 0 


0, Sites Terry (5.3.190) 


). (ttt.) = 0 (5.3.191) 
+ Pay ig ae te) : 
0 s c < f ‘6 
Note that vor t e(Te-t Te]; equations (5.3.189) through (5.3.191) 
hold true with exception that the last term in the left-hand side (time 
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lead) disappears. Note that (5.3.189) now depends only on n.(t) and 
Q.(t), so that this equation isthe modified optimal equation for this 
type of plants. However, (5.3.150) and (5.3.191) contain the r.(t) 
functions implicitly. This will be considered as follows: 


Consider (5.3.68)rewritten as: 9 


Mi 4% 3 
ty, (eT! = Cegnglt) = gh — al 
1 1s 
(5.3.192) 


substituting (5.3.162) through (5.3.164) in (5.3.192) one obtains: 


S| Byn,(t) 3 
Py Cea = [C,n4(t) = 2, Te 2, cata (5235193) 
i 4 
Also (5.3.74) is rewritten as: 
-] Beng (t) 5 
[by Ce) = [Cen,(t) - 2, (ey hip Xe OM (5.3.194) 
1] 6 


Thus (5.3.190) and (5.3.191) are given by: 


ge eCgng(t)Q, (t) + ng(t)A,(t) = Ba q(t) ) Ne (t)] 


i=] 


+ Byhg(t)Q, Ct) + Bemg(trrg)Q, (tery) a0 
O-< t< Tp “ty (5x 3k 95} 
5 
aeeeCgng(t)Q, (t) + ne(t)Ag(t) - Beng(t) 2 Xe (t)] 


+ Bofie(t)Q, (t) + Byng(ttr,)Q, (ttt,) = 0 
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Note that by (5.3.11) and (5.3.12) we have 
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“Syd les sty tama) Ot re (5.3.198) 


so that (5.3.195) and (5.3.196) are in terms of n.(t) and Q, (t). 
Finally substituting (5.3.156) in (5.3.176) and (5.3.177) one obtains 


ael2yny (40, - Bono (t)x¢(t) + no(t)A,(t)] 


+ B ho (t)Q, (t) = 0 C5 5Sersoy 


d : yw 

geettgngtt)0, Mt) + Ng(t)Ag(t) ] + Banged, ke) =\0 
(5.3.200) 

The modified optimal equations (5.3.189), (5.3.195), (5.3.196), 


(5.3.199) and (5.3.200) can be rewritten in the general form 


Geen; (t)O,_(t) + n.(t)A.(t)] + Bins (t)Q, (t) i g(t) ae 


] i 
(S.SreU'l) 
Here the g.'s are given by: 
g(t) = Sac eile OF Ghett <ilgors 
a ae (5.3.202) 
g,(t) ae 6 Wie <.t < Te 
i = 1,233 (5.3.203) 


ga(t) = Beng(ttt.)0, (tttg) Os t<Te- tg  (5-3.204) 
g(t) = 0 Tete < ts T,  (5.3.205) 
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5 
g¢(t) : B7n7(t+r,)Q, (ttt,) a co[Beng(t) L x eeierad a 


7 jace eee 
Orectns eae (5.3.208) 
d 
g¢(t) ct a dee Bengt) 2 Xe (tory) 
1=4 °j 
Tete cotsedpred.. (5935209) 
go(t) = - E18 n,(t)x¢(t)] (5.3.210) 
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CHAPTER VI 


OPTIMAL HYDRO-THERMAL POWER FLOW 


6.1 Background 


In this chapter extensions of the problems considered in the 
previous chapters are made. Here, the active power balance equation 
used in the previous formulations is replaced by the load flow 
equations. This is as far as the electric network variables are 
concerned. One more consideration is the inclusion of a reliability 
objective in the cost functional. Moreover, realistic inequality 
constraints imposed on the electric variables are considered. In the 
final section of this chapter a practical form of the reservoir 
(trapezoidal) is considered. Furthermore the effect of efficiency 
variations with the active power generation at the hydro-plants is 
included. The problem of implementing the optimum generation schedules 
is illustrated by way of an example. 

The extension of the existing economy dispatch solutions to include 
the exact model of the transmission network is due to Carpentier [58]. 
The resulting optimization problem was shown to be one of nonlinear 
programming. Necessary conditions for optimality were derived using the 
nonlinear programming techniques. In their paper [59], J. Peschon and 
his associates presented the general problem considered by Carpentier 
for an all-thermal system. Another important contribution is that of 


H.W. Dommel and W.F. Tinney [60]. Here optimal power flow solutions are 
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obtained for an all-thermal system. The method is based on Tee flow 
solution by Newton's Method [61], and a gradient adjustment algorithm 
for obtaining the minimum is employed. 

A unified approach to load-flow, minimum loss, and economic 
dispatching problems was presented by A.M. Sasson [62]. Here investigation 
of the application of various nonlinear programming methods to the problem 
was considered. In [63], A.H. El-Abiad and F.J. Jaimes presented a 
variational method to solve the optimal load flow problem. It is noted 
that these two works were also concerned with all-thermal systems. 

The problem of power systems reliability motivated the work of R.L. 
Sullivan and 0.1. Elgerd [64]. An effort to define a reliability 
objective in terms of the system's reactive power generations was made. 
The basic idea of their work was to optimally distribute the reactive 
power generation between the system generators. 

The work by C.M. Shen and M.A. Laughton [65] was of the same nature 
as those previously mentioned. The main contribution here was exploring 
the problem of existence and uniqueness of the optimal solution using 
nonlinear programming techniques. The problem of a hydro-thermal system 
with negligible head variations was solved in [66]. Here a discrete 
formulation was adopted and the problem is solved using the nonlinear 
programming techniques. The discretization process makes the problem a 
one of a large dimension. However a method of splitting the problem into 
ones of smaller dimension was proposed. 

A dual linear programming formulation was given in [67]. Here a 
fast solution can be obtained for the problem of an all-thermal system 


under inequality constraints. This was a contribution to the online 
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dispatching problem. The need for including objectives other than 
economy was given in [68]. Here a minimum emission dispatch problem 
was considered. A variational technique is employed to obtain the 
solution. 

A generalized reduced gradient technique is used for obtaining 
optimal-power flow solutions in [69]. This represents the best method 


to date for solving problems of very high dimension. 


6.2 Statement of the Problem 

A hydro-thermal power system is considered. The system is assumed 
to have N generating plants (generator buses). There are Ny hydro- 
plants and (No - Ny) thermal plants. The system's electric network is 
represented by N buses (or nodes) and these are connected by branches 
or lines having conductance gi and admittance pis, Connected between 
bus 1 and neutral is a branch having conductance g'o and admittance ae 
This is required for the equivalent m representation of transmission 
lines. 

At a given bus i, the phasor source current into the bus is given by 


Kirchhoff's Current Law as: 


E(e)p-e'y (E. (ty - Et 12 (6.2.1) 


where 
yid = gi5 + jp' (6.2.2) 


The net power and reactive volt-ampere at the ith bus is 


5.(t) = Ep (rjry(e) (6)2.3) 


S.(t) = P.(t) - JQ; (t) (6.2.4) 
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* 
where E, (t) represents the conjugate phasor voltage. Let: 


E,(t) = E,(t)e!% (6.2.5) 


E.(t) = E(t) + Eg (t) (6n2.6) 


Then (6.2.3) reduces to: 


Es N 2 SS aes oi eee 
Sl Ee Fy A A G a SAY) Pee E. ee E(t)Y (6.2.7) 
= j=0 
j#i 
substituting (6.2.2) and (6.2.6) in (6.2.7) and separating real and 


imaginery parts of S. to obtain P. and Q. as given by (6.2.4) one obtains: 


’ ‘. e 
Po(t)'= Es>(t)e.»- Egi(t) 2 TE Ce -E Alas 
N e 
a aly 2 (je See GG) (6.2.8) 
een j 
iri ‘ 
eetee iJ ig 
-Q.(t) = E,"(t)B; + Rts ed E kh ] 
j#i 
N . 48 
gob dbl adlE gait )B ye org bag t)G 0] (6.2.9) 
es aif 
where 

G, 3h gl (6.2.10) 

g=0 

j#i 

pes 

B. oF) pid (6.25140 

=o 

Vel 


Note that if the E.'s are the phasor voltage to neutral then it is obvious 
‘ 


eee Cee 0. Hence (6.2.8) and (6.2.9) are rewritten as: 
q 
0 0 


| . Ar) 
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; N - 4g 
Pat) = E,°(t)G; - Ey (t)o} [Eg (t)a' -£.. (t)B'] 
Laer eer tgs 2) J 
j#i 
N i,t, aa 
aE g Wb) IEE wr ae (ova 
th pdelt end j 
Ira 
ee el Ae 
(6.2:12) 
N ri ‘ 
0, (t) = £;°(t)B; +E, (t) fF Ley (told - & (t)Btyy 
In el J 
bg 
N * a 
= toy CEG CGE: eyceg (ech aa 
eat Pm) j 
jfi 
1S goats aly 
(65 21g) 
Note that 

E. = Ey “(t) + * 2(t) i= 1,...5N (6.2.14) 


The equations (6.2.12) and (6.2.13) are the load flow equations. 
Each bus is characterized by four variables P.(t), Q.(t), a and 
as In a normal load flow study, two of the four variables are 
specified and the others must be found. Depending upon which variables 
are specified, the buses can be divided into three types [60]: 

1. Generator bus with P and E specified, Q and the phase angle 
[tan7! a unknown. 


2. Load bus with P and Q specified, EY and EY being the unknowns. 
3. Slack bus with EY and EY specified, P and Q unknowns. For 


convenience this shall be the node No and E(t) is taken as zero. Since 


Bee | 
Byist 


mans 


go 


bes). on] t 
ee | 


162 


the slack bus is taken as a generator Ng node, this means hat the 
number of the unknowns is reduced by one. We may assume that E is not 
specified at the (N-1) bus. 

In the economy dispatch problem, the active power generation 
P at the generator buses are sought. These generations are obtained such 
that maximum economy is achieved. This requires minimizing the operating 
costs at the thermal plants. Thus the problem is: 


" 


: 

‘a 

Min "oh [operigeP) (t) eyo cue hy iq denmacenetisy 
ible ic Male etal Ya 


while satisfying the load flow equations. Note that now the powers P.'s 

are no longer specified for the generator buses in these equations. 
Another important objective in the power system's operation is 

its reliability. To improve the system reliability during operation, it 

is necessary to ensure that the reactive generations are minimally 

proportional between the system generators [32,64]. To achieve this it 


iS suggested that the schedules obtained must 


a molh Pa K, 49; (t)]dt (6.2.16) 


where the Ki 's are assumed to be known weighting coefficients. 
bas es here, that both the economy requirement and the 


reliability requirement are of equal importance. Thus one is required 


LOZ 
: Tr Ny Ny 
Min le} [s.P. (t) +, Ps. t)] + lay ” 
gC) Ng (t) 0 i=N +1 1 j ae 
Q. (t)K; Q(t )}dt (622.1%) 
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while satisfying the load flow equations (6.2.12) and (6.2.13). 
There are several inequality constraints that must be satisfied 


at the optimum for a valid solution. These are given by: 


M 
2 2 2 ; 
P. (t) + Q, (t) < S, T= Ty... Ng (6.2.18) 
Wee Ont) a OS 7 =u N (6,2. 19) 
eee ae eo og a: 
peep. (tee! BS} N (6.2.20) 
j = 7 a i Ses G ° ° 


Furthermore, since E at a generator bus is specified then one 
requires 


Ly. 69 Sano ss (6. 20c 1) 


where E.(t) is assumed to be known. 
The hydro-plants active power generation is assumed to vary with the 
rate of water discharge as: 


Pp (t) + Ay(t)ay(t) + By as (EQ, () + C,4;“(t) = 0 


i= 1,...5N, (6.2.22) 
with 
Pa lhe. 5% Seve 
q(t) yt) iad Np (692428) 
Moreover, the volume of water discharge at any hydro-plant is a prespecified 
constant: 
- 
; = i= N (6.2.24) 
f ajledae = 6, 1 aes A A A 
0 


6.3 A Minimum Norm Formulation 


The problem stated in the previous section can be formulated as a 


minimum norm problem as follows: An augmented cost functional J, can be 


obtained as: 
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2(t)] + m,(t)q,(t) + 1h, 


(64o.72) 


dt 
(6.353) 


(6.3.4) 
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Te Ng 
dL. chaberh p foe 
0, J Pai efOsy fe (8) 
NN 
ae : Q; (t)K; ,Q,(t) ide (6.3.6) 
Pew 
erg M 
dob) =f My (enDP;7(t) + a,7Ct) 5,7 Jat (6.3.71 
I, N 
Gy Ae eee cone 6.3.8 
05 , rie ee i i (6.3.8) 
Te N 
g M 
Jo! ) -{ a Ee) =r lat (0. 3-9) 
O 
Te Ng 
Jog(s) “‘h J. ej (t)LQ," - 9; (t)]¢t (6.3.10) 
ie 
0 
vay 
fg M 
Jost) mi J ef(t) LQ, (t) - Q)"at Gea 
O 


Here an is obtained from (6.2.12) and Ue is obtained from (6.2.13). Also 
1 2 


weeecorresponds to (6.2.21), J. wben(6i2.22) and (6.7.230%" J isi the 
0, Oy 0, 
original cost functional of (6.2.17). The inequality constraints (6.2.18) 


through (6.2.20) are included using the Kuhn-Tucker theorem [59] by 


considering J) (.) through J, (.) so that the following exclusion equations 
10 


6 
must be satisfied at the optimum. 
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M,(t)[P,2(t) + Q,2(t) - S,° ] = 0 (6.3.12) 


2.(t)[P.” - P.(t)] = 0 (5.3598) 


ui(t)[P,(t) - p"] = 0 (6.32145 
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ej(t)[Q," - Q,(t)] = 0 (6.3.15) 
e/(t)[0,(t) - a," = 0 (6.3.16) 
for i = Vs... NG c e[0,T,]. Moreover, Mp, ft) Mg ft) four e n.(t) 


and m; (t) are to be determined such that the corresponding equality 
constraints are satisfied. 
The augmented cost functional Jo(-) of (6.3.1) can also be 


expressed as: 


where 


+ e,(t)Q,™ - e1(t)Q,""Idt (6.3.18) 
Note that Soe is explicitly independent of the control variables. The 
control variables here are P.(t) and Q.(t)Li = V,...5NGI; Eg tt EG res 
hie= 4... | and q.(t)Qy (t)Li = 1,...sN J. And 
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h : Ng 
a Pause: FU ii (M;(t) + y.) 
1-1 i=N +1 
P; (t) at (6.3.19) 


Ng Ng 
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(6.3.20) 
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N oN be 
% F REs (ea Gee aya, (eye TE. (t)]dt 
i=] j=l i 71 i j 
j#i 
(6.3.21) 
Let 
ijt) = -[) (t)Gid +r. (t)BY] 
71 j 
is ls (6.3.22) 
get) = Pe (BY =, (He) 
1 i 
teh = 1, Seahhaiz3 (6.3.23) 
ga (t) = [rp (#6, h eee " re, ()1 
Realy (6.3.24) 
a.;(t) = ry (t)G; + Ng, (t)B: 
7 7 
1 = Ng-Tye eel (6.3.25) 
Then (6.3.21) reduces to 
Te oN : ; 
Jo, ‘f bd) aai(tI eg Ger fake 
17Ng 
; ; (t)ass(t)E, (t) 
e ERM (C)ae-(e 
era oy oF 
j#i 
a) (t)a,.(t)E (t) 
+ bE. At)ase 
iby A Ti faa 
j#i 
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Since EG (t) and EG (t) are specified (this is the slack bus) then the 
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And 
Te Ny 
Jo, = U2, (ng (EDA; (t) + m.(t))q,(t) 
O 
Ni Ni 
+ 2 he (t)Q, (t) + i Cin. (t)q, =o) 
48 
Ni By he (t oe 
: L - (t)]dt (6.3.30) 
1= 
Define the control vector as 
ult) J=<col[P(t)[0(t) ,E(t) WGe) (63531) 
with 
P(t) = col.[P,(t)P.(t)] (6.3.32) 
P(t) =col[ Ps( thy. isi (6.3.33) 
P.(t) = col.[P, tilt )aeeeaPy (t)) (6.3.34) 
Die comlOy(t)...2 0 (te (6.3.35) 
EV G = cole. Eo ) E(t) ] (6.3.36) 
E,(t) = col.[E, (t) dy “AN t dy? (6.3237) 
. 6.3.38 
Pet encoluck 3(t)s Fay e ah ay )] ( ) 
W(t) = 20) 5 Mauer (6.3.39) 
W(t) = col.[q.(t).Q,) (t)] (6.3.40) 
71 


The control is a 2[n+Ny +N Ix column vector function as can be seen by 


inspecting (6.3.31) through (6.3.40). 


Define the auxilliary vector L(t) as: 
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at leche oN 
Let the square matrix B(t) be given by: 
B(t) = diag[B, (t) Bo(t) Bp (t).B,(t)] 
BA (t) = diag[Bp (t) .Bp (tJ 
ee) = diag[M,(t)] deli iterhdi 
Bp (t) = diaglB, (t)J 1 = Nits Ng 
Bp t) = M. (t) +, 
By(t) = (Ki, (t)) 
Ka = Ky ie geet i,j =l, »N 
Ki, (t) M.(t) + Kj; ip ies Ng 
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eee eon ayy (t) re BA al IS er 
Bt Beng, (Hipls ciety tepRle) : eb BY (He One et a eee Pont 

Be) =| yy {*) ayg(t)e- eee ayn (t) | “byy (t)-byp(t) bya (t) : 
0 belt) byg(enby Lg) aygltd wees ay 

ay {t) : pels Pon (t) ! a5, (t) ea) Pie Any (t) 

by (t) bya) byg(t) : ayy (t) ayp(t) ayy t) 

or Be Mt) Me cote i,j # NG 
iJ # NAN, (6.3.66) 


Then (6.3.29) reduces to 
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t 
ne 
a T 
Jop(-) = { LE (edb (EEC) + Le (E)E(E)] (6.3.67) 
6) 


Note that B- (t) is non symmetric. However, one can replace Be (t) 
0 6) 
by the symmetric matrix Be (t) such that 
if Nel 
SBE E(t) EB, (t)E(t) (6.3.68) 
where 
B-(&) £2, (ios 0Ct)) 
“t Cij. | 2(N-1)x2(N-1) (6.3.69) 
with 
Porte bt) (6.3.70) 
Qi, re 
aru ye 
bait) = aC; 5 (t) rs C5, (t))) Tied Views eaten 
re 5.5 # Ny and N4N, (6.3.71) 


Using the above definitions, the augmented cost functional Jaf.) of 


(6.3.17) reduces to: 
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JC.) = | Ch (t)u(t) + uo (t)B(t)u(t) Jdt (6.3.72) 


Oo 


Note that the terms explicitly independent of the control u(t) are 
dropped in (6.3.72) 
Let 


v(t) = UT (eB (ty (6.3473) 


then the cost functional of (6.3.72) becomes: 


; 
fF 
) aff ru(t) «Sy Tacentu(e) + Sqae (6.3.74) 


T V 
Here - V_{t)5(t) =E) was dropped since it is explicitly independent of 
Wee) 
Thus the problem is now reduced to that of minimizing (6.3.74) 


subject to satisfying (6.2.24) which is: 
f 
b. af q,(t)dt ieee ean (o.876} 


Define the NL xl column vector: 


bieeco |e bahon. iby | (6.3.76) 
= ] Ny 


and the ees matrix a as: 


= [0 D5 2g Heh 5 (Goons 


with 0. being the NAXNg matrix whose elements are all zero. 9 and 


Or are NAX(Ng) and N,x(2(N-1) zero matrices respectively. And Ky being 
the N,x2N matrix given by: 
CA a (6.3.78) 
Ky seen sky 
1 Nv 


The Kies being N x2 matrices with zero elements everywhere except at the 
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Ist column and ith row location where the element is one. Thus (6.3).7 5) 


‘reduces to: 


. 
b all K'u(s)ds (6.3.79) 
0 


The control vector u(t) is considered an element of the Hilbert 


2(n+N) +N 


space Lo B g) [0,7 ¢] of the 2(N#N) +N) vector valued square integrable 


functions defined on [0, Ted endowed with the inner product definition: 


t, 
) u(t)> “fate )B(t)u(t)dt (6.3.80) 
O 


2(NRANG +N) 


for every V(t) and u(t) i nbsp [O,T els provided that B(t) is 


positive definite. 
N 
The given vector b is considered an element of the Real space R ! 


with the Euclidean inner product definition. 
<X,¥> = Xly (6.3.80) 
for every X and Y in R 
Equation (6.3.79) defines a bounded linear transformation 
T: LZINNNg NTT Jorn, This can be expressed as: 


= Tutt) (6.3.62) 


and the cost functional given by (6.3.74) reduces to: 


Jolu(t = ||u(t) + aft) 2 (6.3.83) 
Finally it is necessary only to minimize: 
stu(t)] = || (t) + Sey (6.3.84) 
subject to b = T[u(t)] (6.3.85) 
Nh 


for a given b in R 


7 eh 
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6.4 The optimal solution 


The optimal solution to the problem formulated in the previous 


section using the results of Chapter 2 is: 


= T'[b + 14) a3. (6.4.1) 


as 
where Te is obtained as follows: 


t. the adjoint of T, is obtained using the identity: 


<@, lu> N = Teese 
Rh Ente *N) rg 7 i (6.4.2) 
Let 
B= oa ea (6.4.3) 
* 
Te = col.[T, Tote Ty (6.4.4) 


where I,» Ty and qT. are of the same dimension as P,Q and E respectively. 


Jy = col.[T sere ly ] (6.425) 
] Nu 
qy, = col.[Ty, Ty, ] Ten oN (6.4.6) 
q QW 
and 
9; = col.[é. 0 | 7 = T,...5N, (6.447) 
; Nh : 
Then in R-" one obtains: N 
ie 
<é,Tu> = fy ¢, W dt (6.4.8) 
he i=l 
and in ae b*Ng*N) 0,1]: 
if 
* 
TE yu> -f [T,»8,(t)P(t) + TyBy(t)Q(t) + TeBe(t)E(t) 
fe) Nh 


Thus the equality of (6.4.2) reduces to: 


UP 
Ease 0 
feeb 
Ty, 7 4 By C8 i 
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eee Lot 3 0] } 
This completely defines VE as given by 


The operator J is evaluated as 
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this is found to be 


or 


Denes a & 
A - diaal({ Cn, Fs Ctrete 
f@) 


so that the operation tee is given by: 
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This yeilds the pseudo-inverse operation given by: 
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ng(t) nj (t) dt 
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Consider (6.3.74) and let: 
V(t) ‘= col.[V,,(t) V(t) Ve (t) M(t) 
Then using (6.3.72) and (6.3.57) one obtains 


By (t) 


Boe =e 
Wye) = Lyi Ct) 
Furthermore, let 
V,(t) : co) [Yp)(t).Np (E)My(t)] 


using (6.3.42) and (6.3.57) then one obtains: 
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vi (t) = E(Lp, (#851 (EDV (LD (85 (2))1 (6.4.31) 
Thus 

ve (t) =o) (eel Lai 

*, Lp =e ) (6.4.32) 

Vp (t) = Lp (t)Bp (t) (6.4.33) 


From (6.3.50) and (6.3.58) one obtains (6.4.32): 


n.(t) + ai (t) - a. (t) - Ap, tt) 
] 


V, (t) = col. 
oa Sg ee ra LS ee ee 


and from (6.3.51) and (6.3.60) for (6.4.33): 


pj + a(t) - aj(t) ay (t) 
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Mpg tt) = Ch TH Cat mgt 
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i = Mpls cso 


(62435) 
The matrix inverses Bp (t) and Be (t) will be assumed as: 

es 
Do (t) = By (t) (6.4.36) 

-1 
D-(t) = Be (t) (6.4.37) 

with 

= 6.4.56 
Da(t) = (dg, (¢)) (6.4.38) 
D-(t) = (d,  (t)) (6.4.39) 


Note that the dy _(t)'s are functions of M.(t) and the ate Pe ere ee 
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Vg (t) = SP) ally Maen ee cig CE}, (6.4.40) 


and g 


V_(t) = pane ae Ve (t)s...5Vp (ta 


1 N 9) IN 
(6.4.41) 
Note nes and YE. are not present in V_(t). 
Ng Ng 
Then (6.3.52), (6.4.24) and (6.4.38) yield: 
Ng 
Veo Ct)e= > bet) =e. (te tea iCry dd ent) 
j = V,-..4N, (6.4.42) 
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Consider (6.4.29) and let 


V(t) = Ol TW (Ey. sy (td (6.4.45) 


Wy.) = by (t)Byy Ct) oe isa ae if (6.4.46) 
if WaaCt) = col Vi, Melee celta 
W Ng, Na, 


then (6.3.55) and (6.3.65) yield: 


Ste 1 i s 
Pee) = ee i= TseeeN, (6.4.47) 
q. WA 
2h, (t) | 
MW (t) = - Bn, (t) 1 = T,... oN, (6.4.48) 
1 


The optimal solution as given in (6.4.1) is written in component form 


Eeioyss T'[b +T Ey) ee (6.4.49) 
Qt) = 7 Vien Ally - =o (6.4.50) 
E(t) = T HES Oy). - ait (6.4.51) 
W(t) a ie rh, - Suse (6.4.52) 


From (6.4.21) through (6.4.23) one obtains for (6.4.49) through (6.4.52) 


v_(t) 

Pea Sa (6.4.53) 
V(t) 

Q, (t) = - hae (6.4.54) 
V(t 

E(t) = - oat (6.4.55) 
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Moreover W(t) in (6.4.52) is rewritten component-wise as: 


a thee rthya. == vreal (6.4.56) 
=j 


using (6.4.24), (6.4.49) and (6.4.50) this yields: 


[m.(t) + n.(t)A.(t)] 


j ah nett 
TF im, (t) mCP (t) 
ge dt 
O 
2 


Furthermore, using (6.4.33), (6.4.34) and (6.4.35) in (6.4.53) one 


obtains: 


[ay + aj(t) - 24(t) = rp (8) 


184 


Also from (6.4.42) 


dy (t) 
emetic Sout. 
co 5 [ gh ‘age + pie ae. 
J#N, 
2(N-1) 
‘ “5 Coy jen a7 't) + Os y atin CE) 45 (2) 
j#N_+N 
g 
en (6.4.62) 
dy (t) 
eee. t t)]d 
q. ( ) y) La yd ) + 5H, )] a nae 
‘ i#N, 
ie b ft heb (t)]d 
- yu +.B t 
jf [b, j-N.+1 j-Njt1,N, ] E(u) a 
j#NIN, 


= Tesh (6.4.63) 


The last two equations are modified by noting that the symmetric 


matrix Be in (6.3.69) can be written as: 


A-(t) C,(t) 
Ba (thes ‘ (6.4.64) 
ce) A, (t) 
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C(t) = (C.. (t 6.4.66 
= "JE °° (N-1)x(N-1) \ 
-aht) 7a, 3(t) 
ij, ¢) - (6.4.67) 
ob. (t) + belt) 
ig, 6") = (6.4.68) 


Note that A, (t) is symmetric but C, (t) is not. Moreover C, (t) = - Ce(t). 


The inverse matrix Bo! (t) is denoted by 


E 
(t) = pe (t) (6.4.69) 


Then using (6.4.67) and the property C,(t) = - Cats one obtains for 


D-(t): 


F(t) He (t) “ie 
0 
Det) 4. 
—£E 
-He (t) F.(t) 
where 
-1 -1 
Fe(t) = [A-(t) + a (t)C,(t)] (6.4.71) 
He(t) = A(t) Ge (t)E p(t) (6.4.72) 
Thus 
Ye Mt) Z Le (ede) : be (EHEC) (6.4.73) 
ye {t) Le CedHEtt) + be (OEE) (6.4.74) 


so that (6.4.65) and (6.4.66) can be written in the alternative form: 
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“dy . N 
“ty (e) ar = [ iby (ay (ty) 5 as ON (t))f; 5 (t) 
& Z O3J g 
J#NG 
N 
By Ong, + Pian, (en 5(H0] 
J#N 
g (6.4.75) 
“dy m N 
Ey eee 2 (ay 5%) aj yy (ED hg (O) 
: i#N, 
N 
Po N ool + EN (6.4.76) 
JANG 


6.5 Implementing the Optimal Solution 


The method suggested for actually implementing the optimal solution 
is best illustrated by way of an example. The example concerns a practical 
power system as shown in Fig. (6.1). Here bus number one is a hydro- 
generator bus, bus number two is a thermal plant's bus (also slack bus) 
and bus number three is a load bus. The classification of the various 
variables in the sample system is given in Table (6.1). 


The reliability matrix K cost functional (6.2.16) is taken as 


I+K, 2 i 
2 2 
ee (6750) 
K, +k, 
ee 2 
K, = 0.02 (6.5.2) 


] 


This is the same reliability matrix as the one given in [64]. Thus the 
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Hydro Thermal 
plant 


{Fr}, 
2 


Load 
bus 


Figure 6.1 EXAMPLE 3-BUS SYSTEM 


Bus Type E E Es p Q qf a(t)dt 
1@) 


] Hydro Optimum Optimum Optimum Optimum Optimum Optimum Specified 


2 Thermal Specified Specified Specified Optimum Optimum 
(Slack) 


3 Load Optimum Optimum Optimum Specified Specified 


TABLE 6.1 CLASSIFICATION OF VARIABLES IN THE 
EXAMPLE SYSTEM. 


wy ha 


ee 
Ve rs 


matrix By(t) in (6.3.61) is given by 


I+K, K Z 
(M(t) +> (- sh 
By(t) ‘ Ky 1+K (65543) 
(- 5—) (M(t) + ? ) 


M(t) + (ky 
*Q *Q 
Dy(t) = 14K, (6.5.4) 
ae M(t) + (=-) 
2hq AQ 
+k, 142k, 
Ao = M, (tM, (t) + (M, (t) + M,(t)) ay (6.5.5) 


The optimum active and reactive powers at buses 1 and 2 are given by 


(6.4.59) through (6.4.61) as: 


Em (t) + 24(t) = a(t) = ap CE) 


P, (t) = - M(t (625264 
E ] 
[Bp + a(t) - tp(t) ~ ay (01 
Ae =- TT tena] (6e5.47.) 
14/64 as 
ay (2) # = mpl ley(t) = eC) + 2g, (00H (8) + UD 
1 
+ (ex(t) = ep(t) +a, (t)) cA (6.5.8) 
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K 
Q(t) = - Zit (ey (t - ex(t) #24 (t)) xt 
+k 
* (ep(t) = ep(t) + ag Ct) (M(t) + HI 
VG 55.9) 


According to (6.3.66), the matrix Be (t) is given by 


(t) = (6.5.12) 


lee Wi ia 


ae an iat 
a ena tr 
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C(t) = (6.51135 


Aa ; a, (t) 5 


E 
2 
; mee. ack) (6.5.15) 
Also using (6.4.72): 
: by3(t) - b3,(t) 
2 
He (t) ~ t 

b3,(t) - by3(t) ; 

é (6.5.16) 


Using (6.5.14) through (6.5.16) in (6.4.75) and (6.4.76) one obtains: 
Ey (t) 
2 in epee asad t ler aonb) ante) 
ea eae (a5 ay2(t))a3, 


: de. (tow as-( @) 
~ (a53(t) + ago(t)) (2, 4) 


EY 


o Pa mene 
we neg 


oer ia 


19] 


#2 bps (Epret be wey geet bias send nag 
(6.5197) 
Ee alta) 
d ele) ig (t) 
a, 77 FE Flagg (t) + a(t) 
E 
¥ (a54(t) + a35(t) ay, (t) 
DAE Dae bt) 
- (boy (t) + by (t)) (A —, 89) 
(6.5.18) 
Ey, (t) Ree e bet) 
: 1 3 
Fa = = ge — [aya (t) + ago(t)) (2 ! 
if (b,, (t) a by 5(t))az4(t) + (b(t) 43 b35(t)) 
(tt (t) 
13 : mba Fy (6.5.19) 
ED AE) 
d bee CE ane e) 
2 31 13 
“Q, (t) aie Ohta [-(a5,(t) + ay o(t)) ( y) ) 
ase (t)ateasn(t) 
+ (byq(t) + byy(t)) (2, 44) 


is (bo3(t) = b5(t))a,,(t)] (6.5.20) 


The optimal expressions for q, (t) and My, 6% as given by (6.4.57) 


and (6.4.58) can be combined to yield: 
d , 
[2c no (t)Q, (t) + ny (tyA, (t)] # Byons (t)Q)) (t) = 0 
dt W ] W, ] ] W, ] Ws 


Q, (0) = 0 Qy. (T_) = by (655721) 
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Here the optimal values for the unknown variables are obtained in 


terms of the functions Le A TA) sy 0 Re ayant ele 
p(t)» App (ts Ap (EDs Ag (EDs Aga (Ee dg. 


and n,(t). These are to be determined such that the equality constraints 


(t) 


fOre.t2) for 1.= 152.3, (6:2.13)-for i = 1.2.3 and (GTZ. ce) toy ie | 
are satisfied. Note that the Kuhn-Tucker multipliers M.(t), &.(t), 
a(t), e, (t) and e.(t) for 1 = 1,2, are determined in accordance with the 
exclusion equations (6.3.12) through (6.3.16). 

In the case when none of the inequality constraints is violated, 
the Kuhn-Tucker multipliers are zeros. Thus (6.5.7) yields 
y(t) - Bo | 
P. (t) = cay La pe (Oo 56e2) 


(6.5.8) and (6.5.9) reduce to: 


] 

5 oi THAR LK kg 1 ig Krg, (td (6¢5¢28) 
] 

pen? THO LK (t) t+ Bs iota told (6.5.24) 


while (6.5.6) yields 


ny (t) = p(t) (6.5925) 


Note that this result can be obtained if P,(t) in the formulation of the 
problem was not taken as a control. This is the case when there is no 


second order term in P, (t) in the cost functional or the associated 


constraints. 


The number of unknown variables in this example is 16, these are 


divided into two categories: 


i, a 
he 
_ 
rl 
Le ae 
} 


_abtaty (a. ») aut 120785 


SUED hh OM AVON GS: 


‘| ’ dh i " ‘ L, \ a A i 
(esta) Lr s + 4, , pith 
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1. Physical variables: these are P, Gea P. (ie Q; (cts) Q. Cele 
a cae Ea Gar ae (t) and M, (it ealihie Biyeiean rare oe 
Pete tad directly by fic load fide tise (6.2.12) Gand? (6.2e13)) for 
i = 1,2,;3 and the hydro power equation (6.2.22). Note here that the 
number of the physical variables is 9 while the constraining equations 
BGen/.. 

2. Dual variables: these are arate X(t) (Giv=bec co) cand n.(t). 
These are related to the physical variables by the optimal expressions 
(6.5.6) through (6.5.9), (6.5.17) through (6.5.20) and (6.5.21). The 
number of the dual variables is 7 and the constraining equations are 9. 
Thus one has 16 equations in 16 unknowns. 

It is possible in this example to reduce the number of equations 
and unknown variables by algebraic manipulations. However, this will not 
be done for many reasons. One of the reasons is that the reduced equations 
for this example are not easily extended to higher dimensions. Another 
reason is that it is important to keep the variables as they are in order 
to satisfy the inequality constraints and the corresponding exclusion 
equations. The third and most important one is that the load flow equations 
can be solved easily with the present state of the art. Thus developing 
a computer program that is an extension of a load flow program is highly 
desirable. 


The Example 3-bus System's particulars are summarized in Table 6.2. 


The results of the optimum load flow study are shown in Fig. (6.2) through 
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TABLE 6.2 


EXAMPLE 3-NODE SYSTEM'S PARTICULARS 


The Thermal Plant Cost: 


g = 4.0 Seo bela tne: 
The Electric Network: (Admittances in mhos). 
al = o.1ayx10 BMS covegx10oe eg ee 
i?j 
Gg’? = 0.0 Bp’? = 0.05x10°°. i = 1.2.3 


The Electric Variables Specifications: (Voltages line to line). 


Slack Bus: Patty =aceOany EF (t)-= "0.0 
qs % 
Load Node: P(t) = -41.5MW Q3(t) =-14 MVAR 
Load at Node 1: Ph (t) = -50 MW Q (t) = +30 MVAR 
] ] 
Load at Node 2: P. (t) = -50 MW Qreee(t moe 4/5 MVAR 
Ds wy 


The Hydro Plant: 


13.35x10°Ft® 


n = 0.708 V 


1/7.2x10 19 g, = 0.0 
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Figure 6.2 ACTIVE Power Variations. 
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Figure 6.3 
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Variation of the Direct Components of the Voltages. 
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Figure 6.4 
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Variation of Quadrature Components of the Voltages. 
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6.6 Trapezoidal Reservoirs and Variable Efficiency Hydro-plants Consider- 


ations 


In formulating the problems in Chapters 5 and 6 it was assumed that 
theefficiency of each hydro-plant remains constant over the operating 
range. Another assumption is that of vertical sided reservoirs at the 
hydro-plants. In this section these two assumptions are relaxed and the 
modifications to both formulation and optimal solution are shown. Here 
the power system problem stated in Section 6.2 is considered. 

The ith hydro-plant's active power generation is given by: 

h (t)G,(t) = h.(t)q.(t) (6.6. h) 
This is precisely equation (5.2.13) except that here the inverse efficiency 
G. is no longer a constant. 

The effective hydraulic head at the ith hydro-plant is given by 

(5.2.2) which is: 


h(t) = y,(t) - yy (t) (6.6.2) 
The forebay elevation isrelated to the forebay volume of water stored 


S;(t) by the relation: 


2 
= : + ait 6,623 
S,(t) oy Yi (t) By Ya) ) ( ) 
where ay and By are constants for the trapezoidal reservoir given by: 
i j 
y = he tang. 
8 = 2. D 


adit ts entorroeat ibaa Ieotvey Yo edd at mgbiqrves’ 
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The geometry of a trapezoidal reservoir is shown in Fig.(6.5). The 


volume of water stored S.(t) is also given by: 


t 
Sst} = S. (0) [00 ue - { ay (e)dc (on 
(0) 


This is the reservoir's dynamic equation in the case of no hydraulic 
coupling between the plants. Here i, (t) and q;(t) are the rates of water 
inflow and discharge respectively. The volume of water stored variable 


S.(t) can be eliminated by ele os (696;3)) 1n) (6.6.4) torobtain: 


j 
t 

a y4@(t) + By i 0 + faxt \do . S. o) -f Til o)do = 0 
(0) 


The tail-race elevation yz, (t) is given by: 
j 


yet) oy) oF Beige t) (6.6.6) 


This is the same equation as (5.2.4) 
Thus the active power generation equation (6.6.1) becomes: 
Py, (t)G;(t) + ¥y_ a(t) + By ay (t) = aj (thy, (t) = 0 
(6.6.7) 
Here (6.6.2) and (6.6.6) were substituted in (6.6.1) to obtain (6.6.7). 
For all practical purposes the variation of the inverse efficiency 


G, (t) with the active power generation can be represented by 


2 
Goce) F236 Gilt aye 
oe ie 95, 95 9; 
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“nee ito shai vides 


201 


Vy 


Figure 6.6 A Trapezoidal Reservoir. 
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This is the equation of an ellipse and is assumed to hold true over the 


operating range ce < P(t) < pM) of the hydro-plant. This is shown 
j 
mn Fig 6.7) 
-] 
6-11.87 


P( KW) 


MIN. MAX. 


Figure 6.7. Efficiency versus Power Output. 
The performance of the ith hydro-plant is completely specified by 
the relations (6.6.5), (6.6.7) and (6.6.8). Thus the only modification 
to the formulation of the problem in Section (6.2) is replacing (6.2.22) 
by the three relations mentioned above. Note here that two more variables 
per plant are introduced. These are y(t) and G.(t). Accordingly the 


functional Jo (.) in (6.3.5) becomes: 
4 


t, 
()=f'6 Pin Pp, (BIG, (t) + yp (8) a5 (t) 
0 


ae 
- fi, (odaap ot (6.6.9) 
oO 


enn bas @®y X ate se 
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Note that the variable Qy, in (6.3.5) is no longer of ngeun Thus 
jn effect only one extra variable is needed. 

Here the control variables are Ph, ft) G.(t), q.(t) and y(t). 
Thus the terms explicitly independent of these variables can be dropped 


from Jo (.) so that one needs to consider only 
4 
a N 
Jeu) =i [ Ping (t)[P, (t)G;(t) + yp a, (t) 
4 i i=] i io 


2(t) - as(t)y,(t)] + m,(t)La, 6,°(t) 


j g, 1 


t) + 8) ys(t)] - ri (t)q,(t)]dt 


(6.6.10) 
As for the augmented cost functional Jo(-) Of (6.347), the only 
changes are in Jo (.) of (6.3.19) and Sy he of (6.3.30), these will be 


given by: 
(t) JP. (t) 


Pj 


+] 


Te Ng 
J Ohya [aio Ms (een (t) 
fe) Np 


; 
g 2 
Ht Ta § M. (t) 1% y;)P. (t) 
h 
Ny 
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+ ‘te (t) - Q. (t) - Mp AP) + my fEIS¢ J 
Nh 
PsA Slade Le m(t)yg, 1P,°(t) 
Ny Ny ; 
+ J ni (t)P,(t)G,(t) + J mi(t)o, @;°(t) 
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N 

+ y ny) (OTe fea HheN ae (6eG. ul) 
T, Ni 

Jo fe) = [Ung (tlyy, ry (t) 19, (t) 


2 
+ ee My MOLey ¥, (t) + By y(t) Ide (6.6.12) 


The control vector defined in (6.3.31) is modified as follows: 


P(t) = CON ERIE At) Saar rye L aan de (6.6.13) 


= col.[q,(t),y.(t)] (6.6.14) 


= 
— 
ct 
~ 
| 


The last two equations correspond to (6.3.33) and (6.3.40) respectively. 
The control vector is a 2[n + TSN, ANG Ix! column vector function in this 
case. 


The auxilliary vector L(t) of (6.3.32) is modified to: 


Ly (t) = col.[L, (t),.-.sbp (t)] (6.6.15) 
h 1 Ni 
Eyl) = Gol Ltn (t)ayabe) + a Satya Mp, (ED) 
im (t)8g J i= 1,...5N, (6.6.16) 
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Finally the matrix B(t) of (6.3.47) is modified to: 


a3 = “agtB, (6) i = Ty... Np (6.6.18) 
n.(t) 
Im; (t) +m (t)yg J > 
B, (t) = 
me n,(t) 
5 eae 
iad ee, oN, (6.6.19) 
Wee) 
By ny (t) uc: 
j 
By (t) = 
: n,(t) 
ry ay, P(t) 


i= T,...5N, (6.6.20) 


These are the changes as far as the formulation is concerned. 

The modification in deriving the optimal solution as obtained in 
Section (6.4) follows easily. First in obtaining Te of (6.4.4) the 
following modification is needed to (6.4.14) using (6.6.20) 

(t) 


Se Penn) 
71 1 
Gis) > vtaeind eat tee ae 
H t = MS B 10) . ra ° e 


is the determinant of the matrix in (6.6.20). Second, the operator J is 


obtained as: 


f 


Payee eae 


sath 


a Luna) » ae 
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J 3 3 
Le } coll(éja, a tty ot) 1 T,...5N, 
[@) 1 


Thus A of (6.4.18) becomes: 


‘ 
eh a, Ce dt aa, yf xt) dt 
O 1 0 1 
i= 1,-.-.N, (6.6.25) 


The expression for V(t) as given by (6.4.25) will not be changed. 


Only the expressions for its components Vp (t) and Vy (t) will be modified. 
h i 
Thus (6.4.34) is modified to: 


i= 

aS 
ct 

~~ 
iT} 


cOT.LV, (t) Vo (t)seeeaV, Ve (t)] (6.6.26) 
Py 1 N 


with: 
V(t) = {fm (t)o, + 25(t) - 25 (t) - a, (t) Im; (t)o 
1 1 


a. 0.8m; (t)nj(t)Bg Weoy (t)3 (6.6.27) 
1 
- 0.5n, (t) 


Ve Ct) = CONE ho mi(t)yg Imj(t)eg. 


[m.(t)d, + 23(t) - R(t) - Mp MED / ay (t)} 


9. 
1 i ‘ 


(6.6.28) 
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ty, (t) = my(t)ag OM, (t) + yg m,(t)] - (n,2(t)74) 


P. , 1 1 
(6.6.29) 
Also (6.4.46), (6.4.47) and (6.4.48) are modified to: 
Wy, () = col.[V,, (t); Vy rs (6.6.30) 
q. Y; 
Vw (t) = {Ln.(t)y % rj(t)la, *;(t) + 0.5n;(t)t, (t) 
q. io j 
hye aie (626.31) 
ia = ee - r,(t)] . acre 
n.(t)P.(t)}/ta, (t)) (6.6.32) 


j 


The optimal solution as given by (6.4.53) through (6.4.56) is 


unchanged. However, component-wise only (6.4.57) through (6.4.59) change to: 


~ 
£ 
q(t) = [-Wy (t)/2] + {Lb, +f [Wy (t)/2]dt]r, (t)/ 
a 9; @) 9; 
Mf (t) 
ay] ay ey (6.6% 3) 
Li 
tye="E-V,, (t)/21 Abe ee Gere idiin (ty 
Ye | Wy 1 j Ma, 1 
Senet) 
ty (0 J ay, CE) dt]} (6.6.34) 
Bea Ge eea'p (t)/2 (6.6.35) 
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(6.6.36) 


(ii Uae ae (t) and G. (t) are to satisfy 
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CHAPTER VII 
CONCLUDING REMARKS 
7.1 Conclusions 


In this thesis a functional analytic optimization technique is 
applied to problems of economy scheduling of hydro-thermal electric 
power systems. Here, the minimum norm formulation is employed to find 
the optimum generation schedules. This investigation shows how the 
powerful minimum norm formulation can be applied to complex problems 
of high dimension. 

In Chapters 3 and 4, some simplified economy scheduling problems 
are considered. The problems are posed and solved using the minimum 
norm formulation. These problems were investigated earlier using other 
optimization techniques. However, the solution obtained here is guaranteed 
to be the unique optimal solution. Moreover, limitations on the unknown 
functions obtained through this particular formulation facilitates the 
practical implementation of the optimal solution. A further simplification 
is the elimination of the multipliers associated with constraints that 
are linear in the control vector. The solutions obtained here are easily 
shown to agree with the previously obtained solutions using other methods. 
This provides a firm ground for investigating more complex problems. 

The problems posed and solved in Chapters 5 and 6 represent one 


main contribution of this investigation. So far as the author of this 
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thesis knows, there has been little or no work in the following areas: 
First, the time delay of flow between hydro-plants on the same stream 
is included in the formulations of Chapter 5. Also, the tail-race 
elevation effect on the operating hydraulic head is considered here. 
Second, the formulation in Chapter 6 is in terms of the exact model of 
the electric network. The reliability objective and the practical 
limitations on the network variables are also considered. The formulation 
concerns itself with a hydro-thermal system with variable head hydro- 
plants. Third, the effects of efficiency variations were incorporated 
and a trapezoidal reservoir is considered. These formulations in addition 
enjoy the advantages cited before. 

An important aspect of the problems considered in this thesis is 
the computational schemes adopted. Due to the nonlinearity of the 
resulting equations, the method is iterative in nature. Here employing 
the modified contraction mapping principles is very useful. Furthermore, 
the transformation of the differential equations into operator equations 
in Chapter 5 guarantees satisfaction of the boundary conditions at each 


jteration step. 


7.2 Suggestions for Further Research 


The minimum norm formulation employed in this investigation has 
demonstrated the capability of solving complex power system scheduling 
problems. Further research with the same technique would be desirable 
in order to explore the possibility of solving more complex problems. 

For example, it may be possible to solve common-flow problems in the case 


where the time delay of water flow is a function of the rate of water 
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discharge. The justification of the constant time delay fe the 
fact that these delays are of the order of a fraction of the optimization 
interval considered. Also efforts in the field of defining an overall 
reliability functional are highly desirable. 

The optimization method used in this research also holds promise 
for related problem areas of the electric power industry. Similar 
system optimization problems arise in the design of new power transmission 
systems and in determining favourable locations for building new generation 
stations. In addition, the system design and optimization of new trans- 
mission line parameters as independent variables in the optimization 


problems are problems that may be handled using this method. 


1 cm 


aD | rapa i pie i 


/ ni nt aiuiaieligy wm ee 
NWR AY be 
iy had ten! ae 


oy ot 


ie | beyr CT 


\ de Wing bug 


212 


BIBLIOGRAPHY 


IEEE Working Group Report on "Present practices in the Economic 
Operation of Power Systems", IEEE Trans.,PAS-90, No. 4, pp. 1768 
-1775. 


F. Noakes, and A. Arismunandar, "Bibliography on Optimum Operation 


of Power Systems: 1919-1959", AIEE Trans., PAS-81, pp. 864-871. 


I. Hano, Y. Tamura, and S. Narita, "An Application of the Maximum 
Principle to the Most Economical Operation of Power Systems", 


IEEE Trans., PAS-85, No. 5, pp. 486-494. 


C.W. Watchorn, “Inside Hydro-thermal Coordination", IEEE Trans., 


PAS-86, No. 1, pp. 106-117. 


E.B. Dahlin, "Theoretical and Computational Aspects of Optimal 
Principles with Special Application to Power System Operation", 


Doctoral Dissertation, University of Pennsylvania, 1964. 


R.J. Cypser, “Computer Search for Economical Operation of a Hydro- 
thermal Electric System", AIEE Trans., Pt. III, October 1954, 


pp. 1260-1267. 


2 >imonoos ant at 259 atten ‘insegy no srogas si ' 
bi bg Lee: My tier y it 
“ear 9 ae ot 0-204 2 enevT at : “enareye 


=n Hh 


ere MiNi re me 
notserog0 mite fo © wigenpot Fata ribet * bre rig 


| ) | t sc ait yd Sah a 
i) ne Deh gay fy | ‘ea 4 ag ‘ oy 
\ Ms ; ‘ : 7 t vad 7 r it i : re F iy ; 
a A Nie ( r y oy) ee nape J P ; ‘| ; en ay mel 
Ma Md eet ey! 
bil 


oe é ed « ae it “enaee aa Pc: 7) coe ur ane - 


Bite vane 7331 L"notdanies) Femods-owoat oe 2 ou vm wi Wee 


Prey Ree. wah f ‘a ae 
re ae ' . p 


-enag0 to atooaeh “ant sumo bas jaa 
i | “snoisere40 mosey? ‘130 ‘ot boseaa inet tom 
lt eer ef sstnsvivannot Yo 1 


12. 


is. 


14. 


ZS 


J.H. Drake, L.K. Kirchmayer, R.B. Mayall, and H. Wood, "Optimum 
Operation of a Hydro-thermal System", AIEE Trans., Pt. III, August 
1962, pp. 242-250. 


A. Arismunandar, Discussion of Ref. 7. 


W.A. Porter, "Modern Foundations of Systems Engineering", MacMillan, 


New York, 1967. 


L.K. Kirchmayer, "Economic Operation of Power Systems", John Wiley 


and Sons, New York, 1958. 


A.F. Glimm, and L.K. Kirchmayer, "Economic Operation of Variable- 
Head Hydroelectric Plants", AIEE Trans., Pt. III, December 1958, 
pp. 1070-1079. 


D.K. Hughes, "Variational and Optimal Control Problems with Delayed 
Arguments", Journal of Optimization Theory and Applications, Vol. 2, 


No. 1, pp. 1-14, 1968. 


F.M. Kirillova, "Applications of Functional Analysis to the Theory 
of Optimal Processes", J. SIAM Control, Vol. 5, No. 1, pp. 25-50, 
1967. 


D.S. Carter, "A Minimum-Maximum Problem for Differential Expressions", 


Canadian Journal of Mathematics, Vol. 9, pp. 132-140, 1957. Also 


| -alashaa to. otSer89d Stier 


-B2er redneo90 MEL ities ane 320i 


16. 


te. 


19. 


20. 


oie 


214 


correction ibid p. 226. 


W.T. Reid, "Ordinary Linear Differential Operators of Minimum 


Norm", Duke Mathematical Journal, Vol. 29, pp. 591-606, 1962. 


A.V. Balakrishnan, "An Operator Theoretic Formulation of a Class 
of Control Problems and a Steepest Descent Method of Solution", 


J. SIAM Control, Vol. 1, No. 2, pp. 109-127, 1963. 


L. W. Neustadt, "Minimum Effort Control Systems", J. SIAM Control, 
Nol? 1NeNoveds  ppest6-31, 1962. 


G.M. Kranc and P.E. Sarachick, "An Application of Functional 
Analysis to the Optimal Control Problem., ASME Trans., Journal of 


Basic Engineering, pp. 143-150, 1963. 


W.A. Porter, "A New Approach to the General Minimum Energy Problem", 


J.A.C.C, Stanford, California, 1964. 


W.A. Porter and J.P. Williams, "A Note on the Minimum Effort Problem", 


Journal of Mathematical Analysis and Applications, 13, pp. 251-264, 


1966. 


W.A. Porter, and J.P. Williams, "Extensions of the Minimum Effort 
Control Problem", Journal of Mathematical Analysis and Applications, 


13, pp. 535-549, 1966. 


q A ANiyie re 
' hy at 


: mh Pat sate 6  bottsw snose50 seansed2 tne ae 
: eae .TSE-207~ 99 S “om tt. 


tty ee aD %.. YES AS 


tas08 | Yo tsar juts 2 
ene aMeA e -watdon * trai 


Gc. 


(aS 


24. 


Zo, 


26. 


2h. 


23% 


Aa 


215 


W.A. Porter, "On the Optimal Control of Distributive Systems", 
J. SIAM Control, Vol. 4, No. 3, pp. 466-472, 1966. 


W.A. Porter, "On Function Space Pursuit-evasion Games", J. SIAM 


control, Vol. 5, No. 4, pp. 555-574, 1967. 


W.A. Porter, "A Basic Optimization Problem in Linear Systems", 


Mathematical System Theory, Vol. 5, No. 1, 20-44, 1971. 


G.F. Simmons, Introduction to Topology and Modern Analysis, McGraw- 


Hill, New York, 1957. 


B.Z. Vulikh, Introduction to Functional Analysis for Scientists and 


Technologists, Addison-Wesley, Reading, Mass. 1963. 


A.N. Kolmogrov and S.V. Fomin, Elements of the Theory of Functions 
and Functional Analysis, Vol. 1, Metric and Normed Spaces, Graylock, 


Albany, N.Y. 1957. 


D.G. Luenberger, Optimization by Vector Space Methods, Wiley, New 


York, 1969: 


M.J. Steinberg and T.H. Smith, "The Theory of Incremental Rates", 
Part I, Electrical Engineering, March 1934, Part II, Electrical 


Engineering, April 1934. 


4 


303 


St 


32. 


335 


34. 


35. 


36. 


216 


M.J. Steinberg and T.H. Smith, "Economy Loading of Power Plants 


and Electric Systems", John Wiley and Sons, New York, 1943. 


M.J. Steinberg, "Incremental Maintenance Costs of Steam Electric 
Generating Stations", AIEE Trans., PAS, No. 34, pp. 1251-1254, 
Feb. 1958. 


0.1. Elgerd, Electric Energy System Theory, McGraw-Hill, New York, 
Tort. 


D.P. Lijesen and J. Rosing, "Adaptive Forecasting of Hourly Loads 
Based on Load Measurements and Weather Information", IEEE Trans., 


PAS-90, pp. 1757-1767, 1971. 


A.M. Sasson, "Economic Operation of Power Systems", Lecture notes 
presented at The Purdue University Short Course in Power System 


Engineering, Lafayette, Indiana, 1971. 


R.D. Chenoweth, "A Study of B Coefficients for Loss Prediction 
and Economic Dispatch", 71CP98PWR, Presented at the IEEE Winter 


Power Meeting, New York, 1971. 


H.H. Happ, J.F. Hohenstein, L.K. Kirchmayer and G.W. Stagg, "Direct 
Calculation of Transmission Loss Formula II", IEEE Trans., PAS, 


yol. 83, pp. 702-707, 1964. 


Pee PAE i 
ee 


- oaget3 mai ba ec) sonsnosin tah Iteomral adits 


<tasttes "t iE 


Wat ht; 


aor oi U1 wera ft | «toneT azeye vera oriaef 


} 


abso. tuo he: pnigense nee avid F2qsbA" .enieoh: ons ngeatt 


1 


eet 3431 eno not] veridsok bri asoonuessh bso fo bee 
A a, poh Aer, XNT-NENE aq) 


zagon | smite vemos vewod 0 nottsy9q0 atwonas* .noees 


, - asteye x9wo nt seo) voile whevevtal, aul oft ae betnse 
en nel <¥R1 epson cottage 


“not totbans a0. +09 zinatat Vad a 40 vb 
| nl sean maa od 36, 6 ono <AMSBRIDT 


38. 


oy. 


40. 


4). 


42. 


HS: 


rad 


J. Ricard, "The Determination of Optimum Operating Schedules for 
Interconnected Hydro and Thermal Stations", Revue Général de 1' 


Electricite, Paris, France, Sept. 1940, p. 167. 


W.G. Chandler, P.L. Dandeno, A.F. Glimm and L.K. Kirchmayer, 
"Short Range Economic Operation of a Combined Thermal and Hydro- 
mlectricuPowensSystem » ATER Mansa N01 4) 7 2, Pisa Ly - pp. 1057 
- 1065, 1953. 


J.d. Carey, "Short Range Load Allocational Hydro-Thermal Electric 
Systems: SAIEBydranS on VO) 2973..Pt. gLLL=B «pp. Lb05-1111,,.1954. 


R.A. Arismunandar, "General Equations for Short Range Optimization 
of a Combined Hydro-Thermal Electric System", M.Sc., Thesis, 


University of British Columbia, 1960. 


R.A. Arismunandar and F. Noakes, "General Time-Dependent Equations 
for Short Range Optimization of Hydro-Thermal Electric Systems", 


ATEESTranses PtedI Ils Voley82snpp. 88-935..1962: 


L.K. Kirchmayer and R.J. Ringlee, "Optimal Control of Thermal- 
Hydro System Operation", IFAC Proceedings, pp. 430/1 - 430/6, 1964. 


B. Bernholtz and L.V. Graham, "Hydro-Thermal Economic Scheduling, 
Part I: Solution by Incremental Dynamic Programming, AIEE Trans., 


Paper No. 60-837, 1960. 


saat FennstT-orbutt Tenolysootth bso ont 3 pore 


( ebttteaore, 0 oI 99 GET) LOW, 


not aesiniae0 ones srodl2 46% enot tsp ts 9p" 


afeon 4 92 a Scsachi ened iw 


| enot sup Inabnege-aaT Tevenaa” pri: 
ema: otrtbeia fanned -orb vo ohana 


44. 


45. 


46. 


47. 


48. 


49. 


218 


E.B. Dahlin and D.W.C. Shen, "Optimal Solution to the Hydro- 
Steam Dispatch Problem for Certain Practical Systems", IEEE Trans., 


Vol. PAS-85, No. 5, pp. 437-458, 1966. 


M.E. El-Hawary and G.S. Christensen, "Optimization of Hydro-Thermal 
Generation Schedules via Functional Analysis", Proceedings, Second 
Southeastern Symposium on System Theory, University of Florida, 


Gainesville, Florida, March 1970. 


M.E. El-Hawary and G.S. Christensen, "Application of Functional 
Analysis to Optimization of Electric Power Systems", International 


Journal of Control, 1972. 


M.E. El-Hawary and G.S. Christensen, "Functional Optimization of 
Electric Power Systems with Variable Head Hydro-Plants", Proceedings 
Third Southeastern Symposium on System Theory, Georgia Institute of 


Technology, Atlanta, Georgia, April 1971. 


M.E. El-Hawary and G.S. Christensen, "Optimum Scheduling of Power 
Systems Using Functional Analysis", Short Paper, IEEE Control Systems 


Society Transactions, 1972. 


A.M. Sasson, F. Aboyte, R. Cardnas, F. Gomez and F. Viloria, "A 
Comparison of Power System Static Optimization Techniques", IEEE 
1971 Power Industry Computer Applications Conference. To appear 


in IEEE Power Engineering Society Transactions, 1972. 


be 


-bnedse er d ural tants onit stv el sits . 
bso 0 o steaaetat irons me aya’ 90° a 
bay ies Ore dora 


ster stoi ho ‘om 


4) ve ett Tae “Z o, : Gh ie as i bile a ain Wg 1 ad bh ur 


| si 


cepa L*atnst-onbyh bso atastiev ag Fw : SFeYe NSN 
10 sdustent aioe? «vosett mot eye no sa zogmn ava3 2698 
on er Teer Dh «stoves 

13008 ‘wee a getten rer ts erewsH-T3 
ennsi | Fonaned 2a et ron . *atevten rats ae mt Fel) emad 


a es iow 


9 of ay 
PM a esd) 
alee 4 


90. 


Sis. 


O25 


53. 


54. 


05. 


56. 


Os 


219 


H.A. Burr, "Load Division Between Common-Flow Hydroelectric 


stations, MaSc., thesis. M07. AO4). 


P.R. Menon, "On Methods of Optimizing the Operation of Hydro- 


Electric Systems", Ph.D. Thesis., University of Washington, 1962. 


R.H. Miller and R.P. Thompson, "Long Range Scheduling of Power 
Production, Part Ii", Paper C72 160-5, IEEE Power Engineering 


Society, Winter Power Meeting, New York, 1972. 


M.N. Ogdztoreli, "Time Lag Control Systems', Academic Press, 1966, 


New York. 


M.E. El-Hawary and G.S. Christensen, "Functional Optimization of 
Common-Flow Hydro-Thermal Systems", IEEE Paper Engineering Society 


Transactions Paper T72 172-0, 1972. 


M.E. El-Hawary and G.S. Christensen, ‘Extensions to Functional 
Optimization of Common-Flow Hydro-Thermal Systems", Submitted to the 
IEEE Power Engineering Society's Trans. 

P.L. Falb and J.L. De Jong, "Some Successive Approximation Methods 


in Control and Oscillation Theory", Academic Press, New York, 1969. 


M.E. El-Hawary and G.S. Christensen, "Functional Optimization of a 
General Hydro-Thermal Electric Power System", Submitted to the IFAC 


Journal, Automatica. 


’ a 


Ct ee Hoe feorbul of) nasuiog, pola 
f } ii ’ att he 


NUR Ted M e2kae r 
\ he eae i) \ 
i 2 iP i y hi ; " oi ay" 


~obYH Foot 0 ise ott anfsin a0 3 
305 he we im 4 y x Be a 


die i a cee a P Ay ieee 
oORT. <soaeni asy to ate i vevinll:, e oheadT His a sPeses0 at 
* oy ‘ i : | Rs nee : pad A F q ‘i 
g at \ . en 4 i ies 
mieye | Oe ee 
"19 i i ive a ie tan aime roeqnodt 4. A, bas, ee “ 


“pai vesnf pad 1909 agSt 2.001 '§ gi9 soa “Sh dred «srokoybo 


ae | Me Sher HOY walt ¢ OH peta wedi M «tatoo 
| | | ne : 2 ar ‘ ia oe nN ah P e) 
nthe Wicd aha ( 
ae! veo? simabe> " , Br at 24 ie. lp ie 006 a6. omer" 7 r {pi08s000, Me al 


an) ae aS ivy Denies ree ho ut 


» 


+9 nottps cinta Tenoitonut dens tatyno ce bas vanaf 


Be 


ieieiov?. gatrpanton 124 464,332 “ezmosey2 ‘emradiyoriogh wort 


oop SRO1 DSSS aa atl 


Ry Bik et Loy ? \ Pin) 


fanattoni' at enotenedx3" Ay sanad revi) ee as via aH 


i are | i \! 
e \ 


at @ _ bad ime 2 dhanetene IegreeT orb wo -niamihan 


i ; *, Li uh Ale A ) ~ : haa ‘ f, 
Pe ae a an ante < de i902 —— irk | 
ae th Ue i ; ¥ i a : sic “hy * y ety 7 


/ Bear ‘iam votseibxoIgal, ay tzesooue, ‘m0, «200% atl , vf 


eg Rt” gy 


on -00eF nat wait aah ate etoad noe 


- 7 
] | , i ; , tl FC f btn ald ae 
ee AW, ae lal 1 ei | Sead | Gas 
y a) i} fl We i m2 rug aaa ani ath Aa pa 
i ; 4 ae 
f 


| Oe, soll sspanaie a ie 


220 


08. J. Carpenter, "Contribution a 1'étude du dispatching économique", 


Bull. Soc. Franc. Elect., Ser. 8, Vol. 3, August 1962. 


99. J. Peschon, D.S. Piercy, W.F. Tinney, 0.J. Treit and M. Cuenod, 
"Optimim Control of Reactive Power Flow", IEEE Trans., Vol. PAS-87, 
pp. 40-48, 1968. 


60. H.W. Dommel and W.F. Tinney, "Optimal Power Flow Solutions", IEEE 
Trans., Vol. PAS-87, pp. 1866-1876, 1968. 


61. G.W. Stagg and A.H. El-Abiad, Computer Methods in Power System 
Analysis, McGraw-Hill, New York, 1968. 


62. A.M. Sasson, "Nonlinear Programming Solutions for the Load Flow, 
Minimum Loss, and Economic Dispatching Problems", IEEE Trans., 


yol.PAS-88, pp. 399-409, 1969. 


63.  A.H. El-Abiad and F.J. Jaimes, "A Method for Optimum Scheduling of 
Power and Voltage Magnitude", IEEE Trans., Vol. PAS-88, pp. 413-422, 


N969. 


64. R.L. Sullivan and 0.1. Elgerd, "Minimally Proportioned Reactive 
Generation Control Via Automatic Tap-Changing Transformers", 
Presented at the Power Industry Computer Application Conference 


(PICA), 1969. 


ca ean Ni aac a spi Me 
aaa! Manette Wo newoe tamego “i als jW brs fannmod "WAR 
‘oe aret-ae ae tee ov | 


: iota woot it Aloria9M eduiqiod'” tal _ i 
BORE caret wei i Ate 
Hore ‘eek ails: yor aHoidutoe pm ia 
“ea ccela | ametdos4 nto | 


\ 4 1 hil 


to + eat ine ro wh bai an \ a | i 


{ . aie 
1 
o bh ; $ i. : 


hha Aa © 4 


65. 


66. 


67. 


68. 


69. 


70. 


221 


C.M. Shen and M.A. Laughton, "Determination of Optimum Power- 
System Operating Conditions Under Constraints", Proceedings 


IEE, Vol. V1. pp.7220-239, 1969. 


M. Ramamoorty and J.G. Rao, "Load Scheduling of Hydro-Thermal 
Generating Systems Using Nonlinear Programming Techniques", 


Proceedings IEE, Vol. 117; pp. 794-798,.1970. 


C.M. Shen and M.A. Laughton, "Power System Load Scheduling with 
Security Constraints Using Dual Linear Programming", Proceedings 


TEE, NO Meet 76, PDs cl teealet a 197 Uy 


M.R. Gent and J.W. Lamont, "Minimum Emission Dispatch", IEEE Trans., 


Vol. PAS-90, pp. 2650-2660, 1971. 


J. Peschon, D.W. Bree, and L.P. Hajdu, "Optimal Power-Flow Solutions 
for Power System Planning", Proceedings IEEE, Vol. 60, pp. 64-70, 
1972. 


M.E. El-Hawary and G.S. Christensen, "Optimal Hydro-Thermal Load 
Flow Using Functional Analysis", Journal of Optimization Theory and 


Applications, 1972. 


MRR 93 


: a ; m Ma " 
t te . mots 
aa i 


gas ar a . 


Past ortane 


| AN, Ve 
| a nk JHE a ae are un 
| Ny oS oe ‘s ¥ : . fe ‘ Ae " Lea i 
amines it - Gi ca | pat hit 4 be Bo ioe Ak bis ihre 
: Ct PW ed PO ated ; Sdn eee 
ATS.) ae up te oat, ante 1p ont ‘sant Foot ea 
TS ed CS 1 + | as 
: Ye a a 24 aM. Fo 4331. ae nt 
Ta | Ba | ae tae EN: ’ V9 Ay: : ie ; a a a 
-: , his " if ney er Cae a i q at a y) ee i ah 
igtw pnt tubsdee: badu meteye “aod” soapy cS “ bits, ro * | 


yy 


ibesaow4 ¢ BA Paiatg ors el iain 9 ud eA nied) ante 2n00, 


Y Saat 4 <nbisé erg” nor 2a f nd mus rie gnome 


a ree re es. hi 
Wied) Mai asthe a 1,83 00. e 


hal 
TOTP OE 2, Wot a9MO aint aqgo” sti 4. 2 it 


| | | a 
an (ote hin’ 8 te aaa, zentbonooy® 4 ‘ 


Lae 


APPENDIX A 


CONVERGENCE CONDITIONS FOR THE ALGORITHM (5.2.227) 


Consider the operator equation (5.2.226) of Chapter 5. Then using 
(5.2.217) through (5.2.220) this is rewritten as: 


V bey 
TEy(t)] =H M(tyc +f 6 M(t,s)LEly(s).s) - W(s)¥(s) Ids 
0) 


(A.1) 
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MN 2 2. 93 357 
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V 
MN, 


TyEv(t)] = fie M(t,s)Ety,s) - W(s)Iv(s)as (A.6) 


(2) coat 
ee 


‘ F 
An expression for Days) - V(s)] can be obtained from (A.3) and (A.4). 
This is given by: 


0 0 0 0 
aF A 2 eagay aa 
ay ls) 8) - V(s)] = (A.7) 
0 0 0 0 
ayy aY5 a¥,  -BYg 
Thus af 
- } g,,(t,s) =— v.(s) 
4 of. 
Vn aF pee ay5 v,(s) 
G (t.s)[5y(y(s).s) Zz V(s)]v(s) fra 4 af 4 
- } g,,(t,s) =~ v.(s) 
521 12 ay; J 
ts 
- ) goo(t,s) ——v.(s 
j=] Ze 3Y; q) 
(A.8) 
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The modified contraction mapping defined in (5.2.227) is 


Y= PY = [I-uy- [tly] - uly)] 
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This has a derivative given by: 


Py(v(t)] = (L-U]7'[Ty[w(t)] - U v(t)] 


In the case when U is given by: 


Then 
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Now the norm of Pytv(t) J is given by 
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The convergence conditions for the algorithm (A.16) are 
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Consider the components of the derivative operator Ty as given by 


(A.12) through (A.15). These yield: 
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Then (A.30) and (A.31) yield: 
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Equations (A.84) and (A.85) together with (A.46, (A.85) and (A.47) define 


the convergence condition (A.23). 
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value of Yo: If one is interested in defining the region of convergence, 
then n must be calculated. However, the crucial convergence inequality 


is (A.46). 
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